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A theory of nuclear forces is developed, based on the hypothesis that the interaction between 
a nuclear particle and the mesotron field is strong. Two types of mesotron fields are considered, 
the charged scalar and the neutral pseudoscalar. The latter, for large enough separation, gives 
forces between two nuclear particles of the same type as those obtained from perturbation 
theory and, hence, with the spin dependence and exchange properties required to fit experi- 
ment. However, at closer approach the forces become ordinary (non-spin-dependent). It is 
found impossible to obtain spin dependent forces which, at the same time, extend to small 
separations and are of sufficient strength to account for the properties of the deuteron. 


I. INTRODUCTION 


T is well known that the interaction between 
mesotrons and nuclear particles cannot be 
treated as small; such perturbation treatments 
are not only inconsistent—they lead to inac- 
ceptable nuclear forces and to a much too large 
scattering cross section for the mesotron.* Oppen- 
heimer and Schwinger' have shown that the 
strong coupling theory affords an explanation of 
the smallness of the scattering cross section ; the 
question of nuclear forces was left open by these 
authors. 

It is not difficult to see that in the limit of very 
strong coupling, the nuclear forces will not be 
right. For this there are two related reasons: 

(1) For strong enough coupling, the mesotron 
field due to the interaction of two nuclear par- 
ticles becomes large compared to the field of a 


*See W. Pauli, Abstract No. 25, Bull. Am. Phys. Soc. 
18, No. 1 (1943). Professor Pauli considers the possibility 
of removing ty terms associated with a point source 
by a subtraction formalism, while reducing the scattering 
cross section by means of radiation reaction. 

and J. Schwinger, Phys. Rev. 60, 


single mesotron ; field fluctuations become unim- 
portant, and the equations of the field may be 
treated classically. The classical solutions of the 
field equations of course exhibit the same sin- 
gularities as arise in the lowest order perturbation 
theory. (2) At the same time, the interaction 
energy becomes large compared to the excitation 
energy of proton and neutron isobars, so that it 
is energetically favorable to excite many isobaric 
states; the total spin of each particle (intrinsic 
spin or isotopic spin, plus spin or isotopic spin 
of the associated mesotron field) becomes 
effectively very large, and behaves also in an 
essentially classical manner. The spins are thus 
free to orient themselves in such a way as to 
reduce the potential energy to a minimum. It 
follows that the interaction energy is an ordinary 
potential, non-exchange and spin independent, 
with singularities as bad as those of the per- 
turbation theory. In the following sections it 
will be shown that these expectations are fully 
borne out by the detailed quantum-mechanical 
calculations. 

It is, of course, true that with fixed magnitude 
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of the coupling constant, the interaction energy 
between two nuclear particles becomes small as 
the distance between them increases; for suf- 


ficiently large separation, the interaction energy . 


becomes smaller than the energy of excitation of 
isobars, and beyond this point the forces change 
over to the type predicted by perturbation 
theory. The question that now arises, and which 
is the subject of our investigation, is whether a 
theory with an intermediate coupling strength 
can give a satisfactory account of the forces and 
whether the coupling constant can be so chosen 
that the region of spin dependent exchange forces 
is sufficiently extensive and the forces in this 
region sufficiently strong. 

As we have just seen, for a theory with spin 
dependent coupling the nuclear forces remain 
spin dependent only for separations large enough 
to make 


< (xa) /g*, 


where g is the coupling constant, a the radius of 
a nuclear particle, and (x)~! the range of the 
forces. The observed scattering cross section 
demands xa<0.1. It should be remarked that 
the condition for the validity of a perturbation 
treatment is that the self-energy of a nucleon 
due to its interaction with the field be small com- 
pared to the mesotron rest energy, g’?/(xa)*<1. 
Weak coupling would thus require g?<(xa)® 
~0.001. Forces derived from these constants are 
100 times weaker than actual nuclear forces. 
Thus we have only to consider g*>(xa)*, or 
strong coupling. 

The second decisive point in our investigation 
is this: Although for g’e~*"/xr<(xa)/g’, the 
forces given by the strong coupling theory are of 
the same type and radial dependence as those of 
perturbation theory, they are reduced by a 
numerical factor f : } for the charged scalar, } 
for neutral pseudoscalar, charged and sym- 
metrical pseudoscalar.2 These factors arise be- 
cause the spin and charge of a nucleon oscillate 
with high frequency in their own fields even in 
the absence of a second nucleon; their com- 
ponents responsible for the nuclear interaction 
have smaller expectation values. 


2 We are indebted to Professor Pauli for the results in 
the charged and symmetrical pseudoscalar cases. 
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If now g?<1, the spin dependent forces are too 
small; if g?>1, the forces are spin dependent 
only for xr>~2, and their maximum depth will 
be less than ~ 10-2uc? or 1 Mev. Thus we see that 
with xa < 0.1, no value of g’ gives spin dependent 
forces large enough to agree with experience.® 

It may be remarked that values of a of the 
order (x)~! not only conflict with the data on 
mesotron scattering, but essentially render 
nugatory a field theory of forces, since for r~a 
these forces are entirely determined by the 
nature of the source, and not by that of the field. 


II. CLASSICAL THEORY 


As has been pointed out in the introduction, 
the main features of the strong coupling theory 
of nuclear forces can already be seen in a classical 
theory. We shall, therefore, begin with an in- 
vestigation of the symmetrical scalar theory in 
the classical limit. 

The heavy particle will be supposedly spread 
over a finite region of radius a with a density 
U(r), fdrU(r) =1. In accordance with our above 
remarks concerning the size of the source we 
suppose xa<1 where x=yc/h and yu is the meso- 
tron’s rest mass. The Hamiltonian of the sym- 
metrical scalar theory when two heavy particles 
are present is 


H=} f 


f de[ (1) 


Here U, and U, are the source functions of the 
two heavy particles, ¢, and ~ their isotopic spins. 
The wave function ¢ is a vector in ¢ space; its 
x and y components are the real and imaginary 
parts of the charged field, its z component is the 
neutral field. The operator w*=(x?—A), whereas 
the coupling constant g is related to the dimen- 
sionless parameter g’ through g=g’(hc)!.4 The 


3 According to investigations of Nelson and Oppenheimer 
(unpublished) strong —— pair theories give a mesotron 
scattering cross section equal to the square of the range of 
the forces, and forces which are therefore inacceptable. 
In addition, in these theories it is impossible to combine 
spin dependence with saturation. 

We adopt rational units throughout: A =c =1. 
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equations of motion are: 


(2) 


hta= f dr¢(r)U.(r), (3) 


These equations will be treated classically ; that 
is, we shall ignore all commutators of the quan- 
tities which appear. 

For strong coupling, g>1, we can find a class 
of solutions which are non-radiating (i.e., involve 
only frequencies smaller than «x, and thus no 
wave zone field). We take 

The real and imaginary parts of the complex 
vector are the components of and 
along two axes which form, with <0, an orthogonal 
axis system. The solution of (2) is then 


$(r) f V(r, + Use") Jar’ 


f V(r, + Use’) (4) 


Where Y is the Yukawa potential : 


e —x|r—r’ 
Vie, 
|r—r’| 


and 


|r—r’| 


This expression for ¢ substituted in (3) gives 


— r’) + (5) 
or, expanding Y, in powers of v to order v?: 


= X41) (v?/2«)(R+S) 


with 


R= f Use) exp (—«lr—r'l)~1, 
(xa<1) 
S= f f exp (—«|r—1'|) ~e-*, 


d being the separation of the two sources. Hence 
we obtain a solution if 


v= 2«/g?(R+S)ro. (6) 


Solutions with +, not parallel to * give high fre- 
quencies of order 4rg* {(f'drdr’ U,(r) U(r’) Y(r, r’) 
in contrast to the quasi-static solution just found. 
The energy associated with our quasi-static 
fields can be calculated from (1) which gives, to 
order 1/g*: 
2x 


(7) 


f f drdt'Us(t) Use’) V(t, ’)~1/a, 

(xa<1), 

Judy f f drdr’U.(t) Y(t, 
(d>a). 


The dominant g’? term contains the self-energy 
of the two heavy particles, —g*J, and the inter- 
action energy between them, —g*J. The small 
term in 1/g* may be rewritten in terms of the 
charge vector Q= {'(¢X)dr, whose z compo- 
nent gives the mesotron charge in units e. Using 
(4) we find Q= —2¢,e*"', in agreement with the 
constancy of the total charge vector Q++.+*. 
The first term in (7) is thus 


kK Q? 
2g7(R+S) 1—Q? 


and represents the classical analogue of the 
heavy particle isobar energy discussed by 
Wentzel.® This non-radiating solution can be ob- 
tained only for values of the charge Q<1. 

The most striking feature of the classical cal- 
culation is that, for the solution of minimum 
energy, there is a strong coupling between the 


5G. Wentzel, Helv. Phys. Acta 13, 269 (1940). 
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directions of the isotopic spins of the heavy 
particles causing them to line up parallel to each 
other. This has as a consequence that the inter- 
action energy —g*J is an ‘“‘ordinary”’ potential, 
completely independent of isotopic spin. It will 
be observed also that for a—0 the singularity in 
the potential is of the same type as yielded by 
perturbation theoretic treatment. 

This treatment of the classical theory is due 
to Professor J. R. Oppenheimer. We are very 
grateful to him for permission to quote his 
results. 


Ill. CHARGED SCALAR THEORY 


We shall begin our consideration of quantum 
field theory with a discussion of the charged 
scalar field, the mesotron theory originally pro- 
posed by Yukawa. Although this simple theory 
cannot be expected to yield the observed nuclear 
forces, all the cogent results can already be ob- 
tained—no essentially new factors appear in the 
more complex fields which have been invoked to 
explain nuclear forces. 


A. One Nuclear Particle 


A treatment of a single heavy particle coupled 
to the charged scalar field has been given by 
Schwinger.* However, we wish to introduce a 
somewhat different mathematical treatment 
than was used by Schwinger, or by the present 
authors in their original work.’ The change 
consists in expanding the fields in terms of the 
Yukawa function w~*U(r), instead of in terms of 
the source function U(r) itself. This more con- 
venient treatment was suggested by Professor 
W. Pauli and it is with deep appreciation that 
we acknowledge his permission to use it in this 
paper. 

The Hamiltonian of the heavy particle plus 
mesotron field is 


H= f 
(8) 


The notation is the same as that used in Section 
¢ J. Schwinger, unpublished. 
7 t M. Dancoff nen R. Serber, Phys. Rev. 61, 394 (1942). 
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II except that ¢(r) is a complex mesotron field 
amplitude and x(r) its canonical conjugate. 


Cx(r), o(t’) ]= (9) 


The operators r_ and 1; change the projection — 


of the heavy particle’s isotopic spin from proton 
to neutron and from neutron to proton, respec- 
tively. We introduce also 


N=4r | dr[ U(r) }?, 


X(r) U(r), 
that is 
(x?—A)X(r) =4r U(r), 


1 
I= (10) 


&(r) =X(r)/TJ, 
so that 


f dré(r) U(r) =1. 
It is convenient to split ¢ and z into ‘‘coupled” 
and “‘uncoupled”’ components as follows: 
11 
a(r) = (42) U(r) +7,(r), 


where we require of the functions ¢;(r) and 7;(r) 


that 


f U(r) = f (12) 


and are canonically conjugate expansion 
amplitudes. They can be determined by 


(4x)! | dro(r) U(r), 
(13) 


f drx(r)é(r). 
m(r) and ¢;(r), however, are not exact canonical 


conjugates, but satisfy the commutation relation 


Cwi(r), U(r) g(r’) (14) 


| 


(10) 


ed” 


11) 


i(f) 


12) 


sion 


13) 


ical 
‘ion 


14) 


In terms of these variables 


+ f +#1(r) (r) 
+ (41) !#o f drm,(r) U(r) 


+ (41) f dr#;(r)U(r). (15) 


Introduce polar coordinates : 


pot DX) ], (16) 

], 
17) 


where 


=—{ f — #1(r)$,(r) J. 


Here (po, do) and 8) are real canonically con- 
jugate variables. The total mesotronic charge is 
given by 


ail f (18) 


The complex variables f:(r), gi(r) satisfy the 
same commutation relation as do m:(r), ¢:(r), 
namely, (14). The Hamiltonian becomes 


(pe? — 1/4) ] 


+ f + 
+(4m)¥(po/2) f 


f dr U(r) 


(19) 


po* is symbolic expression for (1/go)PogoPo. We have 
m4.’ terms of the type 2/go because, as the later work 
shows, = vanishes in the absence of free (unbound) meso- 
trons. Such terms are of higher order in 1/g and are 


negligible if g/pp>>1. 
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We transform to a representation in which the 
interaction energy is diagonal. The eigenstates 
of the heavy particle for this term are not 
proton and neutron, respectively, but are 
derived from the latter by a unitary transforma- 
tion F’=SF where S is the unitary operator 
exp (i7,0/2). When the states of the heavy par- 
ticle are represented as proton and neutron +_, 
7+, and r, are given by the matrices 


After the S transformation’ the Hamiltonian 
becomes H’ = SHS-". In particular the expression 
is transformed into (r_+ 174) = rz. 
If we choose for the states of the heavy particle 
a representation in which 1, is diagonal, the 
interaction energy is found to be diagonalized 
with eigenvalues +ggo. On the other hand, off 
diagonal terms are introduced elsewhere in H, 
since pp—>fe+7,/2 and in this representation r, 
is off diagonal. In the strong coupling limit, we 
may choose the lowest value of the interaction 
energy and neglect terms linear in 7,, as coupling 
widely separated states. The condition implied is 


g/pe>t. (20) 
The Hamiltonian becomes 
H' =(N/4)(bo? + (21) 


The term linear in go is eliminated by a shift in 
origin : 
Qo=go—gl/2. 


Since gI/2 is a constant, po is conjugate to Qo. 
We get 
H’ = 
(23) 

The leading term, —g*J/4, gives the energy, to 
order g*, of the static field bound to the heavy 
particle. ; 

It is convenient to resolve the complex com- 
ponents of the field into real and imaginary parts: 

pi(t) = (1/v2)[p.(r) +7p,(r) J, 


(24) 
qi(r) = (1/V2)[g.(r) —tgy(r) J. 


® For previous use of the S transformation, see W. Pauli, 
Helv. Phys. Acta 12, 147 (1930); also reference 12. 
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Here again (2, q:) and (py, gy) are real pseudo- 
conjugate pairs satisfying the commutation 
relations (14). The Hamiltonian separates into 
two commuting parts: 


= —g*I/4+H.+H,, 


f 


f 
+4f (25) 
f 
X [au (r)w*gy(r) + 


H, gives the energy of free oscillations about the 
equilibrium position go=gI/2. For, if we con- 
struct the functions, 


= (4) + 
= + (1/v2)p.(r). 


Then ¢, and 7, are canonically conjugate and 
H,= (26) 


H, is diagonalized in plane waves and represents 
a field of unscattered mesotrons. 

Meanwhile H, can be transformed into a 
similar form if we first eliminate the term linear 
in p,(r) by applying a unitary transformation. 


F=eUF’e-U= F’+iU, F] 


where 


f dtX(r’)qy'(r’). (28) 


This produces a canonical transformation from 
the variables (py, gy) to (py’, dy’), each pair satis- 
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fying (14). To first order in 1/g: 


pu(t) =py'(t) — U(s)/I 


—X(r) / f (29) 
Qu(t) 


It can be shown that the criterion for neglect of 
succeeding terms in (29) is, once again, g/pe>>1, 
If (29) is substituted into (25), the linear term 
is eliminated, and 


/ dr{X(r)}? 


+4 f (30) 


Here we have substituted for go the equilibrium 
value gI/2 in accordance with the strong 
coupling approximation, a procedure which can 
be proved valid for g>1. The latter two terms 
in (30) would, in analogy to H,, respesent the 
energy of a set of free mesotrons except that ,’ 
and q,’ are not exact canonical conjugates. We 
will not here calculate the eigenstates of H,; we 
will only remark that they can be represented at 
large distances as plane+scattered wave, leading 
to a scattering of free mesotrons by the heavy 
particle.” Since we are interested only in the 
heavy particle in the absence of unbound meso- 
trons, we may assume that none of the unscat- 
tered ‘‘x’’ waves nor any of the scattered “y” 
waves are present. The total energy then reduces 


to 


=— gl / g? f dr{X(r)}%. (31) 


The latter term thus appears as a correction to 
the static self-energy of order 1/g’. Since for the 
states under consideration the expectation value 
of r,=0, it follows that the expectation value of 
the charge of the heavy particle core =}. Hence, 
if Q is the total charge of the system pp=Q--}. 
For small sources, (xa>>1), 


(1/40) f de {X(r)}?2=1/2k. 


10 J. Schwinger, unpublished. 


Suc 


he 
W 
| 
+10, LU, (27) 
_ 


MESOTRON THEORY OF NUCLEAR FORCES 


Hence the isobar energy takes the form components : 
2uc?(Q—4)2/g*. (32) = *{ 04 E4 


Weare restricted by our approximation to isobars w(t) = (4m)! { U4 +207 U9 (8) } +21(8), 
for which (Q—3)/g<1. where we require 


B. Two Nuclear Particles f dr¢,(r) U4(r) = f dr,(r) U(r) =0, 
The Hamiltonian is (36) 


H= f dr[#(r)x(r) f drm,(r)t4(r) = f drm,(r)£4(r) =0. 


The amplitudes ¢o and m» are determined by 


= (4)! f dro(r) U4(r); 


+ +747 (33) 

The procedure parallels closely that for a single 0 = (4x)! f drg(r) U*(r), 
heavy particle. We introduce 

XA(r) rot = f 

X3(r) U(r); 


N=4e | dr[U4(r) | dr[U*(r) mo? = (4r) 4 f drr(r) £4(r). 


(go4, mo“) and (08, mo”) are conjugate pairs. 
M=4r However, 
[mi(r), ]= —i[6(r—r’) 
and define ~ — (38) 
In terms of these variables 
+[1/(I?—J*) +4607 G0") 
J= f drX4(r) = f drX®(t)U4(r). 
+ +179 +747 ] 


I= f drX4(r) U4(r) = f drX4(r) U(r), 


We now introduce functions orthogonal to the 
U's: + f 


drt4(r) = | dré#(r)U4(r) =0, 
f ré4(r) U9(r) f + (4x) f drU4(r)m,(r) 


dré4(r) UA(r) = f dr®(r) =1. + f dr 


Such functions are 
Ix4-—Jx8 —JX4 
~ 


(34) 


+ f (39) 


We expand ¢, x into coupled and uncoupled 
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We introduce polar coordinates for ¢o4, do? : 
$04 exp (—104); exp 


It is more convenient, however, to use instead of 
64, 05 the quantities 


py=por—po?. 
The canonical transformation takes the form: 
04 =qo4e-i@+v); 
X { (bet 
= ; 
#o4 = /2 — 
X 
; 
/2 po? + (40) 
X 
X J 


f 


- The independent conjugate pairs are (po4, go“), 
(po?, qo"), (Do, 6), and (Py. y). pi(r) and qi(r) 
commute with the other quantities, but 


= —1[6(r—r’) — U4(r) — U8(r) 


The total charge is 
f dr[x(r)o(r) —#(r)6(r) ] 


XL) +L = po. 


Before writing down H, we shall consider just 
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the leading, or interaction term. This becomes 
Be-i0-W)) J, 


Again apply an S transformation, where 

S=exp [ir.4(0+y)/2] exp 
The interaction energy becomes 

J. 

This is diagonal in a representation in which 7,4 
and 7,8 are diagonal. The lowest eigenvalue is 
—glqo4+qo®]. Off diagonal terms coupling 
widely separated states are introduced elsewhere 
in H, since under the S transformation 


Pe and 


As before, we neglect terms linear in 7,4 and 7,2 
since these quantities have no diagonal matrix 
elements for the lowest eigenstate of the inter- 
action energy. The Hamiltonian becomes (neg- 
lecting >> as before) : 
= (N/4)[(po*)?+ (bot 
+(N/4)[ (p08)? + 
cos po4 po® + (be? — py?) ] 
+(M/4) sin py)/qo” 
— po® (pet py)/go4 
XLT { (qo4)? + (go®)?} — cos 2y 


+ (4x) f drU4(t) 
X 
x f 


(pot 

x 
+(4m)¥i{ 

x f dr — 


+ f (41) 


1 The argument for the neglect of terms involving = is 
the same as in note (6). Where products of the ferm 
py f() appear, they will always be considered symmetrized. 
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1. Freesing of the Isotopic Spin 


We first concern ourselves with the elimination 
of the linear term. This requires shifts in the 
origins of go“ and 


go4 = Qo4 + (g/2)(I?—J?)/(I—J cos 
go® = Qo? + (¢/2)(I?—J*)/(I—J cos 2y). 
The part of H quadratic in go* and go? becomes 
— (g?/2)(I?—J*)/(I—J cos 2p) 
+[1/(I?— J?) JLT { (Qo*)? + (Qo*)*} 
—2JQo4Qo® cos 2p}. (43) 


The first term is the static self-energy of the two- 
particle system and the dominant part of H in 
the limit of strong coupling. It is a function of 
the separation of the two sources through its 
dependence on J. y is the angle variable whose 
canonical conjugate is py=pe,a—pes, the dif- 
ference of charges. In the limit of large g, y will 
adjust itself so as to make the energy a minimum, 
i.e., cos 2V=1, ¥=Oor In this limit, the leading 
term becomes simply 


—(g?/2)(I+J), (44) 


a result analogous to that obtained in II for the 
classical symmetrical scalar theory. For xa<1, 
the part of (44) depending on the separation, 
—g*J/2, reduces to —g*(uc*)e~**/2xd, where d 
is the separation. Because of its non-exchange 
character, this is an inadmissable potential for 
nuclear problems. 

The condition ¥=0 or y= describes a “‘freez- 
ing’’ of the two isotopic spins with respect to each 
other. It does not follow that this freezing will 
occur for all values of d because the leading 
potential term—that of order g*—falls off 
rapidly with increasing d and may eventually be 
dominated by other terms in the Hamiltonian 
whose minimizing does not require the freezing 
of the isotopic spin. That this is indeed the case 
may be seen by considering the approach of two 
heavy particles. At infinite separation they may 
be thought of simply as two independent systems 
such as described in Section IIIA. They have, 
together, an energy given by (31): 


(42) 


pn) / f 


which can also be written 
H= —g?I/2+uc?(pe? + py*)/g?. (45) 


At this stage the value of ¥ is completely arbi- 
trary. As the particles approach each other, 
corrections appear to both terms in (45). By 
comparison with (43) one sees that the g* term 
gets a correction —g*J cos 2y/2. The isobar 
energy levels are also perturbed, but the cor- 
rection is small down to separations of the order 
of x~'. For the purposes of this argument we will 
not be concerned with smaller separations and 
hence will use the second term in (45) unchanged. 
We can predict that freezing will set in when the 
separation becomes less than a certain critical 
value determined by 


x/g?= —g*J/2, 
2/g*=e-4/xd. (46) 


When the particles have approached to a separa- 
tion equal to the range, xd =1, the isotopic spins 
will still be “liquid” if g is not larger than 1.2. 
For larger values of g, freezing will take place 
at larger separations, according to (46). An 
essentially identical, but somewhat more trust- 
worthy estimate of the freezing radius will be 
given later by (55). 


or 


2. Exchange Forces at Large Separation 


For separations larger than the critical one 
defined above, this theory does indeed lead to 
exchange forces. For here we may treat the 
expression —g*J cos 2¥/2 as a perturbation to 
the zero-order energy, uc*p,y?/g?. The unperturbed 
stationary states are where gives the 
number of charge units by which the two systems 
differ. If the charge difference is zero, the per- 
turbation vanishes in first order (proton-proton 
or neutron-neutron). For one unit charge dif- 
ference, there are two independent states, sym- 
metric and anti-symmetric in the isotopic spin 
coordinate, namely cosy and siny. The cor- 
responding perturbation energies are 


Vanti-s = +2°J/4, 
—gJ/4. (47) 


This differs only by the factor } from the result 
obtained in perturbation theory. The occurrence 
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of this factor is connected with the fact, which 
will be demonstrated presently, that half of the 
field [see, for comparison, H, of (25) ] consists of 
unbound mesotrons which contribute nothing in 
order g* to the self-field of the systems. The inter- 
action (47) has numerous defects. If the anti- 
symmetry of the total wave function of the 
deuteron is invoked, then the state sin y can be 
identified with the ‘‘triplet’’ state of the deuteron, 
and cosy with the “‘singlet.’”’ Equation (47) 
would then predict the singlet to lie lower, in 
contradiction to experience. Furthermore, it 
could not explain the quadrupole moment of the 
deuteron. The strong coupling theory has the 
additional weakness that the exchange character 
of the force may disappear inside a critical freez- 
ing radius. It will be useful, for comparison with 
the pseudoscalar theory which will be treated 
later, to investigate the freezing condition more 
thoroughly and without using the approximation 
introduced above. 


3. The “Almost Frozen” System 

We consider values of g and of the separation 
which are of such a magnitude that the minimum 
of H is determined, to first order, by minimizing 
the g* term in (43). In other words, we will 
consider y to be in the neighborhood of 0 or z, 
so that cos 2) ~1—2y*. In the remainder of the 
Hamiltonian, we keep only the leading terms in 
the expansion of cos 2y, sin 2y, etc. We obtain 


H=(N/4)[(bo*)?+ 
+(N/4)[(po®)?+ (be — 
+(M/2)[popo? + (pe? — /4q04q0" 
+[1/(1?— JLT{ (Qo*)?+ (Qo*)?} 


~2FQo4Qo® (4x)'(po4/2) f drU4(r) 
x 


+ by) /4q0"] f drUA(t) 
x f dr 


+ f (48) 
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Eigenstates for oscillations of the system about 
the position of minimum energy may once again 
be obtained by resolving f; and q; into real and 
imaginary parts: 


bi(t) = (1/v2)[p.(r) J, 
qi(r) = J, 
H=—g*(I+J)/2+H.+H,, 
where 
Hz =(N/4)[(bo*)?+ (p08)? ] + (M/2) po® 
+LT{(Qo*)?+ (Qo?)?} —2JQo4Qo? ]/ (I? — J*) 


f drU4(r)p.(t) 
+ (4m) ¥(po®/v2) f drU®(t)p.(t) 


+4 f (49) 


Introduce the canonically conjugate variables 


m2(t) = (4x)4{3p04 U4 (r) 
U¥(r)} + (1/v2)p.(r), 


H,= fact { } ], 


representing a system of unscattered mesotrons. 
Meanwhile 


H, = (be — py)*/ 
+(M/8) (pe? — +2 IP(I+J)/(I-J) 


+ (pot py) f dtp, (x) UA(r) 
+(4m)}{ f drp,(t) 


+} f (50) 


The terms linear in p, are to be eliminated by a 
unitary transformation of the type (27), (29). 
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ya 
29). 


To first order, g, is unchanged and 
py(t) = (t) — 
XL(pet+ py) (U4 —cX4 (4) —dX¥(r)) 
+ —dX4(r)—cX¥(r))]. (51) 
We have replaced go4 and go® by their equilibrium 
values 
qo" 


The constants c and d [which may be deter- 
mined, e.g., by the condition that p,’(r) be ortho- 
gonal to ¢4(r) and £4(r), a condition that is 
already fulfilled by ~,(r)] are given by 

IR-—SJ JR-—SI 


where 


S= f drX4(r)X4(r). 


The condition for the convergence of the expan- 
sion (51) is po/gK1 and py/g«1. H, becomes 


H,=} f 


(I+J)? 4xpy*? (I—J)? 
R+S 80? R-S 


(52) 


The first two terms give the energy of a set of 


scattered waves, p,’ and q,’ being quasi-conjugate 
quantities satisfying (38). The isobar energy is 


2g*(R+5) S)\I4J 


(53) 
g 


Since 7,4 and 7, are diagonal, the average heavy 
particle charge is the expectation value of 
Consequently, if M 
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is the total charge, the first term can be written 
4x(M—1)? 
2(R+S)g?’ 
giving a parabolic dependence on the total charge 
of the two particle system, the minimum lying 
at M=1. 
The rest of H, describes harmonic oscillations 


of y about the “frozen” position. The zero-point 
amplitude of these oscillations is 


8J(R— 


For sources separated by d and of radius a<d 
with a defined by 


1 
a=4r drU“(r)—U4(r) 


= f (r’) exp [—«|r—r 


[r—r’| 


the various constants are given (approximately) 
by: 
I=1/a, 


J=e—“/d, 
R/4nr=1/2k, 


Since J/I<1, (54) simplifies to 


1 
(55) 
8J(R—S) 


We may take ¥=1 as the maximum amplitude 
for which the isotopic spins can be regarded as 
even approximately frozen. This determines a 
value for g corresponding to any given value 
for d. For xd=1, this gives g=1.0 asa minimum 
value below which the system will start to thaw 
and exchange forces will make an appearance. 
But since g>1 is a condition for the validity of 
the above calculations, it is clear that with the 
strong coupling hypothesis, it is not possible to 
obtain exchange forces in the charged scalar 
theory for separations as small as (x)~'. 
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IV. THE NEUTRAL PSEUDOSCALAR THEORY 


This theory, in the perturbation approxima- 
tion, gives charge independent, but spin de- 
pendent forces through a coupling of the heavy 
particle spin with the gradient of the mesotron 
wave function. We shall treat first the one-source, 
then the two-source problem. 


A. One Nuclear Particle 


The Hamiltonian is: 
H=} f 
f drU(t)é-grad ¢(r). (56) 


6 is the heavy particle’s spin vector. All quan- 
tities are defined as before with the exception of 
the following : 


0 2 
Ox; 
X(r) U(r), 
dX(r) AU 
I= far 2,3 
Ox; ax; 
dU 
§=X/I so that fan = 6;;. 
Ox; Ox; 


We split ¢, 2 as follows: 
o(r) = (4) *¢o-grad 
w(t) = (4m) !xo-grad U(r) +7:(r) 


with auxiliary conditions 


f grad U(r) =0, 


f drm(r) grad é(r) =0. 
It follows that 


$o= (4x)! f drg(r) grad U(r); 


f der(r) grad é(r). (58) 
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(0:, o;) are canonical pairs, whereas 


Cai(r), o:(t’) ]= 


—grad U(r) -grad £(r’)), 
We obtain 


$0 
K 
+ (41) f drm,(r) grad U(r) 


+4 f (59) 


The cross term in xom:(r) may be eliminated by 
a unitary transformation of the type (27), with 


U=+(4n) f grad &(r)/R, (60) 
with 


R= f 7}, i=1, 2,3 


mi(t) =21'(r) — (4x) !x0-[grad U(r) — grad ¢(r)/R] 
= ¢1'(r). 
Substitution gives 


Ha 


+} f (61) 


It is to be observed that (60) also generates a 
change in ¢o which is not taken account of here. 


A more detailed analysis, parallel to that carried 
out for the scalar field, shows that this correction 
is negligible in the limit of strong coupling. The 
effect of the additional terms which would appear 
in H would be to spoil the exact separation of 
(61) into a set (21’, $1’) of unbound waves and a 
set (x0, $0) of bound states. The condition govern- 
ing the validity of this approximation is 
(xa/g)L<“1, where L is the total angular mo- 
mentum of the bound mesotron field defined 


below.” 


A more complete discussion of the strong coupling 
condition in this case is to be found in W. Pauli and S. M. 
Dancoff, Phys. Rev. 62, 85 (1942). 


(r’)]. 


We assume that none of the (7’, ¢:’) states 
are excited, and proceed to treat the system 


defined by 
(62) 


It is convenient to introduce an orthogonal 
system of unit vectors, m1, M2, M3, defined as 
follows: Ns is parallel to ¢9 and n, has the direc- 
tion of the intersection with the xy plane of the 
plane perpendicular to n3. The co-latitude and 
azimuth of n; are a and 8, respectively. The com- 
ponents of the unit vectors are 

n,: (—sin 8; cos 8; 0), 

ny: (—cos a cos 8; —cos a@ sin B; sin a), 

n;: (sin @ cos 8; sin a sin 8; cos a). 
It follows that 


$o=QoNs, 


63 
mo = PoNs / Qo Sin a. 


Canonical pairs are (Po, go), (Pa, a), and (pz, 8). 
The angular momentum of the mesotron field is 


f 


_ where @ is the angular momentum operator. For 


our case, this becomes 
L=[¢0X 0] 

= Palit pens/sin a. 
We can also write 
mo = Ponts + (L 
The energy becomes 
H= (44/2R)[po?+ (L?—1)/qo*] 

(64) 


where L?= p,?+ p,?/sin? a. 

We transform to a representation in which the 
interaction energy is diagonal. The corresponding 
unitary operator is S=exp [teya/2 [ie,8/ 2]. 
The transformation F’=SF is applied to the 
wave function and H’=SHS“ to the Hamil- 


MESOTRON THEORY OF NUCLEAR FORCES 


155 


tonian. In particular, is transformed into 
By choosing for the states of the heavy particle 
a representation in which o, is diagonal, the 
interaction energy is diagonalized with eigen- 
values +(g/x)go. The angular momentum be- 
comes 


L--L—- (n2/ sin a) 
X(e,cosa—o,sina)], (65) 


and its square is transformed into 


pg?/sin? a+} 
+1/4 sin? a—pgo, cos a/sin?a. (66) 


Terms linear in o, and o, have been neglected 
as they couple widely separated states. Aside 
from the additive constant, }, (66) is simply the 
square of the total angular momentum of a 
completely symmetric rigid rotator, the angular 
momentum corresponding to the third Euler 
angle being o,/2. The eigenvalues of (66) are 
therefore j(j+1)—%2, where j is the quantum 
number for the total angular momentum of the 
heavy particle plus mesotron field. 

We choose the lowest eigenvalue of the inter- 
action energy, and get 


4a J?—3 
H pr 2I—gqo/x. (67) 


The linear term is eliminated by the shift 
Qo=qo—gl/kx, 
with the result: 


H= 
(68) 


The leading term in (68) is the static self-energy. 
The two succeeding terms may be described as 
the energy of oscillation of the variable go about 
its position of equilibrium, g//x. We treat the 
case where this degree of freedom is unexcited, 
i.e., when there are no positive energy mesotrons 
in the field. Hence, in the last term of (68) we 
substitute for go its equilibrium value, and obtain 
for the isobar energy 


Hy = +1) — 2). 


d by 
with 
(60) 
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For xa<1, We find 


IX(r4) — 


PR= f d ~4n/3a, = P-J? 
Buc? (69) IX(r8) — J;X(r4) 


in conventional units. The lowest state of this The fields are split as follows: 
system is identified as the proton-neutron (j=). 


a 
B. Two Nuclear Particles o(r) = oi 


The Hamiltonian is the following: 


H=} f 


f grad +-no®-grad U(r®)| 


with the auxiliary conditions: 


+ ¢(r)]. (70) 


face f =0 
. Ox i Ox i , 


Introduce 
X(r4) U(r4), 2,3 
= (r8 
f = f =0. 
Ox; Ox; The amplitudes ¢9 and x» are determined by 
aU(r4) aU(r®) 


A A B B 
1= fa ) dU(r = far ) dU (r*) 


Ox; Ox i 


B 0 
(4x)! f 


Ox; Ox A 


i= = T i= A 
ax; ax; ax; ax, roi= (4) J 
Note that all integrals vanish which involve B a 
derivatives with respect to two different x,’s. moi=(40)> f 


States orthogonal to the U’s are defined as 


follows: $03] = $05] = —16;;, 
0&,(r4) dU(r?) 0&,(r8) dU(r4) a 
fa Ox; Ox; fa Ox; Ox; 0, [m (r), ] -| d(r—r’) ax, U(r4) 


71) 


In terms of these variables 


+h 2 i 


+ (49) f drm(r) grad U(r4) 
+ (4) - f drm(r) grad U(r®) 


+ (g/x)(64- $o4+67- $0"). (72) 


The cross terms involving 7 are eliminated by a 
unitary transformation like (27), (60), but with 


A 0 
U=(4n)! E f Rié(t4) 
— S,é(r*)} +0: f x 


— S;é;(r4)} / (R?—S,?), (73) 


where 


Sim f 
Ox; Ox; 


It follows that 


Rié(r4) — 
= 


— 


U(r8) 
Ox; R? S; 


=¢1'(r). 
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We have, as a result 


A Ri 
H=4r>, 


/ { (box)? 


I 
+ (g/x) +63 bo”) 
+4 f (74) 


As in the one source case, changes generated by 
(73) in ¢o4 and ¢o” are neglected here. Again, the 
(31’, ¢1’) system of states is taken empty so that 
its contribution to (74) vanishes. We introduce 
two systems of unit vectors and 
defined as above through two sets 
of angles a4, 84 and a®, 8? 

$04 =qo4ns4," 

= poAns4 / sin a4, 

$0" = 

xo? = / qu? sin a®. 


Using 
a4, 


we find for the total angular momentum of the 
mesotron field : 


L=L4+L?. 


We now obtain the minimum eigenvalue of the 
interaction energy by applying an S trans- 
formation, with 


S=exp [io,4a4/2]] exp [io,484/2] 
Xexp [io,%a?/2] exp 
The interaction energy becomes 
(75) 


18 Symmetrization with respect to a and functions of a 
is to be understood throughout. 
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We choose a representation in which both o,4 
and o,* are diagonal, and consider the state in 
which each has the eigenvalue —1. The effect of 
the S transformation on L4, L? and their squares 
is indicated by (65), (66). For example, (L4)? 
becomes (J4)?+3 where J4 is the total angular 
momentum of source A plus its associated field. 

So far the orientation of the space axes is quite 
arbitrary and we may assume the 2 axis to lie in 
the direction of the line separating the two 
sources. Then we may write R,=R,=R,, 
S,=S,=S,, Jz=J,=J,; in general J, xJ,, etc. 

The separation of. the static self-energy is 
accomplished by the elimination of the terms 
linear in the g’s. The part of the Hamiltonian 
involved is 


KS 


It is not hard to show that the absolute minimum 
of (76) is attained when the vectors ¢o4 and ¢o? 
coincide in magnitude and in the magnitude 


of their components. If J;>0, then bes, 
while if J;<0, then ¢o:=—¢o: In any case, 
= qo® = go. Equation (76) may be written 

(0:)2 
I+|Ji| 
The direction of ¢o94 which minimizes (77) is 
determined by the relative magnitudes of J, and 
J,. In our case (xa<1), J; is negative and has 
a greater magnitude than J,. Consequently, both 
¢o4 and —¢» (i.e., n34 and —n;*) point in the 
positive z direction for a minimum of (77). The 
minimizing value of go is 

go= —(g/x)(I+|J.|), 

and (77) becomes 


— 


— 2gqo/x. (77) 


(78) 


1. Forces at Small Separation 


If the sources are sufficiently close together, 
the minimum of the Hamiltonian is essentially 
given by the minimum of (76), namely (78). 
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Both ¢o4 and ¢o” are “frozen” and are equal to 


= 
$04 = 


Equation (78) contains, in addition to the sum 
of the self-energies of two separated sources, the 
space dependent part —g*|J,|/«? which in this 
theory is quite singular, varying as 1/d* for small 
separations. Equation (78) is moreover an 
“ordinary” force, showing neither spin de- 
pendence nor exchange properties. We must, 
therefore, investigate the possibility of thawing 
of the spins due to the first term of (74). We 
expand the Hamiltonian about the frozen posi- 
tion ; (74) becomes 


— pot +44 (J4)? 
+(J*)?}Ri/2—J4 / (Ret 


+L{(Qo4)?+ (Qo®)?}I/2 
—Qo4Qo*| 
J,| +J,? 


+| ( 


— (g?/x*)(I+|Js|). (79) 


Here we have used 
pat +024 (pp4+})/sin a4, 
(J4)? = (pa*)? + a4, 
sin a4~a4; sin 
cos a4~1—(a4)?/2; cos a?~ —1+(a?)?/2, 
=go+QoA* 


with analogous definitions for source B. Qo4, the 
displacement from the equilibrium amplitude, 
is canonically conjugate to po4. We may separate 
off the first and third lines of (79) which we call 


(79) 


H, and bring to principal axes by means of the 
canonical transformation : 


pt=(1/V2)(poA + po"), 
gt = (1/v2)(Qe4 


p-=(1/V2)(bo4 — po®), 
= 
A 
(81) 


2(I—|Je|) 


As in the one-source case, all the states of these 
two oscillator systems are taken empty. 

Lines 2 and 4 of (79) give the energy of excita- 
tion of the bound system above its minimum 
position, which we shall call H,. First we sub- 
stitute. 


B=,’ 


for both sources. Then introduce Cartesian coor- 
dinates 
x=sin a cos B~a cos B, 


y=sin a sin B~a sin 8, 
+ (pz)?+ (p,?)?} Ri/2 


— Si] / (Ru? — 
+ { (x4)?+ (y4)? 


I|J.|+Js? 


(yB)2} 


— {xAxB J. (7?—Js*) (82) 


aside from a constant. This consists of two iden- 
tical systems of which it will suffice to consider 
either: we shall treat that involving x4, x*. 
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Normal coordinates are 


v= (x4-+x%)/v2, pr=(p4+p*)/v2, 


(83) 
(x4 —x)/v2, pu= (p4 — p®)/v2, 
4rp;? 
2go*(Ri+S:) — Si) 


20+ 


The integrals of interest are tabulated below for 
the limit «a<1i and also aX<d where d is the 


separation 
I~1/a', 


enw 
«/d), 

enw 
Ri =R,=4na'/3. 


S, and S, are of order xa*. In this limit J/J 
and S/R may be neglected. Consider the (rr, vrr) 
oscillator. The amplitude of zero-point oscil- 
lations is 


RE pot 


Thawing takes place only when ¥>1, 
3. 3\ 

/ ( 1+—+— ) <--, 

kd 3 «a 


or, with g~1, xa~0.1, when «r>3. 


2. Forces at Large Separation 


This problem is handled by a perturbation 
method. In zero order, one takes the sum of the 
energies of two infinitely separated sources [see 


(68) J. 
+4x[ (85) 


where 


(J4)*= (b4)?/sin? 


+1/4 sin? pp cos a4/sin? a4. 
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In this order, the angles a4, 64, a, 88 are com- 
pletely unrestricted. As the sources approach, the 
leading modification of Ho is the addition of the 
space dependent part of the interaction energy, 
which is obtainable from (74), and is 


H,={—J, cos a4 cos a® 
— J, sin a4 sin a® cos (84—B*)}g?/x?. (86) 


Other modifications of (85) contained in (74) are 
less important as long as the separation is appre- 
ciably greater than a. 

We require first the eigenfunctions of (85) or, 
more simply, the eigenfunctions of (J4)? and 
(J®)?. (J)? may be rewritten: 


(J)? = pa? + pp? + (1/sin? a)(pp cos a+})*. (87) 


We assume a solution of the form A(a)e‘. 
Single-valuedness of the wave function in 8 
requires that m be a half integer, since the 
original wave function is related to the eigen- 
function here obtained through the S trans- 
formation: F=S-'F’, with 


S=exp [toya/2] exp [ic.8/2]. 


Since o, has the eigenvalue —1, the above 
requirement on m follows. The equation for A" is 


cosadA 1 


da® sina da sin?a 


X(n cos a+})?A+(E—n*)A=0, (88) 


where E is the eigenvalue of (87). This can be 
transformed into the hypergeometric equation by 
the introduction of a new variable 


z=(1+cos a)/2, 
and a new dependent variable 


dZ 


“R. deL. Kronig and I. I. Rabi, Phys. Rev. 29, 262 


(1927). 
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where 


ao= | n| +3+(E+})}, 
Bo= 
y= |n—3|+1. 


Z is the hypergeometric function, F(ao, Bo, y, 3), 
and the solution for the eigenfunction of (87) is 


einB{ 43 Bo, 2). (90) 


For this function to be regular at a=0, 7, the 
hypergeometric series must terminate, i.e., £, 
must be a negative integer or zero. This means 
that E is restricted to the values (|m| +) 
X(|n|+k+1), where k is a positive integer or 
zero. Each value of || +2 corresponds to an 
isobaric level of the system, each level being 
multiply degenerate, according to the number of 
combinations of m and k which satisfy the above 
relation, namely 27+1. 

We are concerned with the attraction between 
two systems which approach each other in their 
lowest isobaric levels, i.e., with k=0, j=}, 
n=+4. This level is a doublet whose wave 
functions may be written 


(sin a/2)e‘*/?, 
(91) 


(cos 


These two states, corresponding to n= +4, are 
to be thought of as representing the two pro- 
jections of the “empirical” spin of the heavy 
particle (to be distinguished from ¢, the spin 
vector introduced above for the heavy particle 
core). 

We now calculate the expectation value of 
(86). The term involving J, has a finite matrix 
element only if the projections m4 and n® both 
remain unchanged. Its expectation value is 
+(—g’J,/9x*), corresponding to n4==+n*. It 
can be written where 8(Sz, Sy, 5s) 
is the spin projection operator for the doublet 
under consideration. On the other hand, the 
term in J, contributes only if the projection is 
increased by one unit for system and decreased 
by one unit for the other. The expectation value 
is (—g*J,/9x*)(s.4sz2+s,4s,8). The energy of 
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Introduce the unit vector f along the line joining 
the two systems (in the z direction). Then 


9 J. 


The integrals J, and J, are given by (84). Letting 
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xd =x, we have 


Vi (1/3)s4-s®+ {(s4-r)(s*-r) 


3 3 
(92) 
x? x 


Aside from the factor }, this is precisely the 
interaction between two heavy particles derived 
from the same Hamiltonian by perturbation 
methods. 

Acknowledgments are due Professor W. Pauli 
and Professor J. R. Oppenheimer for illuminating 
discussions on all parts of this work. 


PHYSICAL REVIEW 


VOLUME 63, NUMBERS 5 AND 6 


Cloud-Chamber and Counter Studies of Cosmic Rays Underground 


MARCH 1 AND 15, 1943 


C. Witson AND DonaLp J. HuGHEs 
Ryerson Laboratory, The University of Chicago, Chicago, Illinois 


(Received January 20, 1943) 


A counter controlled cloud chamber and two counter coincidence sets were used to study 
the nature of the cosmic rays observable underground. The experiments were performed in 


a copper mine at depths of 71, 141, 582, and 657 meters water equivalent. The data are easily 
interpreted, if one assumes that underground the primary rays are mesotrons and that the 


INTRODUCTION 


T the time of the cosmic-ray symposium! 
held at Chicago in June, 1939, it was 
evident that there was as yet much unknown 
about the nature of cosmic rays underground as 
well as considerable disagreement concerning the 
known material. The shape of the total intensity 
vs. depth curve seemed to be quite well estab- 
lished. When plotted on a log log scale, the data 
fall on a line composed of two straight portions, 
the change in slope occurring at about 250 to 
400 meters water equivalent. However, con- 
cerning the nature of the rays responsible for 
the two parts of the curve, there was a diversity 
of opinion. Several observers had found evidence 
for the non-ionizing character of the rays which 
carry the energy down to great depths ; Wataghin 


1V. C. Wilson, Rev. Mod. Phys. 11, 230 (1939). 


soft rays and showers are electronic secondaries produced by the penetrating mesotrons. 


and Santos? at 250 and 400 m,* Barnothy and 
Forro* at 1000 m. However, Wilson’s® experi- 
ments at 30 m and 300 m indicated that the 
penetrating rays responsible for the observed 
intensity at both these depths are ionizing. 
Clay’s® interpretation was that protons are 
predominant below the break in the intensity 
curve. 

Neither data nor interpretation were clear-cut 
on the matter of shower production and the 
abundance of soft particles. The ratio of soft to 
hard components, as measured by counter 
absorption experiments, varied greatly, and it 
was quite evident that much of this disagreement 

2G. Wataghin and M. Damy de Souza Santos, Ann. 
Acad. Brasil. Sci. 11 (March 11, 1939). 

*m means meters water equivalent of rock calculated 
on a density basis. . 

‘J. Barnothy and M. Forro, Phys. Rev. 55, 870 (1939). 


V. C. Wilson, Phys. Rev. 55, 6 (1939). 
6 J. Clay, Rev. Mod. Phys. 11, 128 (1939). 
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was caused by differences in geometry of 
counters, absorbers, and surrounding material. 
Thus, Wilson had found a 7 percent soft compo- 
nent, Nielsen and Morgan 24 percent, and Auger 
9 percent, all at 30 m depth. Both Wilson and 
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Fic. 1. Cloud chamber and Geiger-Mueller tube arrange- 
ment to study the nature of cosmic rays at 71 m. 


Clay found that the ratio of soft to hard rays 
increases with depth at moderate depths, and 
Clay observed a marked decrease below 400 m. 
Auger’s data showed a constant decrease with 
depth. According to these observers, shower 
intensities vary with depth in about the same 
way as the soft component although it is difficult 
to distinguish between variation in shower 
intensity and possible variation in the nature of 
the showers themselves, such as size, angular 
spread, etc. 

The question of the nature and production of 
the showers is intimately related to the study of 
the nature of the penetrating rays at different 
depths. The experiments of Auger, Wilson, 
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Janossy, Nielsen, and Morgan indicated that 
the showers at moderate depth are ordinary 
electronic multiplicative ones, accompanying the 
ionizing primaries which produce them. Barnothy 
and Forro felt that the great bulk of the ioniza- 
tion at depths below the break in the intensity 
curve is made up of soft particles, secondaries of 
a penetrating non-ionizing radiation. Wataghin 
also found a production by non-ionizing primaries 
at both moderate and great depths, it being the 
only method of production at the great depths. 
The interpretation of Clay, however, was that 
at great depths, the secondary showers, few in 
number, are progeny of the protons which are 
the particles penetrating to such depths. 

Thus, it was realized that much more experi- 
mental work was needed to help clear up this 
diversity of fact and interpretation. Clear-cut 
data were needed to show the kinds of rays 
present at a particular location and the processes 
taking place there, and to show the way in which 
these things vary as the cosmic-ray energy is 


. carried downward to greater depths. It was 


decided, therefore, to perform a series of experi- 
ments which might aid in the resolution of some 
of the difficulties mentioned. The work was 
carried on in the fall and winter of 1940-41, at 
the Isle Royale copper mine near Houghton, 
Michigan. Although it was necessary to termi- 
nate the experiments at the end of December, 
enough data were obtained in the time available 
to give results of definite value. 


CLOUD-CHAMBER EXPERIMENTS 


Since the interpretation of counter tube 
experiments underground is often difficult and 
uncertain, it seemed desirable to supplement 
them by the use of a cloud chamber. The 
immediate phenomena recorded in the chamber 
are more obvious than interpretation of counter 
readings, but quantitative results are, of course, 
much more uncertain. It was felt, however, that 
cloud-chamber results combined with counter 
data would give a good picture of the complex 
of rays present underground. The chamber would 
show directly what was being measured quanti- 
tatively with the counters; the latter would give 
some idea of the frequency of occurrence of the 
processes seen in the former. 
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\. the “1st level,” 


The chamber was of the conventional dia- 
phragm type, 20 cm in diameter, automatic 
in operation, and counter-controlled. It was 
mounted between the poles of a 600-lb. Alnico 
permanent magnet, which produced a field of 
1250 gauss in the plane of the chamber. Illumi- 
nation was provided by two 36-volt coil-coil line 
filament bi-post lamps, which, when flashed in 
parallel on 110 volts and focused with cylindrical 
lenses, produced a brilliant sheet of light of 
chamber depth. It was found possible to let the 
apparatus run continuously in the mine. It was 
visited daily to check adjustments, remove film, 
etc.” 

The cosmic-ray intensity in the mine, and 
hence, rate of expansion of the chamber, was so 


low that in the time available it was felt best to 


get as many pictures as possible at one depth, 
71 m, and use the counter data 
alone to infer the variation of phenomena with 
depth. The counter control and lead absorbers 
used to obtain a general picture of the rays and 
_ Processes present at this depth are shown in 
” Fig. 1. Three general arrangements were used : 


I. Coincidences of tubes 1, 2, 3, and 4 con- 
trolling the expansions, no Pb absorber 
present. 

II. Same counter control as I but with 10 cm 
Pb at B. 

III. Coincidences of tubes 1, 2, 3, and 4 con- 
trolling expansions ; 3 to 10 cm Pb at A+1.2 
cm on chamber; tubes X, in parallel, light 
a neon lamp when they discharge coinci- 
dently with 1, 2, 3, and 4. 


The purpose of I is to give a cross section of the 
rays in equilibrium with the rock at this depth 
undistorted by Pb absorbers. Arrangement II 


_ shows the effect of the usual method of removing 


the soft component in a counter absorption 
experiment—the lead under the chamber allowed 
only penetrating rays to be observed without 
the addition of rays produced in the absorber 
itself. In III the more complicated picture of 
penetrating rays plus secondaries produced in 
the lead above ‘the chamber is observed, the 


*The chamber resembled closely that described b 
H. Jones and D. so Rev. Sci. Inst. 11, 79 (1940). 
The magnet was an enlarged version of the one described b 
ag Phys. Rev. 59, 117 (1941). We wish to tha 
Dr. Herzog also for the use of an electronic timing circuit. 
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‘in counting rate with 10 cm Pb. The last column 
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neon lamp indicating whether secondary produc- 
tion was by ionizing or non-ionizing primaries. _ 
For these three groups, with an expansion rate 
of about two per hour, a total of 1248 pictures 
was obtained. A shower arrangement of tubes 
was also tried, but the counting rate was pro- 
hibitively low, power line failures causing more 
expansions than showers. 

For the single tracks, visible in both of the 
stereoscopic views, passing through the control 
counters, and long enough for curvature meas- 
urement, the pe value (300 Hp) was determined 
by use of the micrometer device described by 
Jones and Hughes.’ There were 364 such tracks, 
divided among the three groups as shown in 
Table I with momentum spectra as in Fig. 2. 
The momentum range shown does not extend 
above 10° ev. The tracks with higher measured 
values are classified with the ‘‘straight”’ tracks. 

The table and momentum spectra show that 
for group I (no Pb absorber), there is an excess 
of negative particles of low energy. These are, 
of course, secondaries and are undoubtedly 
knock-on electrons separated sufficiently from 
their primaries to appear singly. Aside from this 
explicable excess, there seems to be a positive- 
negative equality. A proton with an energy low 
enough to be measured in this apparatus would 


TaBLeE I. Rays photographed with the three counter 
arrangements at 71 m depth. 


Straight —10° 
Group + — (>10*% ev) Total total 
I. No Pb 22 35 117 174 33% 
II. 10 cm Pb below 
chamber 5 3 23 31 26% 
III. 3-10 cm Pb above 
chamber 18 23 118 159 26% 


be quickly absorbed by its great ionization loss. 
Therefore, this is not evidence against the 


presence of positive protons in the penetrating ~~ 


group. Incidentally, except for a few contamina- 
tion alphas, no particles with heavier ionization 
than electronic were observed in any of these 
pictures. 

The expansion rate of the chamber was 2.10 
per hour in I and 1.98 in II, a 5 percent decrease 
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of Table I, which is a rough measure of the 
“shape” of the spectrum, indicates that the rays 
absorbed are of low energy. Comparison of the 
momentum spectra of I and II (Fig. 2) shows 
this even more clearly, for the low energy group 
of I does not appear in II. There are now no 
rays of momentum less than 2.1108 ev, which 
checks well with the fact that the momentum 
value for mesotrons of range 10 cm Pb is 2.3 X 108 
ev. The eight rays still present in II with 
measurable momentum are almost certainly 
mesotrons; electrons of such momentum could 
never penetrate the absorber because of radiation 
loss; nor could protons because of the increased 
ionization loss at these momentum values (a 
proton of momentum 8.2X10* ev has a 10 cm 
range). The uncertainty in curvature measure- 
ment for the higher momentum particles in this 
group would perhaps permit their interpretation 
as protons of actually higher momentum ; how- 
ever, of the three over 5X10* ev, two are 
negative. Of the five most certain mesotrons, of 
momentum under 5X10*, four are positive. 
Protons of such momentum value also would 
show increased ionization, which is not the case. 
Those secondary rays which travel close to 


off fh mom 


Fic. 2. Momentum spectra taken from cloud-chamber 
photographs at 71 m. 


their penetrating primaries are not revealed by 
a decreased counting rate when Pb is interposed 
in a counter experiment. It follows that for this 
particular geometrical arrangement (roof 80 cm 
above upper counter) 5 percent of the total are 
soft rays separated enough from their primaries 
_ to set off a counter train missed by the primary 
rays. The actual soft component is of course 
greater than 5 percent by an amount which will 
be larger the greater the coherence of the 
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secondaries to the primaries. A rough estimate 
can be made by an actual count of the secondary 
rays observable in the chamber accompanying 
the penetrating rays of group II—the resulting 
figure is 13 percent, which should be increased 
somewhat to allow for those secondaries which 
pass out of the illuminated area. Thus, in a rough 
way, we see that an actual 13 percent soft 
component would be measured as 5 percent—or 
even less if the counter train were nearer the 
overlying rock. This point is, of course, closely 
related to the counter absorption experiments 
described below. 

Group III concerns the more complicated case 
in which the Pb block absorbs some particles 
but also creates others, which are observed in 
the chamber. In this sense, the phenomena 
observed are a combination of those in I and II, 
and it is seen in the momentum spectra and 
Table I that the results lie between those of the 
first two groups. Under the Pb block there is 
an approximate positive-negative equality—as 
would be expected from multiplication of those 
soft rays in the absorber which are not absorbed. 
Of more interest in this group than the momen- 
tum spectra are the nature and mode of pro- 
duction of the showers observed. 

The secondaries and showers observed with 
the three arrangements are shown in Table II. 
Again, the relative number of coincidences 
caused by showers to those due to single rays is 
not the ratio of occurrence of these phenomena, 
but depends on the geometry and the surround- 
ing material. Still another reason is seen here for 
the fact that the decrease in counting rate with 
insertion of 10 cm Pb gives a low value for the 
soft component, for an appreciable number of 
showers present in I are removed by the Pb, 
each of which would be counted as only one 
soft ray by the counter absorption method. The 
difference in size distribution of showers for the 
three arrangements is quite marked. The Pb- 
produced showers in III seem to contain a 
greater number of larger size than the rock- 
produced ones in I; this would be expected from 
the difference in material but, of course, can be 
due simply to the proximity of the producer in 
III. The large production of showers in a few cm 
of Pb as shown in III and their complete absorp- 
tion by the Pb in II indicate their nature as 
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¢g ordinary multiplicative showers—at least, no 
evidence is present of an unusual shower type. 
The additional, admittedly subjective, but im- 
portant fact can be added that the showers look 
just like all those photographed at sea level with 
the same apparatus. 

~ The neon lamp data in group III give definite 
information on the mode of production of the 
ionizing rays created in the Pb, which set off the 
counters below the Pb and are seen in the 


TaBLeE II. Secondaries and showers accompanying the 
counter control rays. 


Shower size 
6-10 11-20 21-30 31-40 >40 


Single 
second-| 2-5 
Group Rays aries | rays 


1 174 29%|15 4 2 1 O 
31 13%/11 0 0 0 0 
159 11%] 11 2 1 #2 2 2420 


chamber. A lamp flash on the film signifies that 
coincident with the expansion an ionizing ray 
has passed through the counters X above the 
chamber. On 11 of 251 coincidences, there were 
no neon flashes. Without the cloud chamber 
evidence, we might conclude from this that under 
the Pb 43 percent of the ionizing rays were 
produced by non-ionizing primaries; further- 
more, the primaries were unaccompanied by 
ionizing rays, hence were not photons—an 
important finding. But there would always be 
the question in such a result as to whether 
there might be an error in interpretation. 
Perhaps side showers, scattering, or insufficient 
coverage by counters X could explain the 
apparent production by non-ionizing penetrating 
rays. Here, however, we have the additional 
evidence of the chamber, and, seeing the actual 
ionizing rays in the chamber, we can project 
them upwards to see if their ionizing primaries 
would have fired the counters X. When this is 
done, it is found that in none of the 11 pictures 
showing no neon flashes is there a ray visible 
whose projection passes through X. In 7 of these 
pictures no rays are visible (a ray can miss the 
chamber and fire the counters) ; inclined showers 
are present in 2 of them; and the final 2 show 
single rays inclined so as to hit the ends of the 
control counters but not counters X. So, of the 
162 rays photographed with arrangement III, 
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there is no evidence of production by a non- _ 
ionizing agent which is unaccompanied by 
ionizing rays. Photons are not ruled out by this 
evidence, but production by any penetrating 
non-ionizing ray (which would not be accom- 
panied by numerous secondary ionizing rays), 
as suggested by several experimenters, must be » 
extremely rare. 

A number of knock-on electrons created in the 
gas of the chamber were observed. The great 
majority of these have ranges lying entirely in 
the chamber; hence, their energies can be 
estimated even when scattering prevents curva- 
ture measurement. The probability of ejection 
of a knock-on of energy E electron volts is 
given by® 

z mcdE 
p(E)dE =0.3— —- —dx, 
A E 


provided £E is small as in the cases considered 
here. The electronic rest energy is mc* and dx is 
measured in g per cm?. If the momentum of the 
primary particles is high enough so that 6=1, 
then this probability is independent of the 
momentum and mass of the primaries. As all 
the rays recorded were of electronic ionization 
density (hence with B=1), a study of the 
distribution in energy of the knock-ons in the 
gas will give no information on the mass of the 
primary particles. Nevertheless, it is interesting 
to check the validity of the above formula by 
counting the observed total number of knock-on 
electrons. For 390 rays there are 48 knock-on 
electrons with range greater than 0.5 cm in the 
chamber, i.e., of energy greater than 18 kev. 
Integrating the formula with respect to E 
between the limits 18 kev and infinity, we find 
3.89 per g/cm? of argon or per 516 cm track 
length in the chamber. Thus, for 390 tracks of 
average observable length of 18 cm, or 7030 cm 
track length, we would expect 7030/516X3.89 
or 53 knock-on electrons as compared to the 48 
observed. The number of cases is too small to 
allow any significant comparison between the 
observed and expected distributions in energy. 
Contrasted to statistical data on their energy 
distribution, cases of individual knock-on elec- 


*B. Rossi and K. Greisen, Rev. Mod. Phys. 13, 240 
(1941). 
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|_ can be obtained from the counter data. 


trons, where angles and momenta can be 
measured, furnish excellent methods of mass 
determination. In order to estimate the mass of 
the primary with any accuracy, the collision 
must be almost head-on and the primary 
momentum must be relatively low, as can be 
seen from the formula 


2mc*(uo/k)? cos? 6 
{[(uo/k)?+1 ]#+1/k}*— (uo/k)? cos? 


which gives knock-on energy E (electron volts) 
in terms of the angle of ejection @ of the knock-on 
electron by the primary, of mass km, of momen- 
tum yo mc. For the knock-on electrons observed, 
6, 4o, and E could be measured at least approxi- 
mately, but in most of these cases the values were 
such (@ large, uo large) that k could not be 
determined with any accuracy. The two knock- 
on electrons shown in Fig. 3 illustrate this. 


'~ All in all, the cloud-chamber data show that 


at 71 meters depth the kinds of rays and inter- 
actions which are observed are just those well- 
known at sea level. There are penetrating 
ionizing rays, soft secondaries (with an excess of 
negatives), and ordinary multiplicative showers. 
The penetrating rays of high energy may be 
mesotrons or protons, there being definite evi- 
dence only for the presence of mesotrons. A soft 
component, composed of low energy rays, is 
present. The fraction of soft rays is greater 
than would be indicated by counting rate ratios 
in an absorption experiment. The soft rays 
apparently are largely produced by well-known 
knock-on processes; the number of electrons 
ejected in the chamber gas checking with 
theoretical knock-on probability, and the greater 
number of absorbable secondary rays being 
negative. Under a shower producer of moderate 
thickness the number and size of the showers 
increase just as observed at sea level, and the 
energy of the individual rays is likewise similar. 
The mode of production of showers at this depth 
seems to show nothing unusual; in particular, 


there is no evidence of production by penetrating 


non-ionizing primaries. The evidence is mainly 
qualitative, the more quantitative data and the 
important change of the phenomena with depth 
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ABSORPTION EXPERIMENT 


It was mentioned in the introduction that 
some observers think the penetrating rays are 
ionizing and others believe that below 250 m 
the energy is carried down by a non-ionizing ray 
that produces softer ionizing secondaries. The 
cloud-chamber work shows that at 71 m, most 
of the rays are penetrating and ionizing. There 
might be, however, a small fraction of secondaries 
created by non-ionizing primaries. If, as the rays 
pass to greater depths, the ionizing rays are 
absorbed, the majority of the remaining rays 
observed might be secondaries from these non- 
ionizing primaries. If this is the case, the absorp- 
tion of rays present at great. depths would be 
greater with material between the counters than 
with, material above the counters. Material 
between the counters would give the absorption 
curve for the secondaries produced in the rock, 
and material above would give the absorption 
for the primaries. Figure 4 shows a Geiger- 
counter arrangement used under exactly similar 
conditions at 71 m and at 657 m. Each group 
of three tubes was connected together to form 
one large tube. Lead surrounding the second 
and third tubes shielded these tubes from side 
showers. The upper four tubes were connected to 
one coincidence circuit and the lower four to an- 
other in order to take two sets of data simulta- 
neously. The coincidences were recorded by two 
pens on a moving tape. The ceiling was one 
meter from the top group of tubes, and the same 
vein of rock was overhead at the two locations. 
Figures 5 and 6 show the absorption curves 
obtained. It is immediately apparent that at 
both depths, the observed rays are very pene- 
trating. The comparatively rapid decrease in the 
counting rates in the first 10 cm of lead was 
caused by the absorption of the soft rays. As 
pointed out when discussing the cloud-chamber 
data, this initial drop cannot be used to calculate 
the relative number of soft to hard rays, but if 
this ratio should increase with depth, then the 
initial drop would be greater at 657 m. The 
similarity of the curves indicates that the ratio 
of soft to hard rays is about the same at the two 
depths. This also suggests that the rays at 657 m 
are similar to those at 71 m. Certainly any 
possibility of a large fraction of non-ionizing 
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(a) tb) 


Fic. 3. Two knock-on electrons observed at 71 m. (a) Counter arrangement 1. 
The primary is negative, of momentum 4.7X10* ev (uw=940). The energy of 
the knock-on electron is 215 kev (e=0.43, where energy= emc*) and @ is about 
60°. For proton mass, the measured values for wo and @ give «=0.094, thus 
definitely ruling out this possibility. For k= 200, the resulting ¢ is 0.62, which 
becomes 0.45 for 6=64°, agreeing very well with the measured value. But as- 
sumption of electronic mass results in almost the same values (e«=0.48 for 
6=64°). Here then elimination of the proton is definite but no distinction is 
possible between the mesotron and electron as the primary particle. (b) Arrange- 
ment II. Primary momentum 2.9X 108, «=0.08, 62=60°. For k= 200 the observed 
¢ is obtained for @=78° while for protonic mass, @ would have to be 45°. @ is 
certainly greater than 45° so the primary is not a proton. Also the ionization 
density is too low for a ane which at this momentum value would ionize five 
times as heavily as an electron. As for (a), distinction between k=1 and 200 is 
impossible for large @ values (@= 79° gives the correct « for k=1), but here, the 
fact that the primary penetrates 10 cm of Pb eliminates the electron. Thus, in 
this case, identification as a mesotron is certain, but quantitative mass determi- 
nation is impossible. 


primaries is ruled out. With Pb only in A, a_ points Lio-s (10 cm Pb in A and 5 cm Pb in B) 
non-ionizing ray to be detected would need to and Les3-19 do not fall below the curve, thus 
produce two secondaries—one in the ceiling and further ruling out the non-ionizing primaries. 

one in A. With the same amount of Pb divided 
and placed in A and B, three secondaries would 
need to be produced to cause a coincidence. This Showers detected by counter tubes are pro- 
probability would be much lower. However, the duced in the material surrounding the tubes. 


SHOWER EXPERIMENT 
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Therefore, to make a comparison of the shower 
properties at various depths, the counter tube 
arrangement used (Fig. 7) was completely sur- 
rounded by Pb. The Pb was 15 cm thick on 
the top, 10 cm on the sides and ends, and 2.5 cm 
on the bottom. It was thought that these 
thicknesses would be sufficient to insure that all 
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Fic. 5. Absorption of cosmic rays at 71 m. 
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Fic. 6. Absorption of cosmic rays at 657 m. 


shower particles observed in the cavity would 
have originated in the Pb and not in the sur- 
rounding rock. By moving this entire arrange- 
ment from one depth to another,. any change 
observed in the nature of the showers should be 


Fic. 4. Experimental arrangement used at 71 m and 657 m to measure absorption 
and relative numbers of soft and hard cosmic rays. 


caused by a change in the rays producing the 
showers and not by the geometry of the sur- 
rounding material. 

The five tubes marked X were electrically 
connected in parallel to operate as one tube. 
Similarly, the Y tubes were connected together 
to form one large tube. Figure 8 is a block 
diagram of the circuits, each block representing 
a standard coincidence circuit. The counts from 
each were recorded upon a telephone message 
counter driven by a multivibrator. The threefold 
circuits X12, X23, X34, X45, X24, and X15 
recorded showers of two or more particles. 
Small condensers coupled the plates of the three- 
fold output tubes to the twofold coincidence 
circuits X123, X234, and X345. In this way, 
showers of three or more particles were recorded. 
Similarly, the outputs from X123 and X345 
were combined in X12345 to count showers of 
five or more particles. The outputs from the 
above ten circuits were combined with the pulses 
from the Y tubes in ten more twofold circuits 
(X Y12 X Y12345). This arrange- 
ment made it possible to obtain twenty sets of 
data simultaneously. Combining the outputs not 
only reduced the number of vacuum tubes 
needed, but also reduced the capacity connected 
to each Geiger-Mueller tube, thus short resolving 
times could be maintained. 

This apparatus was first taken to the 29th 
level, 2100 m, and operated from the same power 
line as was used for the experiment. No coinci- 
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dences were observed in fourteen hours. Ap- 
parently surges in the power line would not 
cause false counts. Tests showed that there was 
no intercoupling of the circuits which would 
cause a circuit to count falsely even when the 
other nineteen circuits were activated. Measure- 
ments of individual counting rates and of two- 
fold coincidences were made, and from these 
values, the time constants of the threefold 
circuits were found to be 3X10-* min. The 
calculated threefold chance counts were one in 
2000 days. At the surface, where both the Y 
and the X12 etc. counting rates were much 
greater, the estimated Y—X12, etc. chance 
counts were one in three days. All other chance 
counts were much lower. Data were taken at 
the surface, at 71 m, 141 m, and 582 m. 

The function of the Y counters was to 
furnish some information on the mode of 
production of the observed showers. As men- 
tioned above, it was thought that cascade 
showers could not penetrate the 15 cm of Pb 
above the X tubes, and the showers observed 
would have originated in the lead. Therefore, if 
a Y tube was discharged simultaneously with 
an observed shower, one should conclude that 
the shower was caused by an ionizing ray. Only 
about 40 percent of the area and solid angle 
above the X counters was covered by the 
Y tubes, and no measurements were made of 
the efficiency or sensitive area of the counter 
tubes; therefore, this arrangement was used not 
to determine the actual fraction of showers 
caused by ionizing rays, but rather to see if there 
is any change in this fraction with depth. 
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Fic. 7. Shower apparatus used at various depths. 
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Table IIT shows that the percentage of showers 
produced by ionizing rays does not change with 
depth. Since it is generally assumed that most of 
the showers at sea level are produced by meso- 
trons and since the cloud-chamber work at 71 m 
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Fic. 8. Circuit diagram for shower apparatus. 


depth gave no evidence of a non-ionizing shower 
producer, one may assume that practically all 
showers down to 582 m are produced by ionizing 
primaries. 

Concerning now the size of showers, it is seen 
that the threefold circuits, X12, X23, etc. 
counted all sized showers; i.e., two or more 
particles. The X123, X234, and X345 circuits 
recorded showers of three or more particles. If a 
shower discharged X12 but not X123, one cannot 
say that that shower consisted of only two 
particles, but rather that the shower was small 
and particles penetrated only two of the lower 
counters. Similarly, a coincidence X 12345 repre- 
sented a shower of five or more particles. It is 
more probable that such a coincidence was 
caused by a very large shower than by one of 
only five particles. One can only say that the 
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TABLE III. Percentage of showers accompanied by a 
discharge of a Y counter. 


Depth Percentage 
Surface 36.6+0.4 
71m 31.141.5 
141m 30.8+2.0 
582 m 38.6+6.0 


density of shower particles was sufficiently high 
for at least one particle to pass through each 
tube. We shall call the showers which tripped 
two of the lower tubes but not three (3F-4F) 
small showers; showers which tripped three but 
not five lower tubes (4F-6F) medium-sized 
showers ; and showers which tripped all the lower 
counters (6F) large showers. 

The percentage of small, medium, and large 
showers at the different depths is shown in 
Table IV. The data show quite consistently the 
steady increase in the relative number of large 
showers with depth. Furthermore, as the third 
column shows, the increase in the relative number 
is greater for the large showers than for the 
medium-sized showers. At the extreme position, 
582 m, it is interesting to note that the number 
of showers large enough to set off at least three 
of the lower tubes is actually larger than the 
number of small showers. Since the Y data 
show that there is no significant change in the 
nature of the shower-producing rays, this in- 
crease of size is most probably a manifestation 
of the increase in mean energy of the penetrating 
rays with depth. 

Table’V gives the absorption data for various 
sized showers. 


CENTIMETERS OF LEAD BETWEEN THE TUBES 


Fic. 9. Absorption of the individual particles in showers at different depths. 


Table IV shows that the large showers are 
relatively more abundant at greater depths, 
Table V states the same thing in a different way; 
namely, that the rays producing the large 
showers are more penetrating. This is to be 
expected, since the more energetic a ray, the 
larger is the shower it can produce, and also the 
greater is its probable penetration. 

These tables explain why shower absorption 
experiments performed by various observers can 
give such different results. For example, an 
experiment designed to measure shower intensity 
vs. depth will also be affected by the change in 
shower size with depth. 

To measure the absorption of the individual 
particles in the showers, counts were taken with 
various thicknesses of Pb at A, between the 
X counters and counters 1, 2, 3, 4, and 5. 


TABLE IV. Percentage of small, medium, and large showers 
at different depths. 


Depth Small Medium Large 
Surface 75 19 06 
7im 58 22 20 
141m 53 23 
582 m 45 (M+L=55) . 


*At 582m the X12345 circuit did not record because of a faulty 
vacuum tube. 
counts/hour 
Taste V. Ratio of counts/hour at surface 


Depth 3F 4F 6F 3F-6F 4F-6F 
Surface 1.00 1.00 1.00 .00 1.00 

7im 0.11 0.19 0.37 0.08 0.13 
141 m 0.04 0.08 0.18 0.03 0.05 
582 m 0.003 0.008 — .002 — 
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Figure 9 shows these absorption curves. It is 
apparent that at all three depths the individual 
icles in the large showers are less penetrating 
than in the medium-sized showers, and these, in 
turn, are less penetrating than the particles in 
the small showers. This is probably because in 
the larger cascade showers the energy of the 
original particle has been divided many more 
times, and so the average energy of the individual 
particles is less than in the small showers. Of 
course, there is no way to prevent showers from 
being formed in the Pb absorber A. These, 
moreover, would be confined to a narrow bundle 
in the Pb; and so would have a greater proba- 
bility of being detected as small showers than as 
large showers. These showers formed in A will 
not reduce in number with increasing thicknesses 
of A. This explains why the particles in the 
small showers appear to be hard. The difference 
at any one depth in the absorption curves for 
particles of various sized showers thus can be 
interpreted simply and probably is not of great 
significance. But the fact that these groups of 
curves are very similar at several depths indi- 
cates that the showers are also similar in nature 
at the different depths. At sea level, the showers 
are known to be cascade showers; apparently, 
this is also true at the other depths. 
Showers are produced by knock-on electrons 
or radiated photons. If protons are present and 
are more penetrating than the mesotrons, at the 
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greatest depth the protons alone will remain. 
Clay has pointed out that at high energies, the 
protons do not radiate as efficiently as the 
mesotrons. Thus, if the most penetrating pri- 
maries are protons, as one goes to the greatest 
depths, the relative number of showers will 
decrease. These data do not show any decrease; 
however, there may be a maximum between 141 
and 582 m, or the ratio may fall off at still 
greater depths. Plans have been made to con- 
tinue these experiments to detect the presence of 
protons underground. 


CONCLUSIONS 


The results of the various experiments of the 
present series are all consistent with the inter- 
pretation that the cosmic rays underground are 
mesotrons accompanied by electronic secondaries 
(knock-ons and their progeny). No evidence 
could be found for a non-ionizing primary. The 
data at the greatest depths are not sufficient to 
rule out the possible existence of high energy 
protons underground. 

The writers wish to express their appreciation 
to Dean A. H. Compton for his advice and 
encouragement in planning these experiments. 
It is a pleasure to acknowledge the generous 
assistance of the managers and men of the Isle 
Royale Mining Company. Mr. James Richards 
and Mr. William Boden were especially helpful 
in the maintenance of the equipment. 
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An idealized theory of cavitation in the interior of a liquid is developed as an extension of the 


hydrodynamics of irrotational motion. It is assumed that cavitation occurs whenever the 
pressure sinks to a fixed breaking-pressure and that the pressure then rises at once to a fixed 
cavity pressure. The boundary of the cavitated region either advances as a breaking-front, 
moving with supersonic velocity, or remains stationary as a free surface, or.recedes toward the 


cavitated region as a closing-front. The relevant formulas are obtained. 


AVITATION has been studied most ex- 

tensively at the interface between a liquid 
and a solid, but apparently it may occur also 
within the liquid itself. The latter type of cavita- 
tion should be influenced by the elasticity of the 
liquid and should present certain hydrodynamical 
features of interest. In the present paper a con- 
tribution is offered to the theory of cavitation as 
a hydrodynamical phenomenon. The theory may 
also find application to the motion of liquids in 
pipes, even if the cavitation is actually initiated 
at the walls.! 


1, FUNDAMENTAL ASSUMPTIONS 


In order to develop a tractable analytical 
theory, it will be assumed that the liquid 
cavitates whenever its pressure sinks to a fixed 
breaking-pressure p», and that the pressure then 
rises at once to a definite cavity pressure p,, the 
two pressures being such that 


pe. (1) 


The pressure ~ represents the minimum pressure 
that the liquid can stand; and, as a special case, 
it may equal p,. In actual cases cavitation within 
a liquid takes the form of small bubbles, and 
these commonly contain a variable amount of air 
or other gases in addition to the vapor of the 
liquid; as the bubbles collapse, the foreign gas 
tends not to redissolve entirely, and for this 
reason the pressure in the bubbles is not constant. 
Such complications will be ignored here. It may 
be assumed that the bubbles are indefinitely 
small and contain only vapor of the liquid, so 
that p. is equal to the vapor pressure. The dis- 


1 Cf. Kennard, Phys. Rev. 35, 428 (1930). 
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cussion will be limited to irrotational motion ina 
homogeneous liquid whose density remains almost 
but not quite constant. The one-dimensional case 
is easily treated completely and will be taken up 
first. 


2. PLANE BREAKING-FRONTS 


For one-dimensional motion the usual acoustic 
equations may be written 


—= — pc?—, (2a, b) 
0. 


where p = excess of pressure above the hydrostatic 
pressure, «= particle velocity, ‘= time, x=space 
coordinate, p=undisturbed density, and c=ve- 
locity of sound in the liquid. The derivatives oc- 
curring here will be assumed to be continuous 
functions. The equations are ordinarily assumed 
to hold only for small particle velocities, but they 
may also be applied generally to the relative 
motion of the liquid in the neighborhood of any 
point, that is, provided uniformly moving axes 
are employed relative to which the velocity of the 
liquid near the point is small. 

Cavitation will obviously start where the pres- 
sure first sinks to f», hence at a point at which 


(0p/dt)<0, or (du/dx)>0 (3a,b,c) 


by (2b). From this point, a breaking-front will 
advance toward both sides; that moving toward 
+x or toward the right advances into regions 
where 0p/dx=0, while the other, moving toward 
—x or toward the left, advances into regions 
where 0p/dx=0. Just ahead of each front, con- 
ditions (3a, b, c) must continue to hold. 

To find the rate of advance of the breaking- 
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fronts, consider what happens as a front moving 
toward the right traverses, during a time di, a 
layer of liquid of thickness dx (Fig. 1). At the 
beginning of dt, the pressure in the unbroken 
liquid is p» over the left-hand face of the layer F; 
and is, therefore, 


bot (p/dx)dx 


over the right-hand face F2, whereas at the end of 
dtit has sunk to ~ at F; also. Thus the change in 
pat F, during dt is (—dp/dx)dx; or, if U, denotes 
the speed of the front relative to the liquid, so 
that dx = U,dt, the change is — U,(dp/dx)dt. If, 
for the moment, axes moving with the unbroken 
liquid are employed, this change can also be 
written (0p/dt)dt, the derivative being evaluated 
at F,. Equating the last two expressions and 
using (2b), we obtain 


Us= pc?(du/dx)/(Ap/dx). (4) 


Here U; is the velocity of the front relative to the 
liquid, which may itself be in motion; but in 
evaluating the space derivatives axes moving in 
arly manner may be employed. 

_ At the front there is a discontinuity in the 
pressure. To find its impulsive effect, consider 
the motion of the layer of liquid of thickness dx 
(Fig. 1). On the right-hand face F; the pressure is 
nearly equal to the breaking-pressure p,. It sinks 
to this value as the front comes up to the face, 
whereas on the left-hand face there is acting the 
steady cavity pressure p,.. The momentum in the 
layer changes, therefore, during dt by (p.—p»)dt 
per unit area. Let denote the particle velocity 
just ahead of the front, u, that just behind it. 
Then the velocity changes impulsively from 1% to 
u, as the front passes, and since dx = U,dt, 


(De— Po)dt = (Ue — Up) pdx = (ue — up) pU rat, 
be— Po Op/dx 
du/ ax’ 


(5) 


from (4). 

The same equations are obtained for a breaking- 
front advancing toward the left. In view of (3c), 
Eq. (4) makes U; positive where 0p/dx>0 and 
negative where dp/dx <0, the two signs referring 
to fronts traveling in opposite directions. The 
sign of u,— up is likewise opposite in the two cases, 
the liquid being accelerated in the direction of 
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advance of the front. At the initial point, where p 
first sinks to dp/ax=0 and U,= +. 

It remains to be shown, however, that the 
motion of the liquid represented by u, as defined 
by Eq. (5) is a hydrodynamically possible one. In 
doing this, an additional condition for the 
propagation of a breaking-front is obtained. The 
argument takes very different courses according 
as po Or pro. 

If p.= ps, by (5) ue= the liquid being left at 
each point with that particle velocity which it 
happens to possess at the instant of arrival of the 
front. It is necessary that the liquid be not left 
with a motion of contraction; hence the values of 
u. left behind the front can only vary in such a 
way that 


(du./ dx) =0. (6) 


Now, if u has the value u; at the left-hand face F; 
of a layer of thickness dx (Fig. 1), at the instant 
when a front moving toward the right passes this 
face, then, at the same instant, at the right-hand 
face Fs, 


U= Uy + (du/dx)dx. 


Let uniformly moving axes be employed such 
that momentarily u; is zero or small. Then, by 
the time that the front has reached F;, u will have 
changed at F; to 


(Ou /dx)dx+ (du/dt)dt. 


dx = Updt 

fe) 
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Fic. 1. Motion of a breaking-front. 
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After the passage of the front, u does not change 
with time because of the absence of a pressure 
gradient in the cavitated region. Hence the 
difference, 11, is also the difference between 
the values of u, at F; and at Fs, or du,, so that 


due = = (Ou /dx)dx+ (du/dt)dt, 
or, by (2a) and the relation dx = U;dt, 
due Up(du/dx) — (1/p)(dp/dx) jdt. 
The condition expressed by (6) now requires that 
Us=(1/p)(dp/dx)/(du/dx), (7) 


since du/dx>0. Multiplication of this inequality 
by Eq. (4) gives 
(8) 


and division of this latter inequality by Eq. (4) 
squared gives 


12 (1/p*c)[(dp/dx)/(du/dx) (9) 


whence 
| dp/dx| = pc(du/dx). (10) 

The same results are obtained for fronts 
moving toward the left. 

If, on the other hand, p. > p», the argument just 
given does not apply, because by Eq. (5) u. up 
at the front. The situation appears to be es- 
sentially different. For a front running to the 
right (du/dx>0, dp/dx>0), so that the 
liquid behind the front tends to overrun that just 
ahead of it, which is physically impossible. The 
condition that prevents this happening can be 
found as follows. 

While the front is traversing the layer of 
thickness dx (Fig. 1), during a time dt, the face F; 
of the layer advances a distance ud? (to the first 
order), whereas the other face F; advances a 
distance u.dt. The volume occupied by the layer, 
per unit area, thus decreases by 


fo) /(eUs) dt (11) 


by (5). During the same time, however, the 
liquid in the layer contracts, as its pressure rises 
suddenly from ~ to p.. The elasticity being 
pc*, the contraction decreases the volume actually 
occupied by liquid by the amount 


(1/pc*)(be— po)dx = (1/pc*)(pe— pe) Urdt. (12) 
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Hence the liquid can continue to fit into the 
space allowed for it, with or without the occur. 
rence of cavitation, provided 


(1/pc*)(pe— po) bo) /pUs], Us? 


Thus expressions (8) and (10) are obtained again, 

In case p.>p», however, it is not possible to 
show that necessarily du./dx2=0. Apparently it 
can happen that behind a front du,/dx <0; in 
this case, any cavitation bubbles that may have 
been formed at the instant of passage of the front 
will proceed to decrease in size and ultimately to 
disappear again. 

According to Eq. (8) as written it can happen 
that U,=-c. It is easily seen, however, that if 
U,=-c for a finite length of time, cavitation 
does not actually occur. Thus a true breaking. 
front travels always at a speed exceeding that of 
sound in the liquid. The conditions found for its 
propagation, including (10), may be summarized 
as follows: 


(du/dx)>0, 
|dp/dx| <pc(du/dx) ; 


(13a, b, ¢) 


in addition, it is necessary that 0p/dx2=0 for. 
propagation toward +x or 0p/dx=0 for propa- 

gation toward —x. In these expressions the 

derivatives are to be evaluated in the unbroken 

liquid just ahead of the breaking-front. These 

values are not affected in any way by the ap- 

proach of the front; for no hydrodynamic influ- 

ence of the first order can be propagated through 

a fluid at a speed exceeding that of sound. 

Once started, a breaking-front will travel until 
it arrives at a point beyond which the necessary 
conditions are not satisfied. At a point where 
dp/dax=0, the front may meet another one 
coming from the opposite direction, whereupon 
both fronts will disappear. At a point where 
|ap/dx| becomes equal to pcdu/dx, with larger 
values beyond, the front, approaching at speed ¢ 
relative to the liquid (provided du/dx>0), must 
suddenly stop advancing. What happens next at 
such a point will be considered presently. 

A simpler rule than that just stated might have 
been supposed to be preferable, namely, that the 
front travels as far as pressures so low as » occur. 
This statement, however, although true, is less 
convenient than the conditions just stated, be- 
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cause the occurrence of cavitation results in- 
directly in a contraction of the region within 
which the pressure sinks to p». This results from 
the fact that the pressure is higher in the 
cavitated region than it would have been if 
cavitation had not occurred, and, after the 
breaking-front has ceased advancing, propaga- 
tion of this higher pressure into adjoining regions 
prevents the occurrence of pressures so low as pp» 
at some points where such pressures would other- 
wise have occurred. It is impossible, in fact, for 
the breaking-front to touch the surface which, 
in the absence of cavitation would be the bound- 
ary of the region of pressures below ),; for, at 
this surface, when p= pp», 6p/dt=0, else p would 
sink to p, at points outside the surface as well 
(provided 0p/dt is continuous), and then by (2b) 
du/dx=0 and condition (13c) is not satisfied. 


3. THE CAVITATED REGION 


In the cavitated region the pressure is assumed 
to be uniformly equal to p., so that there is no 
pressure gradient; hence the particle velocity u, 
at each point remains as it was left by the passage 
of the breaking-front. Let n denote the fraction of 
the space which is occupied by bubbles. If 
p.= po, nis zero immediately back of the breaking- 
front. If p.>», however, n starts from an initial 
value no representing the difference between the 
space freed by compression from pressure p» to p. 
and the space filled by incoming liquid. When a 
front travels toward +x, the intial free space in a 
layer dx thick just back of the front equals the 
difference between expressions (12) and (11), so 
that, since dx = U;dt, 


nodx = (1/pc*)(pe— po) Usdt—[(be— Us jdt, 
no=[(be— pc? (14) 


The same result is obtained for fronts traveling 
toward —x. Thereafter » changes at the rate 
du./dx. Hence, at any time ¢ subsequent to the 
passage of the front at time &, 


n= f (du,/dx)dt 


The value of » may vary from point to point in 
the cavitated region. The initial value mo will be 


small, in accordance with the assumption that 
the pressure remains within the linear range; but 
subsequently » may increase to any magnitude. 


4. MOTION OF PLANE BREAKING-FRONTS 
IN TERMS OF WAVE TRAINS 


Since no influence can be propagated with 
supersonic velocity past a breaking-front, the 
pressure and particle velocity ahead of the front 
will be determined by initial or boundary con- 
ditions elsewhere; hence, in the above formulas 
the values of p and u may be regarded as given. 
When, however, the boundary of the cavitated 
region moves more slowly than the speed of 
sound, the presence of cavitation is able to 
influence the values of » and wu in the nearby 
unbroken liquid. In such cases analysis in terms 
of fundamental wave trains is convenient. It will 
be worth while to restate, also, the conditions for 
the propagation of a breaking-front in terms of 
such trains. 

Any one-dimensional disturbance in unbroken 
liquid can be resolved into a train of waves 
traveling at speed c toward +x and another 
traveling toward —x; if p1, p2 denote the corre- 
sponding pressures and 4; and w%, particle 
velocities, then 


(17a, b) 
pr=}(p—pcu). (18a, b) 


Here the significance of ; and p» is really only 
mathematical, hence “:, “2, and u need not be 
restricted to be smaiJ. Substitution from these 
equations for u in (13b) gives 


(dp1/dx) — (dp2/dx) >0, 
and substitution in (13c) squared gives 
[(dp1/dx) + (dp2/dx) <[(dp1/dx) — (dp2/dx) }*, 
<0. 


The first and last of these inequalities can be 
satisfied simultaneously only if 8p:/0x is positive 
and dp:/dx negative. Hence the conditions for 
the propagation of a breaking front and Eqs. (4) 
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and (5) for U, and u, can be written, since w=, 


Pitpe2= ho, (0p: /dx) >0, 


(19a, b, c) 
(Op2/dx) <0; 
and, for travel toward +x, 
(dpi/dx) = — (dp2/dx), (20a) 
or, for travel toward —x, 
= — (dp2/dx) ; (20b) 


+ 


1 
— pot (pe— po) 
(Op1/dx) + (dp2/ =") 
(8p1/dx) — 


Equations (19a, b, c) lead to the useful con- 
clusion that a wave of tension can give rise to 
cavitation only where it overruns a wave of 
decreasing pressure traveling in the opposite 
direction. . 


5. PLANE RESTING BOUNDARIES AND 
CLOSING-FRONTS 


The next question to consider is the behavior 
of the boundary of a cavitated region when it 
cannot advance as a breaking front. It may stand 
still, relatively to the liquid, or it may advance 
into the cavitated region, leaving the liquid 
unbroken again behind it; in the latter case, the 
boundary may be called a closing-front. Which 
alternative occurs will depend, in general, both 
upon conditions in the cavitated region and upon 
the magnitude of the wave that is incident from 
the unbroken side. | 

Let p’ denote pressure in the component wave 
train that approaches the boundary in the 
unbroken liquid and ” that in the receding 
train; let « denote the particle velocity on the 
unbroken side and u, that on the cavitated side, 
but let the positive direction for velocity be taken 
now always toward the cavitated region. As 
before, let », denote the uniform pressure in the 
cavitated region and 7 the fraction of the space 
that is occupied by bubbles. 

Whether the boundary remains at rest or 
advances is found to depend upon the magnitudes 
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Fic. 2. Layer of cavitated liquid dx thick. 


of p’, u-, du./dx, and n. These quantities may be 
regarded as given; p” is then determined by con- 
ditions at the boundary, and p and wu are given by 
the relations, obtained from (17a, b) and (16a, b), 


pb=p'+p", u=(1/pc)(p’—p”), (23a, b) 


whence 
p+pcu=2p'. (23c) 


Four cases may be distinguished. 


(1) Intrinsic Closing-Fronts 
n=0, odu./dx<0, 


(24a, b, c) 
—du,/dx=cdn/dx>0. 


Here x denotes distance measured from the 
boundary into the cavitated region, and the 
relations refer to conditions at the boundary. 

Because of (24a, b), the cavitation will proceed 
to disappear in the layers next to the boundary; 
hence the boundary will advance into the 
cavitated region as a closing-front, leaving 
unbroken liquid behind it. A front of this type, 
being propelled by conditions inside the cavitated 
region, may be called an intrinsic closing-front. 

To find its speed of advance, consider a layer 
of the cavitated liquid dx thick (Fig. 2). If the 
particle velocity is u, over the face toward —«x or 
F,, that over the other face F» is 


(du./dx)dx. 
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Hence during a time dé the volume of the layer 
changes by (du./dx)dxdt, for each unit of area of 
its faces and, since the pressure is constant, the 
change must occur in the cavitated space and not 
in the liquid itself. Hence, 


(du,/dx)dxdt = (0/dt)(ndx)dt, 


(25) 
dn/dt=du,/dx, 


in which the time derivative is to be evaluated at 
a point moving with the liquid. 

On the other hand, if U, is the velocity of 
advance of the boundary relatively to the liquid, 
dx= U.dt. At the beginning of dt, n»=0 at F,, 
hence at Fs »=(dn/dx)dx; whereas at the end of 
dt, n=0 at F,. Hence, during dt, » changes at F, 
by dn=—(0n/dx)dx. This can also be written 
dn=(dn/dt)dt. Hence, using dx = U.dt, we obtain 


(dn/dt) = — U.(dn/dx), (26) 
and 


U.= —(n/dt)/(dn/dx) 
= (27) 


by (25). Here the positive directions for x and u, 
are taken toward the cavitated region and the 
signs of the derivatives are such as to make U, 
positive. 

So far as can be seen, the value of U, as given 
by (27) might be either greater or less than c. If 
it is greater than c, everything is satisfactory. 
The advancing front leaves behind it a region in 
which p= p, and u=4,, u, being the local value of 
the particle velocity as encountered by the 
advancing front. The two component wave trains 
in the unbroken liquid are thereby progressively 
built out in accord with the equations 


p’=3(p+ pcu) =3(pe+ pcue), 
=3(p— peu) = }(pe— pete). 


(28a) 
(28b) 


The part of the p’ wave train that is already 
established and is moving toward the front is 
unable to overtake it. Even if U.=c, there is no 
difficulty, although in this case a fixed discon- 
tinuity in p and u may occur at the front. 

If, however, the value of U. as defined by 
Eq. (27) is less than c, a front moving at speed U. 
would be overtaken continually by values of p’ 
coming out of the unbroken region. In such cases 
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conditions at the boundary must be determined 
in part by the arriving values of p’, and they will 
usually not be such as to satisfy Eqs. (28a, b). A 
different type of action must then occur at the 
boundary. Thus intrinsic closing-fronts, when they 
occur, move through the liquid at velocities not less 
than that of sound. 


(2) Resting boundary 
2p’=p.+ (29) 


Let it be assumed that the pressure is continuous 
at the boundary and equal, therefore, to p,.. Then 
a wave train must be reflected from the boundary 
of such magnitude that, by (23a, c) 


b"=p.—p', u=(1/pc)(2p'—p-). 


These equations represent a possible motion pro- 
vided u=u,, so that the unbroken liquid does not 
collide with the cavitated liquid; and this con- 
dition, in turn, is met provided p’ is such that 
(29) holds. It may be concluded that, whenever 
(29) holds but the boundary cannot advance as an 
intrinsic closing-front moving at supersonic ve- 
locity, the boundary will remain stationary and 
the incident waves will be reflected from it in 
accord with Eq. (30a). 

If 2p’ is actually less than p.+pcu,, so that by 
(30b) u<u,, separation of the two portions of the 
liquid must occur, with formation of a crevasse or 
a special layer of expanding bubbles. If 2p’=p, 
+ pcu., the unbroken liquid remains in peaceful 
contact with the cavitated liquid. 

If, however, (29).does not hold, no solution of 
the type just described is possible and a different 
process must occur. 


(30a, b) 


(3) Forced Closing-Fronts 


2p’>petpcue, n>0. (31a, b) 


When condition (31a) obtains, the incident waves 
are so strong that the boundary advances into the 
cavitated region, closing up the bubbles. The 
action is one of successive impacts, a finite 
discontinuity of pressure occurring at the moving 
boundary if »>0. A boundary moving in this 
manner may be called a forced closing-front. 

Let U. be the speed of advance of the front 
relative to the cavitated liquid ahead of it. 
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Then, during a time dt, a layer of cavitated liquid 
of thickness dx becomes consolidated. The argu- 
ment is now similar to that in Section 2. 

The mass (1—7)pdx of liquid per unit area is 
accelerated from u, to u; the added momentum is 
supplied by the action of a pressure p on one face 
and p, on the other face. Hence, 


(1—)p(u—u.)dx = (p—p.)dt, 
and since dx = Ut, 
(32) 
where » and wu denote pressure and particle 


velocity in the unbroken liquid just behind the 


front. 
Furthermore, since one face of the layer moves 


at speed u and the other at u,, its thickness 
changes during dt by (u.—u)dt. Part of this 
change is taken up by compression of actual 
liquid from p,.to p, which changes the total 
thickness of liquid, initially equal to (1—7)dx, by 
the amount 


](1—n)dx, 


pc? representing the elasticity. In addition, 
cavitated spaces of total thickness ndx are closed 
up. Hence, 


(ue—u)dt= ndx, 
pc?(u— ue) ]Ue. 
In addition, there is Eq. (23c) or 
b+ pcu=2p'. (34) 


Equations (32), (33), and (34) suffice to fix U., 
p, and u. The solution can be written in a con- 
venient explicit form by introducing the auxiliary 
quantity 


(33) 


v= (1/pc)(2p’ — pc) — 


According to (31a) v is positive; it may be re- 
garded as known. Then 


U.=cv/[v+(c—»)n], (35) 
(1—1) pcv? 
20+ (c—20)n" 


(36a, b) 


In these equations, asin all equations not referring 
to breaking-fronts, u and u, are taken positive 
toward the cavitated region. 

Since v>0 and 12=>0, these formulas make 
p>p- and u>,, as is obviously necessary for 
physical reasons. Furthermore, from (35) it is 
clear that U,<c. Or, if we consider the velocity of 
the front relative to the unbroken liquid behind 
it or U.’=U.—(u—u.), we find that U,’'<¢. 
Thus the front cannot run away but must remain 
subject to influences propagated up to it from 
behind at speed c. To sum up, all forced closing. 
fronts move at a speed less than that of sound 
relative either to the unbroken liquid behind them or 
to the cavitated liquid in front. 

The requirement, 7>0, has been included in 
the definition of forced closing-fronts because 
boundaries at which »=0 appear to constitute a 
singular class. One type of such boundaries has 
already been discussed. 


(4) Boundaries at which n=0 


If 2p’>p.+pcu,, a little consideration shows 
that the boundary will advance as a closing-front 
of theintrinsictype if it can do so, otherwise as one 
of the forced type; in the latter case, its velocity 
of advance is initially equal to c, so that any 
tendency to form a front closing by intrinsic 
action but advancing at a speed below c is cut 
short. The front will usually advance at once into 
regions where n>0. 

If, however, 2p’<p.+pcu,, a resting boundary 
as described by Eqs. (30a, b) should certainly 
occur, accompanied by a finite discontinuity in 
the particle velocity and by the formation of a 
gap or a special layer of bubbles lying between 
the unbroken and the cavitated parts of the 
liquid. It appears, however, that the opposite 
side of this gap may also advance into the 
cavitated region as a closing-front, so that the 
original boundary splits into three, two resting 
boundaries separated by a gap and a closing-front. 

The principal type of boundary at which 7=0 
is the boundary that has just ceased moving asa 
breaking-front. For there, usually, the inequality 
expressed by (13c) has just failed to exist because 
its members have momentarily become equal; 
hence by (4) at this instant, U;2?=c* and by (14) 
n=no=0. At this instant, furthermore, Eq. (5) 
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CAVITATION 


becomes, after reversal of the signs of u. and u% 
because of the difference in the direction chosen 
here for positive velocities, 


Ue = Un — (Pe— Po) / pe, 


and insertion of the value of ™ given by this 
equation for u and of p=, in (23c) gives 


2p’ Pet 


The slightest increase in p’ will thus cause the 
boundary to start back as a closing-front. 
Otherwise it remains as a resting boundary. 


6. THREE-DIMENSIONAL CAVITATION 


Insofar as the three-dimensional theory can be 
formulated in compact form, it is easily written 
out as a generalization of the one-dimensional 
theory. 

Equations (2a, b) are replaced by 


du/dt=—Vp, dp/dt=—pc?divu. (37a, b) 

It is clear that wherever 0u/dx occurs in the one- 
dimensional formulas as the equivalent of dp/dt, 
it is to be replaced by div u. The theory of 
boundaries as given above then becomes appli- 
cable to three-dimensional cavitation, provided x 
in the one-dimensional equations is interpreted as 
distance along the normal to the boundary, and 
provided u is replaced by the component of the 
particle velocity taken in the direction of this 
normal. In addition, it is to be assumed that, as a 
boundary moves, no change occurs in any com- 
ponent of the particle velocity tangential to it, 
the impulsive changes being confined to the 
normal component. 

Cavitation will begin at a point where a local 
minimum of pressure sinks to the breaking pres- 
sure ~», and from this point a breaking-front will 
spread out in the form of a closed surface sur- 
rounding the point of initiation and enclosing a 
cavitated region. In the unbroken liquid this 
surface, at which p=», constitutes a surface of 
uniform pressure, hence the normal to the surface 
lies everywhere in the direction of Vp. The con- 
ditions to be satisfied in the unbroken liquid in 
order that such a front may advance can at once 
be written down as a generalization of Eqs. 
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(13a, b, c): 


p=p, divu>d, 
OS(0p/dn)<pedivu, (38a, b, c) 


where m stands for distance measured into the 
unbroken liquid along the normal to the bound- 
ary. As before, the front cannot travel to the 
boundary of the region in which, in the absence of 
cavitation, the pressure would sink to p», at least 
provided all derivatives of » are continuous 
functions of the space variables. For U;, the 
velocity of propagation of the front along its 
normal relative to the unbroken liquid ahead, 
Eqs. (4) and (8) are replaced by 


In the same way the theory of the cavitated 
region, as developed in Section 3, can at once be 
generalized. Those components of the particle 
velocity which are tangential to the breaking- 
front are unaltered by the passage of the front, 
whereas the component perpendicular to the 
front, taken positive toward the unbroken liquid, 
is changed from to Ucn where 


Ucn = Urn + (Pe— Po) / pUr. 


The fractional cavitation , or the fraction of the 
space that is occupied by bubbles, is left by the 
passage of the front at the value stated in Eq. (14) 


or 
no=[(be— }L1 — (c?/ Ue?) J; 


after the lapse of a further time t—4 it becomes 


(40) 


(41) 


f div 


(42) 
Here u, is the particle velocity in the cavitated 
region, which remains constant in time as long as 
the cavitation lasts. 

From the physical standpoint, the behavior of 
resting boundaries and of closing-fronts, also, is 
essentially the same as in the one-dimensional 
case. For intrinsic closing-fronts the mathe- 
matical theory likewise is similar and can at once 
be written down. Equation (27) for the propa- 
gation velocity of such a front, relative to the 
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cavitated liquid, is replaced by 
U.= — (div u,)/(On/dn). (43) 


For a forced closing-front, the analogs of 
Eqs. (32) and (33) can at once be written: 


P—pe=(1—1) pU-(Un— ten), (44) 
pc? (tn — Uen) = [pc*n + (1 — n)(P—pe) ]Uc. (45) 


Here p and wu denote pressure and particle ve- 
locity on the unbroken side of the front and the 
subscript m denotes the component normal to the 
boundary, taken positive when directed toward 
the side of cavitation. Components of velocity 
tangential to the front are unaffected by its 
passage. Difficulty is encountered, however, in 
attempting to obtain a third equation in order to 
determine p, u, and U., if 7 and u,, are assumed 
known. Even the condition for the distinction 
between resting boundaries and forced closing- 
fronts cannot easily be formulated in the general 
case. The disturbance in the unbroken liquid 
might, indeed, be resolved into an infinitude of 
plane waves; but the treatment of these waves at 
a curved boundary is rendered very complicated 
by the occurrence of diffraction. All that can 
readily be stated is that existing conditions in the 
unbroken liquid, which in turn depend in part 
upon conditions at distant boundaries, in combi- 
nation with the existing distribution of the values 
of » and of the particle velocity in the cavitated 
region, will determine whether the boundary 
advances or remains at rest; and, if the boundary 
advances, those conditions in the unbroken liquid 
combine with Eqs. (44) and (45) to fix the values 
of p, un, and U,. 

It is unfortunate that the three-dimensional 
theory cannot be carried further in simple form. 

Aconcluding comment may be added, however, 
concerning the special case of spherical waves. © 

Let 


1 
pi=-f(ct—r) 


represent the pressure in a train of waves diverging 
from a fixed center, r being the distance from the 
center and f any differentiable function. Let 
these waves be superposed upon other waves of 
small amplitude represented as a whole by 


p2(x, z, t). 


In this case relation (38c) may conveniently be 
written, by means of (37b) and the relation, 


p=piths, 
Opi oh Opi = 


—-—-—>c 


ot oat 


If the normal along which m is measured makes 
an angle @ with the direction of r, 


Opi Opi pr 
6 ( 
on Or 


hence, the preceding equation, expressing the 
necessary condition for the propagation of a 
breaking-front, may be written 


0 0 


From this result it is clear that a breaking-front 
cannot follow a diverging spherical wave of 
tension into quiescent liquid (p:<0, p2=0); for 
in sucha case @=0, since over the front p= pp = 
To cause cavitation, a spherical wave of tension 
must overrun other waves of such character that 
the right-hand member of Eq. (46) has sufficiently 
large negative values. 


7. ENERGETIC AND THERMODYNAMICAL 
CONSIDERATIONS 


The preceding equations have been based 
solely upon the conservation of matter and of 
momentum. An investigation, of which no details 
will be given, reveals, however, that, if p,<p,, 
the passage of either a breaking-front or a forced 
closing-front is likely to leave behind it a greater 
amount of energy than is accounted for by the 
elastic and kinetic energy of the cavitated water. 
Upon closer examination of the cavitation process, 
it appears probable that this excess energy is 
either dissipated or radiated away in consequence 
of local oscillations associated with the impulsive 
production or destruction of bubbles. Thus the 
theory involves no conflict with thermodynamics. 

A difference between the breaking-pressure pp 
and the cavity pressure p, implies, however, an 
instability of the liquid which might seem to 
open the way for other types of action than those 
assumed above. There exists in fact, no hydro- 
dynamic or thermodynamic reason why cavita- 
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tion should not occur at any pressure below p.. An 
analogy might be expected with the behavior of 
moderately supercooled liquid, in which freezing 
must be initiated by external causes but, once 
started, spreads rapidly until the temperature of 
the mass has been raised to the freezing point. It 
might be anticipated, therefore, that a breaking- 
front, if unable to advance further at the pressure 
of spontaneous cavitation, p», would continue to 
advance at higher pressure pp’ until it encountered 
a pressure equal to p<. 

Against such an assumption the following ob- 
jection may be raised. It can be shown that even 
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on this alternative hypothesis a breaking-front 
must advance through the liquid at supersonic 
velocity. But the advance into regions of pressure 
higher than ~, would necessarily occur, not as a 
consequence of automatic changes in the unbroken 
liquid ahead, but as an effect propagated out of 
the region already cavitated; and it seems 
unlikely, although certainly not impossible, that 
such an effect could be propagated at a speed 
exceeding that of sound. 

Which assumption corresponds more nearly 
with the behavior of actual liquids remains to be 
discovered by experiment. 
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With a ‘“‘two-bulb” thermal diffusion experiment in which is used NH; having about a 15 per- 


cent N® content to increase the accuracy of the mass spectrometer analyses, it has been found 
that with decreasing temperature the thermal diffusion constant a of ammonia changes from + 
to — values at about room temperature. The value of a varies linearly with the logarithm of the 
absolute temperature, the rate of decrease being, however, nearly eight times that for neon 
and argon. A qualitative discussion is presented, attributing this effect largely to the strong 


N general the phenomenon of thermal diffusion 
operates to produce an increased concen- 
tration of the heavier component in a binary gas 
mixture in the colder portion of the apparatus, 
while the lighter component tends to concentrate 
in the warmer part. For these cases there is a 
positive thermal diffusion constant a in the 
equation Dr/D=acic2, where Dr is the coef- 
ficient of thermal diffusion, D is the coefficient 
of ordinary diffusion, and c, and ¢2 are the 
relative concentrations of the light and heavy 
components, respectively. In fact, until just 
recently no instance where thermal diffusion 
proceeds with a negative thermal diffusion con- 
stant had been reported. Grew! has now found 
a reversal in the sign of a with change in the 
composition of a neon-ammonia mixture. Be- 


'K. E. Grew, Nature 150, 320 (1942). 


first-order dipole-dipole intermolecular forces which are proportional to 1/R*. 


tween 0 and 75 percent neon, the heavier mole- 
cule, neon, tends to concentrate at the upper end 
of a Clusius-Dickel column. This indicates that 
the a is negative. Above 75 percent neon, the 
neon concentrates at the lower end, so that in 
this range @ is positive. The theoretical possi- 
bility of such a reversal had been pointed out by 
Chapman.? 

A definite temperature variation of the positive 
a for neon was found by Nier,* the value of a 
decreasing with the lowering of the mean tem- 
perature. Jones‘ has shown theoretically that a 
variation similar to that reported by Nier is to 
be expected for either the Sutherland or the 
Lennard-Jones 9,5 molecular models. This theo- 


retical investigation indicated that, for tem- 


2S. Chapman, Proc. Roy. Soc. Al77, 38 (1940). 
3A. O. Nier, Phys. Rev. 57, 338 (1940). 
*R. Clark Jones, Phys. Rev. 59, 1019 (1941). 
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peratures near the critical temperature for a gas 
whose molecules are spherically symmetrical, the 
thermal diffusion constant might be negative. 
Recently Stier’ has made measurements for a 
for argon over a number of temperature ranges 
in the two lowest of which a should be negative 
according to the Lennard-Jones 9,5 model. The 
observed a’s are positive, however. 

The present investigation with ammonia was 
initiated with the object of finding an experi- 
mental case where, with decreasing temperature, 
a reversal in the sign of the thermal diffusion 
effect occurs. Ammonia was considered to be a 
good possibility because (1) predictions from 
viscosity data indicate that a should be negative 
in the lower temperature ranges, ® (2) it is possible 
to work well below the critical temperature, 
132°C, and (3) it is known that long-range inter- 
molecular forces are particularly strong where 
the molecules are symmetrical tops with a rela- 
tively large dipole moment.’ Since the ammonia 
molecule is not spherically symmetrical and 
since there is the possibility of an interchange of 
energy between the external and internal degrees 
of freedom, most of the theoretical equations 
which have been developed for the thermal 
diffusion effect do not strictly apply. But thermal 
diffusion depends strongly on the intermolecular 
forces between the molecules involved, and hence 
any measurements of controlled variations of 
this phenomenon should contribute to our better 
understanding of these forces. 


EXPERIMENTAL PROCEDURE 


The experiments have been performed with an 
all-glass ‘“‘two-bulb” apparatus, both bulbs being 
of large volume, 2 liters, in order to have a suf- 
ficient amount of gas for a number of cycles of 
operation in spite of the low initial pressure of 
4.6 cm of Hg set by the temperature of the cold 
bulb, —76°C. The two bulbs were connected by 
a 1-cm bore tube 34 cm long, containing a three- 
way stopcock also of 1-cm bore. Appropriate 
temperature baths maintained the bulbs at the 
required constant temperatures during each 
experiment. After some trials good data were 


5 L. G. Stier, Phys. Rev. 62, 548 (1942). 
6 H. Brown, Phys. Rev. 57, 242 (1940). 
7H. Margenau and D. T. Warren, Phys. Rev. 51, 748 


(1937). 
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obtained for the temperature intervals 197°K-~ 
298°K, 197°K-373°K, and 298°K-457°K. These 
results are presented below. 

Both to overwhelm any residual N“H mass 
peak (all NH; samples were converted to Ne and 
analyzed for the N'%/N* ratio) in the mass 
spectrometer and to develop larger changes in 
concentration in the thermal diffusion experi- 
ments, ammonia gas was prepared from samples 
of NH,NO; in which the N of the ammonium 
radical was enriched in the N® isotope. This 
“heavy” NH,NO; was obtained from the 
Eastman Kodak Company and was diluted in 
every case with ordinary ammonium nitrate so 
that the prepared ammonia contained from 13 
to 15 atom percent N. The NH; was obtained 
by slowly adding a concentrated solution of 
NaOH, the reaction taking place in a vacuum, 
The gas was condensed in a dry-ice trap and then 
thoroughly dried by vacuum distillations and 
drying agents. A sample of this ammonia was 
then taken for mass spectrometer analysis, and 
the two-bulb apparatus filled with the gas toa 
pressure of 4.6 cm of Hg. The hot and cold bulbs 
were then interconnected and allowed to come 
to equilibrium. 

It can be shown that the relaxation time 7 for 
this process is 


+T-)], (1) 


where V is the volume of one of the bulbs, / is the 
length, and A the cross section of the connecting 
tube, D is the coefficient of ordinary diffusion, 
and T, and 7) are the temperatures of the cold 
and hot bulbs, respectively. For 7.=197° and 
T, = 373°, r is 3.2 hours. To insure that equi- 
librium was reached, the first thermal diffusion 
was continued for about five times this value of r. 
The cold bulb was then isolated, and the am- 
monia pumped out of the hot bulb. Next, the 
connection between the two bulbs was reestab- 
lished and the process repeated. Since D is 
inversely proportional to the density, r decreases 


TABLE I. Values of a and Rr for ammonia. 


Te°K Ta°K cc? a Rr 


298 457 0.1379 0.8621 +0.0017 +0.0105+0.0051 +0.41+0.20 
197 373 0.1525 0.8475 —0.0009 —0.0039+0.0035 —0. 
197 298 0.1419 0.8581 —0.0016 —0.0100+0.0063 —0.39+0.25 
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a vs. InTr 


Fic. 1. Variation of the thermal diffusion constant a for am- 
monia with the logarithm of the absolute temperature. 


proportionally with the pressure drop each time. 
This cycle of operations was repeated eight times 
in each experiment. For greater accuracy it 
would have been desirable to have repeated 
this process more than eight times, but there was 
a limitation on the initial amount of gas set by 
the vapor pressure of NH; at —76°C, and in the 
conversion to Nz before analysis the pressure of 
the N2 is only half that of the NH; converted. 
Before analysis we found it necessary to have 
this nitrogen in contact with a liquid-air trap 
for about 12 hours in order to condense all water 
vapor evolved in the conversion and any remain- 
ing non-oxidized ammonia. 

The N"“/N* ratios in the original and changed 
samples of gas were obtained by identical pro- 
cedures using a modified Nier-type 60° mass 
spectrometer. Our modifications of Nier’s original 
specifications* have been made to increase the 
stability of the instrument and to facilitate its 
operation. In order to be able to replace a 
burned-out filament more quickly, the filament 
stem is introduced into the side of the ion source. 
The electrometer tube is placed in a metal case 
in the evacuated system, with the ion collector 
placed at the top of the grid of the tube. The 
connections to the elements of the electrometer 
tube are brought into the case through the base 
of a 6L6 tube soldered into the metal case. Since 
this case may not be raised to as high a tem- 
perature as the other parts of the system during 
the preliminary baking out, its walls are cooled 
with water, ice, or dry ice to a constant lower 
temperature while in operation. 


*A. O. Nier, Rev. Sci. Inst. 11, 212 (1940). 


RESULTS 


The experimental values of a may be calcu- 
lated from the equilibrium concentrations of the 
two isotopic molecules in such a two-bulb experi- 
ment from the relation 


a= x (2) 
In (T,/T-) 


where C2 is the concentration of N™ in the cold 
bulb, c2® and c,;° are the concentrations of N™ 


and N*, respectively, in the gas before the 


thermal diffusion began, is the number of 
cycles of operation, and 7, and T, are the 
temperatures of hot and cold bulbs, respectively. 
Rr, the index of intermolecular force, is taken 
to be the ratio of this experimental a to the hard 
sphere value }. 
The latter is 0.0254 for this 
mixture. The measured concentrations and cal- 
culated a and Rr values are given in Table I. 

The magnitude of the maximum possible error 
assigned to each of these a’s needs explanation. 
For each temperature range only one thermal 
diffusion experiment is represented ; the concen- 
trations given are the average of at least ten 
determinations of the N“/N™ abundance ratio 
in each case. The average deviation from the 
mean for this ratio in each case is about 0.0004, 
which makes the possible error for each c2.—c2° 
difference +0.0008. Improvement in this accu- 
racy could have been made by averaging in data 
from additional runs, of course, but we feel that 
our values of a are considerably more accurate 
than would seem to be indicated by these tabu- 
lated possible errors. For in Fig. 1 in which these 
a’s are plotted against In T,, where 


is the correct mean temperature® for the tem- 
perature interval 7,—T., it is evident that a 
linear relationship holds just as for neon and 
argon.® If Fig. 1 is changed into a Rr vs. In T, 
plot, the equation fitting the resulting straight 
line is 

Rr=1.9 In (T,/293). (4) 


*H. Brown, Phys. Rev. 58, 661 (1940). 


183 
°K- O15 
293°K 
$9 
2s in 
peri- =.010 
iples 
ium 
a 
the 
d in 
te so 
n 13 
ined 
n of 
then 
and 
was 
and 
toa 
yulbs 
t for 
(1) 
s the 
cting 
sion, 
cold 
equi- 
ote (3) 
of 
am- 
, the 
stab- 
D is 
eases 
RT 


7: 


184 W. W. WATSON AND D. WOERNLEY 


TABLE II. Comparison of Rr values. 


Rr Rr(n) Rr(Suth.) Rr(Suth.) Rr(L-J 9,5 
Tr obs. viscosity C =500° C =370° e/k =333° 


366° —0.15 0 —0.27 
268° =—0.15 —0.17 —0.30 —0.13 —0.41 
239° = —0.39 —0.33 —0.22 —0.42 


Of chief interest are the observed reversal in 
the sign of Rr at 293°K, the linear variation of 
Rr with In 7, even though NH; is far from 
being a spherically symmetrical “‘hard”’ molecule 


. like neon, and the size of the slope, 1.9, in Eq. 


(4) as compared to 0.25 for neon and argon.’ 
These observations indicate a very rapid ‘‘soften- 
ing” of the NH; molecule as the temperature 
decreases. 


DISCUSSION 


The temperature dependence of Rr as pre- 
dicted by viscosity data and by the Sutherland 
and Lennard-Jones 9,5 molecular models has 
been considered by Jones.‘ Since the predictions 
on these models have been shown to be rather 
poor for neon and argon,’ they could hardly be 
expected to hold for a polyatomic molecule 
such as ammonia. There is reason to believe, 
however, that the 9,5 model might be a fair 
approximation to the real intermolecular forces 
existing in ammonia. Comparison of such pre- 
dicted Rr values with our observed values is 
made in Table II. For the theoretical relations 
upon which these calculations are based reference 
should be made to the paper by R. C. Jones.‘ The 
exponent n of T to give the coefficient of viscosity 
is taken to be approximately 1.1 in this entire 
temperature range.'® From this a negative value 
of Rr close to that observed for our intermediate 
T, is predicted. The Sutherland constant, 
C=500°, which appears in the equation for the 
viscosity, is also taken from reference 10, but 
Edwards’ " experimental value, C2370°, for the 
range 15°C to 184°C yields values of Rr agreeing 
somewhat better with the observations. As for 
the Lennard-Jones 9,5 model, the prediction that 
Rr=0 for a temperature of 1.54 T. or 625°K 


10M. Trantz and R. Heberling, Ann. d. Physik 10, 155 


(1931). 
1 R, S. Edwards and B. Worswick, Proc. Phys. Soc. 38, 


16 (1925). 


is.in marked disagreement with the observed 
T,=293°K for the reversal point. 

To describe adequately the intermolecular 
forces between dipole molecules such as ammonia 
no simple model would suffice. For there are in 
addition to the attractive forces arising from the 
dispersion effect which are common to all] 
molecules and which vary as 1/R’, R being the 
distance between interacting molecules, (1) the 
dipole induction effect (fa1/R"), (2) the align- 
ment effect which depends on the rotational 
quantum number and hence on the temperature 
distribution of the molecules (fa+1/R’ for large 
R but a1/R‘ for close approach), (3) resonance 
forces due to the possibility of exchange of a 
quantum of radiation if the interacting molecules 
have rotational quantum numbers differing by 
+1(fa+1/R‘), and (4) the first-order forces 
between symmetrical top molecules carrying a 
dipole along their figure axis (fa+1/R‘). It has 
been shown by Margenau and Warren’ that for 
NH; (u=1.5 X10-"* e.s.u.) these first-order forces 
(4) are greater than the second-order effect for 
linear dipoles over the main part of the range of 
the latter. And since for the repulsive part of 
these first-order forces the index is less than the 
critical value 5 for the Maxwellian case, for 
which the thermal diffusion effect is zero, they 
would have the tendency to make the thermal 
diffusion constant negative. Furthermore, all of 
these dipole-dipole forces (1)—(4) are temperature 
dependent, becoming more effective as the tem- 
perature is lowered because of lower molecular 
speeds and the concentration of the molecules 
into the few lower rotational quantum states. 

Qualitatively it is thus clear why the thermal 
diffusion constant of ammonia becomes negative 
and why the molecule “‘softens”’ at such a rapid 
rate with decreasing temperature. More such 
measurements might yield considerable informa- 
tion about certain long-range van der Waals 
forces. It would be of interest, for example, to 
have similar data for a linear dipole molecule 
such as HCN for which the dipole moment is 
large but for which the symmetrical top effect 
is, missing. 

We wish to thank Professor L. Onsager for 
stimulating discussions about thermal diffusion. 
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LTHOUGH the problem of the propagation 

of electromagnetic waves along a plane 
separating two dissimilar isotropic homogeneous 
media can be at present considered completely 
solved through the work by Sommerfeld and his 
followers, the more difficult case of three or more 
different media seems not yet to have been 
successfully attacked. Such an_ investigation 
would, however, be of considerable interest, both 
theoretical and practical. In particular, the 
theory of the coastal refraction, i.e., of the 
peculiar way in which electromagnetic waves are 
reflected and refracted when passing from the 
sea to the shore, relates to the case of three 
different media (air, water, soil, see Fig. 1). 

We appreciate the difficulty of this last 
problem when we recall that even the much 
simpler problem of the diffraction of electro- 
magnetic waves by a simple wedge of finite 
conductivity has not yet been solved exactly. 
We therefore cannot expect that the much more 
difficult problem of the coastal refraction, which 
in its simplest form (when the boundaries 
OA, OB, OC of the different media are planes) 
corresponds to the case of two such wedges with 
different properties, can be solved exactly.! 

In such circumstances it seems to be advisable 
to seek for some approximate method of approach 
to the problem under consideration. Such a 
method can be based on the well-known fact 
that an electromagnetic field varying with the 
time can penetrate but little into a sufficiently 
good conductor. This leads to the possibility of 
formulating a certain approximate boundary 
condition at the surface of separation of two 
media, one of which has a much higher con- 
ductivity than the other. 

1 It might be of interest to note that the corresponding 
electrostatic and related problems allow of an exact 
solution (in quadratures) in the case of an arbitrary 
number of such wedges, the vertices of which are in contact 
with each other. See G. Griinberg, Tech. Phys. U.S.S.R. 
5, 525 (1938); G. Griinberg, J. Phys. U.S.S.R. 2, 213 
(1940); G. Griinberg, J. Phys. U.S.S.R. 3, 401 (1940); 


G. Griinberg and I. Zuckermann, J. Tech. Phys. (russ.) 
12, 36 (1942). 
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Let us derive this boundary condition as 
follows: Consider the propagation of electro- 
magnetic waves along a plane separating two 
dissimilar homogeneous isotropic media and let 
x, y, z be a set of rectangular coordinates whose 
origin lies in the plane of separation and whose 
z axis is directed upwards. Then the equation 


=— (1) 
ot 


holds where E;'® denotes one of the components 
of the electric vector in the second medium. 
Confining ourselves to the case of simple har- 
monic oscillations and putting consequently 
Ej/®=E;e-“', we find Eq. (1) assumes the 
form 


+ + 
Ox? dy? 02" 


=—k,*E;, (2) 


where 
(€2w +4 ot) pow 


2 

If the medium (2) possesses a sufficiently high 
conductivity, the electric field must rapidly 
decrease in the direction of the negative z axis, 
its rate of decrease being greater for greater 
conductivity o2 of the medium. On the other 
hand in the directions x and y, i.e., along the 
surface of separation, the field varies relatively 
slowly.? It follows, therefore, that under the 


(1) (air) 


(soil) (3) 

| (2) | A 
Fic. 1. 


2 E, and E, vary continuously on passing through the 
surface of separation, and E,™ and E,® are proportional 
to each other [see Eq. (5) ]. Therefore, their variation 
along the earth’s surface is in the main determined by the 
wave-length \2xc/w in the air. 
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conditions stated the derivatives 0*E;®/dx? and 
@E;/dy*? will be small in comparison with 
@E; /dz*. Neglecting them, therefore, in Eq. (2) 
we get approximately 
hence 
E;™=A exp exp (tkez). 


For one .must have E;®=0; conse- 


Z 


ATUL LT y 


Fic. 2. 


quently if k2 is so chosen that R(ike) <0, then 
B=0 and 


E;=A exp (—tk22) 
= (E;),-0-exp (3) 
This solution has been derived under the 
assumption that and are 
relatively small in comparison with 0?E,;® /dz?. 
But insofar as orders of magnitude alone are 
concerned, one has 
= (2x/d)2E;, 
/dy? (24/d)2E,, 


where A denotes the wave-length in air. Hence, 
using Eq. (3) we find that the condition for the 


validity of Eq. (3) is 


( (4) 
d| 


We are now in a position to formulate the 
boundary condition in question, which we shall 
write down only for the vertical component £, 
of the field. We start with the two exact bound- 


ary conditions 
(exw +4201) = (5) 
dE, /dz = dE, /dz, (6) 
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which hold at the plane z=0 of separation, the 
first of them expressing, as usual, the continuity 
of the normal component of the total current 
and the second one being an immediate conse- 
quence of the fact that in both media div E=9 
and at the surface of separation itself JE,“ /dx 
=0E,”/dx and dE, /dy=dE,/dy.* Inserting 
now into Eq. (6) the value /dz),_, 
= —tk,(E,),-0 from Eq. (4) and using Eq. (5), 
we at once get: 


, 
oz z=0 z=0 
, (7) 
(exw +4201) 


which is the condition required. Its validity is 
connected with the degree of accuracy of Eq. (3), 
which in its turn depends on the smallness of the 
quantity (27/|k2|)?. If the displacement current 
in the ground can be neglected, inequality (4) 
assumes the form: 


1/2nd\? 
-(—) 
2\ A 

where d=c/(2mp20qw)!. In this case 


Qs (1 


and E,® =(E,),.9-exp [(i—1)z/d], so that d 
determines the thickness of the layer of the 
ground within which the field differs sensibly 
from zero. For example, when A=3 X10‘ cm, 
o2=5 X10", w2=1, the ratio of the displacement 
current to the conduction current is about 0.01 
and d220.02\, For greater 
wave-lengths and the same conductivity of the 
soil, conditions for the application of Eq. (7) are 
even more favorable. It should also be noted 


that for eg= © Eq. (7) gives 
(E, (8) 


E,Y=E,® and E,=E,® at all points of the 
plane z=0. It might be of interest to note that if the 
surface of separation were not plane, then introducing 
appropriate curvilinear coordinates and transforming 
equation div E=0 into them. we could find, taking into 
account the continuity of the tangential components of 
the field, the corresponding generalization Eq. (6). 


A similar condition holds for the normal component H, of 
the magnetic field. 


1 nee 


cto 
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and this is readily seen to be the exact boundary 
condition for perfectly conducting ground. 


In order to apply the results obtained to the 
solution of the problem of the coastal refraction, 
let us consider the following idealization of the 
problem in question: Given an island (f) on a 
perfectly conducting sea, the surface of the island 
being strictly horizontal and at sea level. At the 
surface of the island condition (7) holds while at 
the sea surface (0E,/2),.0.=0. Find the electro- 
magnetic field produced in the upper half-space 
(air, z>0, see Fig. 2) by a given distribution of 
oscillators placed within this half-space. 

Such an idealization evidently can be used 
only when the surface of the island is sufficiently 
flat, horizontal, and does not protrude much 
from the sea. Moreover, the boundary condition 
(7) certainly cannot be correct in the immediate 
neighborhood of the coast line. Still it seems one 
may be allowed to use it even there in conse- 
quence of the supposed smallness of the ratio 
d/\, which shows that already at a small distance 
off the coast line (of the order of a small fraction 
of the wave-length A) condition (7) can be 
safely applied.‘ 

We are now in a position to derive the funda- 
mental integral equation of the problem. To 
that end let us apply Green’s theorem 


ae 
f f (9) 
(v) ce) \ On on 


to a hemisphere of radius R-~, lying in the 
upper half-space, its center coinciding with the 
origin of coordinates (Fig. 3). 

Putting here 


g=E,, w=exp (ikir)/r, 


excluding the point M(x, y, ) (where r=0) from 
the range of integration by means of a small 
sphere of radius p—0, and taking into account 


‘The thickness of the layer of the ground, beneath 
which the field is peety zero, is of the order of d. 
It is, therefore, to be expected that at a distance of the 
order of a few d, say 3d or 4d, condition (7) holds suffi- 
ciently well. If, as in the case considered above, d/rx< 0.02, 
then 3d¢ 0.06, which proves our assertion. 
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that AE,+k’E,= —4nf,, where f, is a given 
function determined by the distribution of the 
“primary” currents in space, we easily get 


( 


A 


f f. exp (tkir) 


r 


exp (ikyr) ) 
On r 


where the subscript M denotes that the corre- 
sponding quantity pertains to the point M. The 
integral on the right in (10) is extended over the 
whole plane AB,‘ and E£,° is the “‘primary”’ field, 
which would be produced by the given sources 
in the absence of the surface of separation, i.e., 
if the air were to fill in the whole space. On 
applying formula (10) to a point M which lies 
on the earth’s surface and thus makes Z—-0, we 
readily find 


/exp (ikir 
(AB) 


“an r 


4 


Fic. 3. 


When we introduce this into Eq. (10) and use 
the boundary conditions (7) and (8), we obtain 


tk 
=2 (E,°),-0+ af (1 1) 


This is the integral equation in question. From 
the derivation just given it follows that Eq. (11) 
applies just as well to the case of several islands 
present, if (f) denotes their total area. It should 


5 The surface integral extended over the outer hemi- 
sphere of radius R->« disappears (at least if we take into 
account the conductivity of the air). 
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be noted that for perfectly conducting ground, 
i.e., when a=0, Eq. (11) gives 
The normal component of the actual field at the 
surface of a perfectly conducting ground has, 
therefore, only in the absence of the ground, 
double the value of the normal component of 
the “primary” field due to the given sources. This 
corresponds evidently to the usual “reflection” 
rule for a perfectly conducting plane. 


3. 


Let us apply Eq. (11) to the case when the 
region (f) occupied by the dry land is a half- 
plane (Fig. 4). It suffices to consider the case 
of a primary field (Z,°),.o having the form 
(E,°),.0=f(x)e**¥ [where Y is the direction of the 
coast line (Fig. 4) and s is a constant], for an 
arbitrary primary field can always be resolved 
into such components by means of a Fourier 
integral. 

Putting consequently (E.,)..0o=F(x)e"*” we 
obtain Eq. (11). 


F(x) =2f(x) +a f F(&)dé 


exp J 
dn, (12) 
or using the formula 
exp 
be 
J], (13) 


where H,™ denotes the Hankel function of the 
first kind and of zero order, we find 


F(x) =2f(x) +a f (14) 


m= 


This equation allows for an exact solution 
(in quadratures).6 This shows that in those 
cases when the approximate boundary condition 
(7) proves to be sufficiently accurate, the solution 


®* By a method of Hopf and Wiener; see Titchmarsh 
Introduction to the Theory of Fourier Integrals (Oxford, 
1937), p. 339, and generalized by E. Reissner, J. Math. 
and Phys. [5], 20, 219 (August 1941). 


G. GRUNBERG 


Z P 
y 
Fal 

Fic. 4. 


of the problem of the coastal refraction can be 
given in quadratures for an arbitrary primary 
field. We intend to give this general solution in a 
subsequent paper as soon as possible. Here we 
shall confine ourselves to the consideration of the 
case of a plane wave propagated from the sea, 
its direction making an angle @ with the X axis. 
Then f(x) and if we put F(x) =2e™Y(x) 
in Eq. (14) we obtain: 


0 
=1+rai| f (x — 8) 
0 


+f FL 2) (15) 


Instead of solving this equation exactly, basing 
on the remark that ¥(x) is a function which in 
general varies but slowly in comparison to e™, 
we shall give here its approximate solution.’ 
This enables us to put approximately® 


7 For perfectly conducting ground y(x) =1. If the ground 
is not perfectly conducting but possesses a sufficiently 
high conductivity, then at a distance of the order of a few 
wave-lengths off the coast line, the wave must be propa- 
gated much in the same way as for perfectly conducting 
ground, its amplitude and phase varying but slowly in the 
direction of propagation. This proves our assertion if only 
@ is not near 90°, i.e., if the incident wave (primary) is 
not a glancing one. It should also be noted that in the 
neighborhood of the coast line ¥(x) cannot differ appreci- 
ably from unity, the difference being of the order of the 
(small) quantity a/m. 

8 For e™*H,)(mz) is a nearly twice as rapidly varying 
function as e*"* and, moreover, rapidly ee is 
can be seen from the asymptotic expression Ho“(mz) 
es (2/a4mz)te*™=-*"), which gives a good approximation at 
arguments of the order of a few units and which shows 
that e*™ (mz) ~ In contradistinction, 
it is not possible to use this kind of argument for a similar 
simplification of the first member on the right of Eq. (15), 
for asymptotically 

e~im (mz) (2/amz) 


which is not at all a rapidly varying function. 
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f “gin Ho Em(E—x) 


242) f "ime (16) 


the last result being a consequence, for x= ~, 
of the formula 

Zz 1 2 
Fo (mu)du -—| +—+mxe*im= 

0 m T 


X [Ao (mx) (mx) ), 


which may be easily verified by differentiation.® 

According to Eq. (16) the second member on 
the right-hand side of Eq. (15) is approximately 
equal to (2at/m)y(x). This is usually not more 
than a small fraction of the member (x), which 
stands on the left of (15). For instance when 
4=3X104 cm, o2=5X10", we=1, we find that 
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d/A=2 X10 and m= (w cos cos 8, 
2ai/m=(d/d cos 0)(ti—1), which, for @< 60°, 
yields cos 4X10-*. It is therefore clear 
that even a relatively high error involved in the 
approximate evaluation of the integral (16) can- 
not affect considerably the determination of (x) 
from Eq. (15) if, instead of this last equation, 
its approximate form 


f 
0 


U=1—(2ai/m), 


K(z) =e-*™* (mz), 


is used. This is an integral equation of Volterra 
type, which can be solved readily by means of 
Laplace transforms. Its solution can be shown 
to be 


o+in 


B= —-rai/Um, 


where the path of integration is a straight line, 
lying in the right half-plane z, parallel to the 
imaginary axis; and the values of the radical 
and the logarithm are uniquely determined by 
the condition that they both are real and positive 
on the positive part of the real axis. 

A detailed investigation of the solution (19) 
will be given in a paper which is to appear in the 
Journal of Physics of the U.S.S.R. (under the 


* This formula is a special case of a series of general 
formulae relating to Bessel functions of any kind, which 
we intend to publish elsewhere. 


{1—(2é/m) Ig 


(z*+1)! 


o>0, 


same title as above). It should only be noted 
that at a sufficiently great distance from the 
coastline formula (19) gives approximately : 


Al(cos 6)! 


(x2), (20) 


which in particular shows that at such distances 
the direction of propagation of the “‘refracted”’ 
wave coincides asymptotically with that of the 
primary one, its amplitude being a slowly 
varying function of x in accordance with the 
assumption made above. 
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Calculations of the electronic energy bands in body- 
centered iron by the Wigner-Seitz-Slater method are re- 
ported. There are found to be two filled and four partially- 
filled bands. The lowest two filled bands are responsible for 
the cohesion; the next three partially-filled bands for the 
ferromagnetism; and the highest band is responsible for the 
electrical conduction. The width of the occupied levels is 
about 0.6 Rydberg but the density of states is high for only 
about 0.4 Rydberg in good agreement with the data ob- 
tained from x-ray emission bands. The density of states 
is calculated as a function of the energy and compared with 


(Received December 19, 1942) 


similar results obtained by Greene for face-centered iron, 
The density of states at the highest occupied energy 
level is 17 electrons/atom-Rydberg compared with 1} 
electrons/atom-Rydberg for face-centered iron. This higher 
density of states accounts for the fact that body-centered 
iron is ferromagnetic, but face-centered iron is not. It also 
indicates a greater electronic specific heat for body- 
centered iron than for face-centered iron. This higher 
electronic specific heat is responsible for the high tempera- 
ture change from a face-centered to a body-centered 


structure. 


INTRODUCTION 


HERE is considerable scientific interest in 
metallic iron because it is one of the few 
elements which show ferromagnetism and be- 
cause of the phase changes from body-centered to 
face-centered and back to body-centered as the 
temperature is raised. 

The method of calculating the electronic 
energy levels reported in this paper is that de- 
veloped by Wigner and Seitz! and extended by 
Slater? and others.* The essential feature of this 
method. is that the wave function within the 
polyhedron surrounding any one atom is ex- 
pressed as a sum of terms (fourteen for the body- 
centered lattice) each of which consists of a 
spherical harmonic multiplied by the appropriate 
radial function. The radial functions are found by 
a modified Hartree procedure and the arbitrary 
coefficients in the linear combination are chosen 
so that the wave functions satisfy the Bloch‘ 
periodicity condition and join smoothly in 
passing from one atom to its eight nearest and six 
next-nearest neighbors. It turns out that these 
conditions impose a determinantal relationship 
between the energy of an electron and the com- 
anscomnadl of the momentum or wave number 


Deceased June 1, 1942. 
1E. Wigner and F. Seitz, Phys. Rev. 43, 804 (1933); 
Phys. Rev. 46, 509 (1934). 
. Slater, Phys. Rev. 45, 794 (1934). 
. I. Chodorow and M. F. Manning, Phys. Rev. 52, 
731 (1937). For a list of other applications of this method, 
see M. F. Manning and M. I i. Chodecon, Phys. Rev. 56, 
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Bloch, Zeits. f. Physik 52, 555 (1928). 
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vector k. For a general direction of k the deter- 
minantal relations become very complicated, but 
for a number of symmetrical directions, simpler 
relations are obtained. For the body-centered 
lattice the results of the expansion of the de- 
terminant for the special directions are given in 
reference 3. The directions of k for which solu- 
tions have been obtained are indicated in Fig. 1, 
which is a drawing of the first Brillouin zone 
corresponding to a body-centered structure. 


RESULTS 


In the calculations the first step was the de- 
termination of an effective field for the calcula- 
tion of the radial factors of the different wave 


Fic. 1.’The first Brillouin zone for the 
body-centered lattice. 


5 The notation is that of L. Bouckeart, R. Smoluchowski, 
and E. Wigner, Phys. Rev. 50, 58 (1936). 
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o7 


wo a 


Fic. 2. Energy ¢ vs. momentum X for the different directions of propagation. 


functions. In determining the effective field for 
the metal, the charge-densities for 1s to 3p 
electrons were assumed the same as found for the 
iron atom by Manning and Goldberg.* The 
charge-densities for the 4s and 3d were re- 
normalized to the volume of the atomic sphere in 
metallic iron and other readjustments made. 
Further details are given in Appendix I. 

By use of the field so obtained, the radial wave 
functions were found by numerical integration as 


TasLe I. Values of ¢ for which certain important 
functions vanish. The notation is that of Chodorow and 
Manning (reference 3). 


Function 


351'go' 
S; 


dy’ 
D,B;+B, 
ds 

d, 


TasLe II. Maximum and minimum values of « for 
each _ and values of ¢ at the symmetry points (refer- 
ence 3). 


0.365 
—0.160 
0.365 
—0.015 
—0.015 
—0.160 


0.470 
—0.160 


0.470 
—0.160 
0.470 
—0.015 
— 0.008 
—0.160 


*M. F. Manning and L. Gold , Phys. Rev. 53, 662 
(1938). 


in the Hartree self-consistent field method for 
atoms. The values of the functions and their 
radial derivatives at the half-distances between 
nearest neighbors and between next-nearest 
neighbors were then determined by standard 
methods of interpolation and numerical differ- 
entiation. (For iron the lattice parameter at room 
temperature is 2.86A. The corresponding distance 
between nearest neighbors is 2.345 Bohr units 
and the distance between next-nearest neighbors 
is 2.71 Bohr units. The corresponding values of 
p(=log, 1000r) are 7.760 and 7.904.) 

The values of the functions and derivatives are 
then used to calculate the coefficients in the 


transcendental equations relating E and k given 


in reference 3. The procedure in solving these 
equations was very similar to that used in the 
calculations for metallic tungsten.? A few com- 
ments about slight modifications in the procedure 
are given in Appendix II. The results of the 
calculation of the energy as a function of the 
momentum are given in Fig. 2. It does not seem 
necessary to tabulate the numerical values from 
which the figures were plotted, but Table I gives 
the energies at which certain key combinations of 
functions and derivatives vanish, and Table II 
gives the maximum and minimum energies for 
each band and the energies at a number of special 
points in the Brillouin zone. In these tables and 
elsewhere in this paper ¢ is the negative of the 
energy expressed in Rydbergs. 

Although energy-momentum relations are ob- 


7M. F. Manning and M. I. Chodorow, Phys. Rev. 56, 
787 (1939). 
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Fic. 3. n(e) vs. ¢ for the separate bands. m(e) =number of 


tained directly for only a limited number of lines 
in k space, very plausible constant energy con- 
tours based on these known lines can be drawn 
for the principal planes in k space. By considering 
the contours in these planes as known, energy 
contours for a number of planes K,=0, 0.1, 0.2, 
etc., can be interpolated. Each of the planes 
K,=constant can be divided into a region in 
which the energy is greater than a particular 
.value and a region in which the energy is less 
than this particular value. For each plane the 
area of the region with energy less than a par- 
ticular energy was obtained by using a planimeter. 
This procedure was repeated for all of the planes 
and by numerical integration the volume in k 
space with energy less than a particular energy 
was determined. Since this volume is propor- 
tional to the number of electrons with energy less 
than the energy considered, the corresponding 
number of electrons can be found by using the 
normalizing relation that there are two electrons 
per atom in each energy band. In the present 
paper this quantity is designated by N(e). The 
values of N(e) for each of the energy bands and 
for the sum of the six bands are given in Table ITI. 
Since there are eight valence electrons per atom 
in iron, the highest occupied level at absolute 
zero in non-magnetic iron occurs at e=0.105. 
From the values of N(e), the values of 
n(e)(=—dN/de) can be found by numerical 
differentiation. (n(e)de is the number of energy 
levels per atom with energy between ¢« and 
¢+de.) Since numerical differentiation is a process 


states per unit energy range. 


in which considerable accuracy is lost, the values 
of n(e) found by differentiation were adjusted so 


that 


as found by a planimeter checked the corre- 
sponding increment in N(e). The resulting values 
of n(e) for the separate bands and for the total 
are given in Figs. 3 and 4. 


TABLE III. Values of N(e) as a function of «. 


€ I II III IV ¥ VI Total 

0.70 
0.003 

0.65 0.017 0.017 
0.032 

0.60 0.058 0.088 
0.090 

0.55 0.122 0.122 
0.199 

0.50 0.275 0.275 
0.456 

0.45 0.623 0.306 0.025 0.003 0.957 
1.003 0.599 

0.40 1.328 0.902 0.118 0.025 2.373 
1.614 1.138 

0.35 1.970 1.494 0.288 0.063 0.033 3.848 
2.000 1.651 

0.30 2.000 1.802 0.460 0.156 0.103 0.034 4.555 

1. 

0.25 2.000 1.998 0.968 0.336 0.170 0.078 5.550 

0.20 2.000 1.477 0.541 0.245 0.121 6.384 

0.15 2.000 1.895 0.904 0.372 0.188 7.359 

0.10 1.939 1.302 0.552 0.255 8.048 

0.05 1.967 1.677 0.945 0.406 8.995 

0.00 1.982 1.937 1.405 0.468 9.792 

—0.05 1.992 1.971 1.747 0.589 10.299 

—0.10 1.997 1.993 1.904 0.816 10.710 

—0.15 2.000 1.998 1.979 0.882 10.859 
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Fic. 5. The average Fermi energy é as a function of the 
number of electrons per atom. The upper curve is for 
bod y-centered iron, the lower one for face-centered. 


From the results depicted in these figures, 
some general remarks about the electronic struc- 
ture of metallic iron can be made. Bands I and II 
are completely filled in the metal and hence do 
not contribute to either the electrical con- 
ductivity or the ferromagnetic properties. Band 
VI is a wide band with a considerable admixture 
of s functions and probably accounts for most of 


the conductivity. Bands III, IV, and V are 
narrow bands with a high density of states and 
are the important bands in connection with 
ferromagnetism. 

This model for the electronic structure of 
metallic iron combines many of the features of the 
“d” band model suggested by Mott® and the 
homopolar group of levels suggested by Pauling.® 

The general features of the energy bands for 
body-centered iron are very similar to those 
found for tungsten.’ The only important differ- 
ence is that for tungsten the energy for which 
35,'g2'+2se'g;'=0 lies above the energies for 
which d,’ and d,’ are zero, while for iron this con- 
dition is fulfilled for an energy lower than the 
energies for which d,’ and d,’ vanish. (See Fig. 6.) 
The general results are also quite similar to those 
obtained by Greene’ for face-centered iron. 
However, the lowest band in face-centered iron 
begins at the energy for which s’=0 and all five 
of the d bands start at the energy for which d’=0. 
(Actually this is due to an approximation in the 
theory as Bethe" and Chodorow” have pointed 
out that even at the center of the momentum cell 
the d level should split into a group of two levels 
and a group of three levels.) 


8. F. Mott, Proc. Phys. Soc. 47, 571 (1935). 


*L. Pauling, Phys. Rev. 54, 899 (1938). 
(1948) B. Greene and M. F. Manning, Phys. Rev. 63, 203 


1H. Bethe, Ann. d. Physik 87, 55 (1928). 
2M. I. Chodorow, Phys. Rev. 55, 675 (1939) and 
M.I.T. thesis (1936). 
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APPLICATIONS 
Energy Relations 


A quantitative calculation of the binding 
energy of body-centered iron and its stability 
relative to the face-centered structure would re- 
quire a consideration of the exchange and corre- 
lation energies which is not possible by the 
method used here. It is, however, interesting to 
compute the average Fermi energy as a function 
of the total number of electrons per atom and 
compare with the same results obtained for the 
face-centered structure. 

The average Fermi energy is defined by the 
equation : 


tm 
N{én)(?—«0) = — f n(e(e—erdde, (1) 


+ 2 5; 


Fic. 6. The energy for which certain functions vanish vs. log, rf. 
r is the interatomic distance in Bohr units. 


where € refers to the lowest level and e,, to the 
highest level occupied at absolute zero. This 
integral could be evaluated from the information 
in Fig. 4, but since n(e) is obtained only as the 
derivative of another function, it is preferable to 
integrate Eq. (1) by parts to obtain: 


f (2) 
N(€m) to 


Although the results calculated are based on 
iron, they should be approximately valid for 
neighboring elements in the same row of the 
periodic table if the appropriate number of 
electrons per atom is used. Similar calculations 
have been made by Greene" for the face-centered 
structure and the results are plotted in Fig. 5. In 
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the neighborhood of iron the difference in mean 
Fermi energy is about 0.02 Rydberg/atom. 
Since the extra magnetic energy of the body- 
centered lattice is of this same order of magnitude 
and there are a number of other corrections which 
may be different for the two structures, it is easy 
to see why iron might show phase changes. In the 
neighborhood of ten electrons per atom the 
average Fermi energy for the body-centered 
structure rises much more rapidly than for the 
face-centered structure, thus indicating that for 
nickel the face-centered structure should be more 
stable. This gives no information about cobalt, 
although here the face-centered structure is 
favored over the body-centered structure by the 
magnetic energy. 

The considerable concentration of states near 
the bottom of the d band has been found for all 
of the transition elements which have been 
studied'* *}7 and may be considered responsible 
for their relatively high binding energies. This 
possibility has been emphasized by Manning and 
Chodorow’ and by Seitz.'> 

The behavior of the lowest energy in each type 
of band is qualitatively similar to that found by 
Greene'® for face-centered iron. However, the 
method used by Greene does not show the 
splitting of the d band at k=0, and the lowest 
band for the face-centered lattice begins at the 
energy for which s’=0 while for the body- 
centered lattice the lowest band begins at the 
energy for 3s;'g2’+2se'g,'=0. It turns out that 
for the portion of the s’=0 curve inside of the 
minimum (it seems probable that this is always 
the important part for transition metals) the 
second condition corresponds to a curve shifted 
slightly toward smaller values of interatomic 
distance but otherwise similar. The behavior of the 
energies for which d;’=0 and 3s,’go’+2se'g:'=0 
as a function of interatomic distance is shown in 
Fig. 6. It is interesting to note that for tungsten 
with its smaller number of d electrons and a 
recently-filled f shell, the normal interatomic 
distance corresponds to the 3s,’g2’+2se'g,' curve, 
being higher than the d’ curves. 


%H. M. Krutter, Phys. Rev. 48, 664 (1935). 

4 J.C. Slater, Phys. Rev. 49, 537 (1936). 

°F. Seitz, Modern Theory of Solids (McGraw-Hill, New 
York, 1940), p. 430 ff. 

°F. Seitz and R. P. Johnson, J. App. Phys. 8, 186 (1937). 


Ferromagnetism 


The relation of the ferromagnetic behavior of 
iron, cobalt, and nickel to the structure of their 
electronic energy bands has been discussed by 
Slater in a series of papers.*!7'!8 There are 
unsurmounted difficulties in the formulation of a 
complete theory of ferromagnetism, but certain 
general features can be discussed in terms of the 
present model. 

In the non-magnetic state there are equal 
numbers of electrons with each spin, whereas in 
the magnetic state there are more electrons with 
positive than with negative spin. In the mag- 
netized state the Fermi energy is of course 
greater, because in order to reverse the spins it is 
necessary to transfer electrons to higher vacant 
levels of the opposite spin. For substances which 
are ferromagnetic this increase in energy is more 
than compensated by the increase in the number 
of exchange integrals between electrons with 
parallel spin. Slater'* showed that this change in 
exchange energy is proportional to the square of 
the number of uncompensated spins. The change 
in Fermi energy is also proportional to the square 
of the number of uncompensated spins if m(e) can 
be regarded as constant. The actual relation is 


Aes =u? /2n(em), (3) 


where yu is the number of uncompensated spins 
expressed in Bohr magnetons per atom, (€,) is 
the number of energy levels (both spins) per 
Rydberg unit of energy at the highest occupied 
level for the non-magnetic state, and Ae; is the 
change in Fermi energy upon reversing 4/2 spins. 
For ferromagnetism Ae,, the change in exchange 
energy accompanying the change to the magnetic 
state, must be larger numerically than Aey. Hence, 
it follows that only those substances are ferro- 
magnetic for which n(e) is large. In this connec- 
tion it is interesting to note that Greene's results 
for face-centered iron indicate a much smaller 
density of states near ¢,, than reported here for 
body-centered iron; therefore, it is quite under- 
standable that face-centered iron should not be 
ferromagnetic.'® 


i“ . C. Slater, Phys. Rev. 52, 198 (1937). 
above the Curie 
temperatures, no direct evidence about its ferromagnetism 


sad . C. Slater, Phys. Rev. 49, 931 (1936). 


ince face-centered iron is stable only at temperatures 
int and cannot be retained at lower 
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If the number of states per unit energy range is 
not constant, the change in Fermi energy ac- 
companying the change to the magnetic state is 
no longer given by Eq. (3) but can be found by 
using Eq. (1) for each direction of spin. In 
computing the change in Fermi energy, Bands I 
and II can be neglected because they are com- 
pletely filled. Also, it is usually considered that 
Band VI can be neglected because of the low 
density of states and because the s and p func- 
tions which predominate in this band are concen- 
trated further from the nucleus than the d 
functions which predominate in the other bands 
and hence contribute less to the exchange energy. 
The actual result of the computation of Fermi 
energy as a function of » showed that the Fermi 
energy increased slightly more rapidly than pro- 
portional to the square of yu, but the difference 
from strict proportionality was not enough 
greater than the computational error to be 
significant. 

For iron there are eight valence electrons per 
atom; four in the two lowest filled bands, and 
0.25 electron in Band VI. There are thus 3.75 
electrons per atom in the three bands which are 
assumed to contribute to the ferromagnetism. 
With only the three bands contributing to the 
ferromagnetism, the saturation magnetic mo- 
ment is determined by one of three conditions :*° 


1. The number of negative spins per atom 
cannot be negative. 

2. The number of positive spins per atom 
cannot be greater than three. 

3. If Ae, increases less rapidly with yu? than Ae,, 
the difference may vanish at some value of 
u less than that determined by either of the 
above conditions. 


There does not seem to be any element or alloy 
in the iron group for which the first condition is 
the determining one. The second would give a 
saturation moment of 3.00—0.75=2.25 Bohr 
magnetons per atom. The experimental value is 
about 2.20 or 2.22 magnetons per atom. This 
small difference cannot be regarded as conclusive 


can be obtained. However, any reasonable extrapolation 
of the high-temperature susceptibility of 
face-centered iron and the Curie points in the face- 
centered portion of the Fe-Ni alloy system indicate that 
face-centered iron would not be ferromagnetic. 

”F. Seitz, 1eference 15, page 430. 
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Fic. 7. Magnetic moments of Fe-Ni-Co alloys as a function 
of composition. (After Borelius.) 


evidence against the assumption that the second 
condition is the important one. In order to say 
definitely about the third, it would be necessary 
to have more definite information about the 
variation of the Fermi and exchange energies 
with 

For the alloys of iron with manganese, chro- 
mium, and vanadium the saturation moments 
decrease with addition of alloying element*! and 
the Curie point also decreases except for a slight 
initial rise in one or two cases. Since, according to 
Fig. 4, the m(e,) values increase with a decrease 
in the number of electrons per atom, the most 
likely explanation of the decreasing saturation 
moments and Curie temperatures is that they are 
due to a decrease in the exchange energy gained 
upon reversal of spin. This decrease in exchange 
energy is due to the fact that for atoms with 
smaller nuclear charge, the d functions are con- 
centrated further from the nucleus and hence 
contribute a smaller exchange integral.” 

For the Fe-Ni (Fig. 7) system the saturation 
moment increases slightly with the addition of 
Ni, and then decreases to the phase boundary* 
where the saturation moment is about 2.2 
magnetons. Since the corresponding number of 
electrons per atom is 8.4, it is necessary to assign 


21 Data sant by J. C. Slater, reference 17, page 935. 


later, reference 17, page 934. 


27.C. 
der Metallphysik (Leipzig, 


%G. Borelius, Handbuch 
1935), Vol. I, p. 291. 
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0.6 electron per atom in order to account for the 
observed moment. 

For cobalt-iron alloys the saturation moment 
reaches a maximum of 2.4 magnetons at 8.4 
electrons per atom and then decreases to 2.0 
magnetons at the phase boundary which occurs 
at 8.8 electrons per atom. For both cases the 
corresponding number of electrons per atom in 
Band VI is equal to 0.80. 

The number of electrons per atom in Band VI 
for nickel and cobalt alloys is much greater than 
the values extrapolated from the present calcu- 
lations. Two possibilities of accounting for this 
difference can be suggested. It may be that for 
these elements the energy at the top of Band VI 
is much lower than for iron and the density of 
states is so much higher, that the number of 
electrons in Band VI can be as large as 0.80. (It 
is also possible that there is such an effect in iron 
and that the saturation moment is determined by 
a balance of exchange and Fermi terms.) It is 
also conceivable that Band VI contributes to the 
ferromagnetism. It would not be ferromagnetic 
by itself, but a net spin for this band might be 
produced by a coupling with the net spins in 
other bands. 


Electronic Specific Heat** 


The electronic contribution to the specific heat 
is given by the formula”®** 


(Co)e= 0.209 X 10-*n(en) RT, (4) 


where (C,), is in cal./g-atom C°, n(e,) is the 
number of energy levels per atom per Rydberg 
unit of energy at the top of the Fermi distribu- 
tion, R is the universal gas constant, and T is the 


“In a previous paper (reference 7) Manning and 
Chodorow reported calculations of the electronic specific 
heat of metallic tungsten and tantalum. For both of these 
metals the computed electronic specific heat checked the 
excess of C, above 6.0 cal./g-atom C°, at high —_% 
tures within the accuracy of the C,—C, correction. There 
were not then and are not yet available any low tempera- 
ture data for tungsten. For tantalum the electronic 
specific heat had been computed from data on the effect 
of a magnetic field on the superconducting transition 
temperature and a result of approximately 27 X 10-*T had 
been obtained. More recent data by Keesom and Desirant 
EN Roy. Soc. Ams. 42, 536 (1939) ] and by Mendelssohn 

Nature 148, 316 (1941)] yield a value 14X10~“T com- 
pared with the calculated value of 6X 10~T. 

** E. C. Stoner, Proc. Roy. Soc. A154, 656 (1936). 

**H. Bethe, Handbuch der Physik (Julius Springer, 
Berlin), Vol. 24, Part 2, p. 430. ; 
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absolute temperature. For body-centered iron 
n(€m) is 17 for the non-magnetic state. Hence the 
corresponding electronic contribution to the 
specific heat is 6.9X10-‘T. This result is valid 
only at temperatures well above the Curie tem- 
perature. Any possibility of checking this against 
an excess of C, above the Dulong and Petit value 
at high temperatures is made impossible by the 
phase change. At higher temperatures where the 
body-centered phase is again stable, there are no 
satisfactory data for determining the dilatation 
correction. 

In one respect there is satisfactory agreement 
with experiment. As shown in detail in a later 
section the body-centered structure can become 
more stable than the face-centered structure at 
high temperatures only if the specific heat at 
constant pressure for this structure is higher in 
this region than the specific heat?’ for the face- 
centered structure. It is very probable (see 
reference 40) that the difference is due chiefly to 
electronic contributions and hence it is a quali- 
tative check on the theory that the value found 
for body-centered iron is about one and a half 
times the value found by Greene for face- 
centered iron. 

At low temperatures the distribution of elec- 
trons is different and the specific heat for each 
direction of spin must be computed separately. If 
the observed value of the saturation moment and 
the computed n(e) curves are used to determine 
the correct values of n(e,) to use in Eq. (4), the 
computed electronic specific heat is found to be 
4.6X10-T. Measurements by Duyckaerts** and 
Keesom and Kurrelmeyer®® at liquid helium 
temperatures yielded a term 12X10~T in the 
specific heat. 


X-Ray Emission Bands** 
The K£;-emission band of iron has been studied 


27]. B. Austin, Ind. and Eng. Chem. 24, 1225 (1932). 

28 G. Duyckaerts, Physica 6, 401 (1939). 
(1999), H. Keesom and B. Kurrelmeyer, Physica 6, 633 

*® Bearden and Snyder [Phys. Rev. 59, 162 (1941) 
have reported studies of x-ray emission and absorption 
metallic tungsten and compared their results with the 
nm(e) curves calculated by Manning and Chodorow. 


Bearden and Snyder compute from their data that the 
range of occupied levels in the metal is about 6.5 or 7.0 ev. 
In comparing with the calculated values they use 10.0 ev. 
as the calculated width. The present author's estimate of 
this width is 6.0 ev. This is in considerably better agree- 
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by Beeman and Friedman* and the La and L8é 
bands by Farineau.* Beeman and Friedman ob- 
tained a width of about 7 ev and Farineau 
obtained a width of about 5.5 ev. The width of 
the calculated n(e) curve is about 5 or 5.5 ev 
depending upon how much weight is given to the 
low energy portion of the m(e) curve. The agree- 
ment is probably better than should be expected, 
as the experimental lines are broadened by the 
finite width of the lower levels and the Auger 


effect. 


Electrical Conductivity 


At room temperatures the conductivity of iron 
is only about one-eighth of the conductivity of 


ment with the pee results, but the question of 
addition of line widths probably needs special consideration 
when the m(e) curves have such irregular shapes and 
when the relative proportions of the different types of 
wave function change with the energy. 

The reason for the discre ne 4 between the 6 ev and 
the 10 ev has been explained to the author by Dr. Snyder 
in personal conversation. It arose because the meaning of 
the term “atomic unit of energy” was not stated carefully 
enough. In both that paper and this one the unit of energy 
is the Rydberg unit of energy—13.54 ev—and not twice 
this quantity which is frequently referred to as an atomic 
unit. The author regrets any misunderstanding which this 
may have caused. 

31W. W. Beeman and H. Friedman, Phys. Rev. 56, 
392 (1939). 

% J. Farineau, Comptes rendus 206, 1011 (1938). 


Fic. 8. Specific heat of Fe vs. temperature. (After Austin.) 
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copper and only a small part of this difference can 
be accounted for by differences in Debye temper- 
ature and atomic volume. The lowest two bands 
can make no contribution to the conductivity 
and it is probably correct to assume that Bands 
III, IV, and V make little contribution to the 
conductivity because of the high effective mass of 
the electrons in these narrow bands. Most of the 
conductivity is therefore due to the electrons in 
Band VI. Since there is only 0.25 electron per 
atom in this band, a lower conductivity than for 
copper should be expected. Since the conductivity 
varies as the one-third power® of the number of 
conduction electrons per unit volume, the smaller 
number of conduction electrons in iron than in 
copper is responsible for an appreciable, but nota 
major part, of the difference in conductivity. 
Mott** has proposed that the low conductivity 
of the transition metals is due to the fact that the 
’ vacant states in the poorly-conducting bands 
increase the probability of scattering for the 
conduction electrons. This theory is able to ac- 
count qualitatively for a number of features such 
as the thermoelectric power and high tem- 


%F. Seitz, Modern Theory of Solids (McGraw-Hill, 


New York, 1940), p. 537. 
*N. F. Mott, Proc. Roy. Soc. A153, 699 (1936); A156, 


368 (1936). - 
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ture variation of resistance of transition 
metals.***5 However, Wilson*® has carried out a 
more detailed investigation in which he shows 
that transitions from one band to another can 
take place only if states represented by the same 
coordinates in the reduced zone scheme have 
energy differences of the order of magnitude of 
thermal energies. This is equivalent to assuming 
that the constant energy surfaces for the different 
bands intersect or are close together at a number 
of points in the reduced zone in order for Mott’s 
model to be applicable.*7 The energy contours 
obtained in the present calculations have been 
examined to see if this condition is satisfied for 
any appreciable number of points in k space and 
it was found that the energy differences between 
corresponding points in different bands** were 
much larger than thermal energies at all points 
of the momentum cell. Hence either Mott’s 
model is not applicable or the condition discussed 
by Wilson is not valid. Seitz*’ suggests that where 
d bands are involved the Ak selection rule used by 
Wilson may not be applicable. 


Phase Changes in Metallic Iron 


The phase changes which take place in iron as 
the temperature is raised are rather unusual. At 
low temperatures the stable structure is body- 
centered cubic, between 1179° and 1773°K the 
stable structure is face-centered cubic, and be- 
tween 1773°K and the melting point the stable 
structure is again body-centered cubic. 

The thermodynamic condition which deter- 
mines the relative stability of two phases at 
constant pressure is that the Gibbs*® free energy 
U-—TS+PV shall be a minimum. For phase 
changes in the solid state the PV term can 
usually be neglected. In the following discus- 
sion the body-centered phase will be denoted 
by the subscript 1 and the face-centered phase 
by the subscript 2. We then have at any 


4 N. F. Mott and H. Jones, Properties of Metals and 
Alloys (Oxford, 1936), p. 267. 

* A. H. Wilson, Proc. Roy. Soc. 167, 580 (1938). 

37 F. Seitz, Modern Theory of Solids, page 540. 

8 This result can also be obtained by reference to Fig. 2. 
The condition could be satisfied only for energies near the 
energies. 

% J.C. Slater, Introduction to Chemical Physics (McGraw- 
Hill, New York, 1939), p. 23. 
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temperature : 


—Gx(T) =[U:(0) — U2(0)] 
T 
C,—C.)dT —T 1/T)(C,—C2)dT, (5 


where C; and C, stand for the specific heats at 
constant pressure. 

If the right-hand side of this equation is nega- 
tive at a given temperature the body-centered 
structure is stable, and if it is positive the face- 
centered structure is stable. The empirical spe- 
cific heats*® as chosen by Austin?’ are shown in 
Fig. 8. Since the face-centered structure cannot 
be retained at low temperatures by quenching, 
the specific heats for this structure are not ob- 
tained directly but by a process of extrapolation 
from a series of alloys for which the face-centered . 
structure can be retained. At low temperatures it 
seems definite that the face-centered phase has 
the lower Debye temperature and hence the 
higher specific heat. At higher temperatures the 
specific heat associated with the disappearance of 
ferromagnetism gives the body-centered phase 
the higher specific heat. At the transition temper- 
ature the free energy of the two phases must be 
equal. To see how this happens Eq. (5) can be 
written 


Gi(B) —G2(B) =[Ui(0) — U2(0)] 


+ (C,:—C.)dT 
TA 
-T, f 
0 


Te 
+ f MT 


A 


TA 
+(Ts—Ta) (1/T)(C2—C,)dT. (6) 


In this equation T, is the temperature at which 


the specific heats of the two phases become equal, 


# For more recent data and for a review of the literature 
mag to Austin’s 4 4 see J. W. Awberry and 
Ezer Griffiths, Proc. a . A174, 1 (1940). Although 
it seems probable that these data supersede those quoted 
by Austin, the data quoted by him are sufficiently accurate 
for the present arguments. 
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and Tz is the temperature of the transition from 
the body-centered to the face-centered structure. 
The sum of the first three terms on the right-hand 
side is just the difference in free energies at the 
temperature 7,4. Since the body-centered struc- 
ture is stable at 7, the sum of these three terms 
must be negative. The fourth term on the right is 
also negative, so that it is the fifth term which is 
responsible for making the free energy become 
positive above T3,. It is interesting to note that 
this term is merely the entropy at 7,4 multiplied 
by the temperature interval between 74 and T3. 

If T¢ denotes the transition point at which the 
body-centered structure again becomes stable, 
we have: 


Gi(C) —G2(C) = 


TB TB 
a/ncar- f 


(7) 


In the bracket the first term represents the 
entropy of the face-centered structure and the 
second integral represents the entropy of the 
body-centered structure at Ts. Since the body- 
centered structure is stable at 0°K, and the 
entropy term increases more rapidly with tem- 
perature than the energy, the only possibility of 
making the face-centered structure stable at Tz 
is for it to have a greater entropy. Hence the first 
term on the right-hand side of Eq. (7) must be 
positive. For the body-centered structure to be 
stable above T¢ the second term must be nega- 
tive and this‘can be true only if the specific heat 
at constant pressure is higher for the body- 
centered phase. Part of this difference in the 
specific heat might come from a difference in the 
dilatation correction, but it is probable that 
most of it comes from the higher electronic 
specific heat“ in agreement with the calculations 
of this paper and the following one.'® 

A few remarks on the effect of alloying elements 
on the transition temperatures can be made. For 
some substances such as nickel and cobalt the 


'The importance of the electronic specific heat in 
making the body-centered structure stable at high temper- 
atures has been pointed out by Seitz in Modern Theory of 


, page 488 
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face-centered structure is more stable than the 
body-centered structure at absolute zero. Hence 
if we assume that the energy difference between 
the phases varies monotomically with composi- 
tion it follows that adding these elements to iron 
should widen the temperature range in which the 
face-centered or y-phase is stable. In general the 
energy difference term will favor the face. 
centered structure whenever the added element 
has more electrons per atom because the spread- 
ing of the energy bands favors the face-centered 
structure. The Debye temperature will also be 
affected by alloying elements but general remarks 
are difficult to make. The variation of the Curie 
point with composition will also have a marked 
effect on the lower transition temperature be- 
cause without the excess specific heat of de- 
magnetization the face-centered structure would 
become stable at a much lower temperature. It 
is also true that because of the greater density of 
states at the top of the Fermi distribution for 
body-centered iron it becomes ferromagnetic and 
at absolute zero there is a magnetic term in the 
free energy which favors the body-centered struc- 
ture. It thus follows that an element such as 
manganese which produces a marked lowering of 
the Curie point will, because of this effect, in- 
crease the range of stability of the y-phase. 
When there is such a close balance between 
opposing factors as has been indicated, the 
electronic contribution to the specific heat can 
have an appreciable effect on the transition tem- 
perature. This point has been discussed by 
Smoluchowski® who pointed out that if the 
addition of an alloying element shifted ¢,, to a 
portion of the m(e) curve where the density of 
states was decreased for the body-centered struc- 
ture and increased for the face-centered structure, 
this would increase the stability range of the face- 
centered or gamma-structure. The m(e) curves 
obtained in this paper and the following paper do 
not yield conclusive evidence in either agreement 
or disagreement with this theory. This may indi- 
cate either that the method used does not give 
sufficient detail to answer questions involving the 


# R. Smoluchowski, Princeton Meeting of the American 
Physical Society, December, 1941. Abstract: Phys. Rev. 
61, 390 (1942). The author wishes to thank Dr. Smolu- 
chowski for sending him a brief abstract of his talk and 
for other interesting correspondence on the subject. 
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TABLE IV. Values of 2rZ, for metallic iron. 


NNN NNN 


change of m(e) with ¢ or it may be that some of the 
other effects mentioned above may be important. 
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APPENDIX I 
Determination of Effective Field in the Metal 


Since the method of finding the effective field to use in 
computing metallic wave functions is slightly different from 
that for the free atom, it will be described briefly. Atomic 
calculations are normally carried out in terms of quantities 
which Hartree denotes by Z and Z, (Z/r? determines the 
effective force on an electron and Z,/r determines the po- 
tential energy of an electron at a distance r from the 
nucleus). Z is determined by numerical evaluation of a 
definite integral and Z, by solution of the differential 
equation: 

dZ,/dp= Z,—Z (8) 


which is integrated inwards subject to the boundary con- 
dition that Z, approaches zero at large values of r. The 
effective field of a valence electron in the metal was as- 
sumed to be due to the nucleus, the /s to 3p core, and to one 
$ type and six d type functions. The field due to the inner 
electrons and the nucleus was taken to be the same as found 
by Manning and Goldberg® for the iron atom. As a first 
approximation to the charge distribution of the other 
valence electrons the atomic charge distributions were used 
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except that they were normalized to the atomic volume 
characteristic of the metal. In computing the field due to 
such a charge distribution the ordinary method of inward 
integration of Eq. (8) is not satisfactory because of the 
difficulty in determining initial values of 2Z,. If the 
electronic charge due to the valence electrons within a 
sphere of radius r is designated by ¢, then 2¢/r? is the force 
on one valence electron due to the remaining valence 
electrons, and the inward integration of this quantity with 
respect to r yields the potential energy due to the other 
valence electrons. This method is not applicable to inner 
shell electrons because of the rapid variation of 2¢/r*. 

By use of the potential field thus found, integrations for d 
functions for a range of energies were carried out by the 
same methods as used in self-consistent field calcula- 
tions.*~** The energy for which the radial factor of the d 
function became zero and also the energy at which the 
radial derivative became zero at the half distance between 
nearest neighbors were thus determined. Since the method 
by which the first approximation was found concentrated 
the charge too close to the nucleus, the charge densities for 
these new wave functions would be concentrated too far 
from the nucleus. An attempt was made to allow for this in 
a manner similar to that used in atomic calculation by 
choosing for the next approximation a field obtained by 
averaging the initial field with one found by using wave 
functions based upon that field. A weighting of six-tenths 
initial to four-tenths final has been used. For the final field 
six different charge distributions corresponding to energies 
between that for which d,’=0 and that for which d,=0 
were used. Since previous experience with transition metals 
indicated a considerable concentration of levels near d’ =0, 
additional weight was given to levels in this region. Two 
more approximations were made and at the end of the last 
approximation the total width of the d band checked the 
width at the end of the previous approximation within a 
tenth of a unit. It was therefore decided that approxima- 
tions made in using this field would be less than other 
errors inherent in the general method. 

Since the effective potential for metallic iron may be 
useful for other calculations, the values are given in 
Table IV. 


APPENDIX II 
Computational Procedure 


After the wave functions and their radial derivatives at 
the half-distances between neighbors and next-nearest 
neighbors have been obtained, the next step is the solution 
of the transcendental equations which relate the energy and 
the components of the wave number vector k. The pro- 
cedure was similar to that used in the calculations for 


“@D. R. Hartree, Proc. Roy. Soc. Al41, 282 (1933); 
and W. Hartree, Proc. Roy. Soc. A149, 
“A. A. Bennett, W. E. Milne, and H. Bateman, Nu- 


merical of Differential ions, Bulletin of 
the National Research Council, No. 92 (1933). 
(193 3 F. Manning and J. Millman, Phys. Rev. 53, 673 
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tungsten, but one change was made which seems worth 
recording. A typical example of the transcendental equa- 
tions to be solved has the form: 


A, +A5=0, (9) 


where 7,=tan (ka/4), T:=tan (ka/2) and the A’s are 
tabulated functions of the wave functions and their 
derivatives and hence of the energy. One method of solution 
is to evaluate the expression for two values of ka/4, 
interpolate or extrapolate for a value of ka/4 which makes 
the expression vanish; try another value of ka/4 and so 
forth until the desired accuracy is obtained. Another 
method is to express 72 in terms of 7; and thus obtain a 
polynomial in 7;? which can be evaluated very rapidly by 
the method of synthetic division.“* This procedure was used 
for all of the transcendental equations except that for F; 
which has twelve different tangent factors. However, a 
table of these tangent factors as a function of the argument 
made by Dr. Chodorow in connection with the calculations 
for tungsten shortened the computations for this direction 
considerably. Details of the computational method for 
obtaining m(e) and N(e) are given in the appendix of the 
paper on tungsten.* 


APPENDIX III 


In reporting calculations such as these, some estimate of 
the approximations made is in order. The most obvious and 
probably the most important approximation is the neglect 
of all exchange effects. This precludes any detailed dis- 
cussion of ferromagnetism. Correlation effects are also 
neglected. The neglect of exchange and correlation effects 
means that calculation of the binding energy cannot be 
made and remarks about the relative stability of different 
structures can be only qualitative. The net effect of ex- 
change and correlation would be to make the energy bands 
narrower and hence n(em) larger. It is to be expected that 
the effects of exchange and correlation will be greater for 
iron than for tungsten because of the smaller number of d 
electrons in tungsten and also because the incomplete 
shielding of the f electrons in tungsten makes the electro- 
static contribution to the potential energy greater. 


“I. S. and E. S. Sokolnikoff, Higher Mathematics for 
Engineers and Physicists (McGraw-Hill, New York, 1934), 
Chapter IT. 


Another approximation is involved in the Slater method 
of fitting boundary conditions. Only fourteen functions are 
used and the continuity condition on functions and deriva- 


tives is satisfied only at fourteen points rather than over the - 


entire surface of the Wigner-Seitz polyhedron. An attempt 
to estimate the magnitude of this error has been made by 
Shockley,‘ but it is not certain to what extent his results 
are applicable to the present problem. More elaborate 
methods of calculating energy levels in metals have been 
developed by Slater,!? Chodorow,” and by Herring.** ** The 
present calculations yield a good general picture of the 
energy levels, and can serve as a starting point for more 
elaborate calculations. 

Even within the method itself certain approximations are 
made. Without much more elaborate calculations it is not 
possible to determine a completely satisfactory charge 
distribution because the criterion for self-consistency js 
more difficult to apply in the case of a solid. Another 
approximation is made in using only known lines and 
interpolating contours. Similarly, when using sections 0.1 
unit apart in finding N(e) by a volume integration in k 
space, some detail in determining the N(e) curves is lost. 
The method of determining n(e) always leads to a further 
smoothing effect so that the net effect of the last three 
approximations mentioned is to make the n(e) less irregular 
than might be the case if a more elaborate method of 
calculation were used. 

In summary, it seems probable that the general features 
of the electronic structure of metallic iron as obtained by 
this method are essentially correct and that detailed 
numerical agreement is not possible without much more 
elaborate methods of calculation. It seems probable that 
the classification of the energy bands into two responsible 
for the binding, three responsible for the magnetic behavior, 
and one for conduction is fundamental to the problem. It 
is interesting to note that the markedly greater value 
(17 vs. 11) of n(em) for the body-centered structure than for 
the face-centered structure is necessary in order to account 
for the observed differences in magnetic behavior and in 
electronic specific heat, and also to account for the high 
temperature phase change. 


47 W. Shockley, Phys. Rev. 52, 866 (1937). 
48 C, Herring, Phys. Rev. 57, 1169 (1940). 
4° C. Herring and A. G. Hill, Phys. Rev. 58, 132 (1940). 


7¢ 
78 


fo 


| 
i 
f 
¢ 
t 
Cc 
a 
c 
V 
St 
| a 


40). 


VOLUME 63, 


PHYSICAL REVIEW 


NUMBERS 5 AND 6 


Electronic Energy Bands in Face-Centered Iron 
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Jack B. GREENE* AND MILLARD F. MANNING** 
University of Pittsburgh, Pittsburgh, Pennsylvania 


Calculations of the electronic energy bands for face- 
centered iron by the Wigner-Seitz-Slater method are re- 
ported. There are found to be two filled and four partially- 
filled bands. The density of states is calculated as a function 
of the energy; the curves of the density of states vs. energy 
appear quite similar to those calculated by Slater for copper 
but resemble more closely those for body-centered iron. 
The density of states at the highest occupied energy level 
at absolute zero is 11.4 states per Rydberg unit of energy 
per atom. These results have been used to calculate the 


(Received December 19, 1942) 


average Fermi energy as a function of the total number of 
valence electrons; to investigate the ferromagnetism of Ni 
and Co by calculating the change in Fermi energy as a 
function of the number of uncompensated spins; and to 
calculate the electronic specific heat at an elevated temper- 
ature. The calculated electronic specific heat is found to 
agree fairly well with the difference from 3R of the experi- 
mentally-measured specific heat. There are no experimental 
data for the specific heat of face-centered iron at low 
temperatures. 


INTRODUCTION AND PROCEDURE 


T ordinary temperatures metallic iron has a 
body-centered structure (a-iron), but from 
870°C to 1390°C the face-centered structure 
(y-iron) is the stable one. Above 1390°C the 
body-centered structure is again stable. y-iron 
is of scientific interest principally because many 
of its properties are quite different from those of 
a-iron. The most important of these differences 
is that y-iron is not ferromagnetic. Although 
y-iron is stable only above the Curie tempera- 
ture, all extrapolations to room temperature 
indicate that it would not be ferromagnetic. 
y-iron does exist at ordinary temperatures in the 
form of austenitic steel, since the steel has been 
quenched from a temperature where y-iron is 
the stable form. 

The Wigner-Seitz-Slater cellular method of 
calculating energy bands in metals'? has been 
applied to sodium,’ to lithium,’ to copper,‘ to 
calcium,® and to one transition metal, tungsten.® 
With modifications it has also been applied to 
some insulators.’ [n this paper it has been 
applied to face-centered iron. In applying this 
method one starts with the radial part of the 


* Now at the University of Illinois, Urbana, Illinois. 

** Deceased June 1, 1942. 

1E. P. Wigner and F. Seitz, Phys. Rev. 43, 804 (1933). 

2]. C. Slater, Phys. Rev. 45, 794 (1934). 

*F. Seitz, Phys. Rev. 47, 400 (1935). 

‘H. M. Krutter, Phys. Rev. 48, 664 (1935). 

5M. F. Manning and H. M. Krutter, Phys. Rev. 51, 
761 (1937). 

®°M. F. Manning and M. I. Chodorow, Phys. Rev. 56, 
787 (1939). 

7M. F. Manning and M. I. Chodorow, Phys. Rev. 56, 
footnote 2, p. 787 (1939). 
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wave function and the radial derivative of the 
wave function at a distance half-way between 
nearest atoms in the face-centered lattice. For 
face-centered iron the lattice parameter at 900°C 
is 3.63A and the corresponding distance r 
between nearest neighbors is 4.85 Bohr units. 
The corresponding value of (p=log, 10007) is 
7.79, where r is in Bohr units. These radial 
wave functions were found by a modified Hartree 
procedure in which the Schroedinger equation 
was integrated numerically, with the same 
effective field as that used by Manning® in the 
work on body-centered iron and explained fully 
by him in Appendix I of his report. These 
integrations were carried out for s, p, d, and f 
functions over a range of energy values wide 
enough to give six complete energy bands. 


Fic. 1. The polyhedral cell for the face-centered lattice. 
8’ M. F. Manning, Phys. Rev. 63, 190 (1943). 
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Fic. 2. Energy vs. momentum for the different directions of propagation. 


The next step was to express the wave function 
within the polyhedron surrounding an atom in 
the face-centered lattice (see Fig. 1) as the sum 
of twelve terms (since there are twelve sides to 
this cell), each term being the product of a 
radial wave function times a spherical harmonic. 
The arbitrary coefficients of these terms must 
be chosen so that the wave functions satisfy the 
Bloch® periodicity condition and join smoothly 
in passing from one atom to its twelve nearest 
neighbors. The application of these conditions 
leads to a determinant, the solution of which 
gives the momentum of the electron in certain 


°F. Bloch, Zeits. f. Physik 52, 555 (1928). 


directions in terms of the radial wave functions 
and radial derivatives; this momentum of the 
electron is expressed by the wave number vector 
K. Since the radial wave functions are in turn 
functions of the energy, a relationship is estab- 
lished between the energy and the wave number 
vector or momentum in that particular direction. 
The solution of the determinant in the general 
direction is too laborious to attempt, but for a 


TABLE I. Maximum and minimum energies for each band. 


Band I II III IV Vv VI 


Max. e 0.688 0.443 0.443 0.443 0.443 0.443 
Min. 0.327 0.240 —0.0475 —0.124 —0.124 —0.525 


number of symmetrical directions the solution 
simplifies a great deal. For the face-centered 
lattice the results of the expansion of the 
determinant for a number of special directions 
have been given by Shockley.'° 


RESULTS 


The results of the calculation of the energy as 
a function of the momentum are given in Fig. 2. 
In Table I the maximum and minimum energies 
for each band are given, and in Table II the 
energies at a number of special points in the 
first Brillouin zone (the momentum cell) are 
tabulated. Figure 3 shows the first Brillouin 
zone for the face-centered lattice with the loca- 
tion of these special points indicated. (The point 
O is the center of the cell but is not shown in 
this figure.) In these tables and elsewhere in this 
paper ¢ is the negative of the energy expressed 
in Rydbergs. The curves for the 011 direction, 
for the 100 direction, and for the 111 direction 
are quite similar to those for copper.‘ These 
were the only directions for copper available for 
comparison. 

Using the energy-momentum relationships for 
the known lines in K space, one can draw fairly 
good energy contours for the principal planes in 
K space. One can then interpolate from these to 
get energy contours in the planes K,=0.1, 0.2, 
--+1.0. In each of these planes the area enclosed 


TABLE II. Values of ¢ for some points in the Brillouin zone. 


II Ill IV Vv VI 
Point 
O 0.69 0.44 0.44 0.44 0.44 0.44 
A 0.517 0.44 0.44 —0.0475 —0.125 —0.125 
B 048 044 —0.0475 -—0.125 —0.125 —0.125 
Cc 0.46 0.31 0.18 0.01 —0.05 —0.45 
D 0.325 0.29 0.29 0.02 —0.125 —0.55 
E 046 0.31 0.18 0.01 —-0.05 —0.45 


by each constant-energy contour, within which 
the energy was less than (i.e., € greater than) the 
value of the energy for this contour, was found 
by means of a planimeter. This was done for 
each of the ten planes and then the results 
integrated to get the volume enclosed by the 
surface of constant energy. This volume is then 
proportional to the number of electrons per 
atom or states per atom with energy less than 


” W. Shockley, Phys. Rev. 51, 129 (1937). 


ENERGY BANDS IN FACE-CENTERED IRON 


Fic. 3. The first Brillouin zone,for the face-centered 
lattice, showing the points referred to in Table II. 


the energy considered. The number of electrons 
with energy less than a given energy can be 
found by the normalizing condition that the 
total volume of the momentum cell represents 
two electrons per atom in each energy band. 
This quantity is called N(e). The values of N(e) 
for each of the energy bands and for the sum of 
the six bands are given in Tables IIIa and IIIb. 
It will be noticed that in the first three bands, 
values of N(e) are tabulated every 0.025 unit 
of «. For Bands IV and V they are tabulated 
every 0.05 unit of «, and for Band VI every 
0.10 unit of «. This was done because the last 
three bands extended through a wider energy 
range and would have required more contours 
than it was felt the accuracy of the method 
justified. 

From the values of N(e) the values of 
n(e)[ = —(dN/de)] can be found by numerical 
differentiation. [n(e)de is the number of energy 
levels per atom with energy between ¢ and 
(e+de).] The resulting values of n(e) for the 
separate bands and the total value of n(e) for all 
six bands as functions of ¢ are given in Fig. 4. 

Bands I and II, being completely filled, would 
not be expected to contribute either to the 
electrical conductivity or to ferromagnetism. 
Bands III, IV, and V are the ones to be investi- 
gated to see why face-centered iron is not 
ferromagnetic. Band VI looks more like a band 
predicted from the ‘‘free electron” picture than 
any of the others. It probably accounts for the 
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TABLE IIIa. Number of states having energies less than a given energy. 


Band III 


Band IV 


€ € € Nie € € N(® 
0.675 0.001 0.425 0.496 0.425 0.165 0.40 0.045 0.40 0.016 0.40 0.006 
0.65 0.009 0.40 0.905 0.40 0.349 0.35 0.146 0.35 0.056 0.35 0.024 
0.625 0.029 0.375 1.24 0.375 0.502 0.30 0.274 0.30 0.104 0.30 0.049 
0.60 0.057 0.35 1.52 0.35 0.615 0.25 . 0.430 0.25 0.182 0.25 0.082 
0.575 0.084 0.325 1.72 0.325 0.731 0.20 80.625 0.20 0.284 0.20 0.126 
0.55 0.126 0.30 1.86 0.30 0.850 0.15 1.03 0.15 0.387 0.15 0.174 
0.525 0.206 0.275 1.99 0.275 0.982 0.10 1.40 0.10 0.520 0.10 0.238 
0.50 0.399 0.25 2.00 0.25 1.21 0.05 1.62 0.05 0.684 0.05 0.304 
0.475 0.850 0.225 1.64 0.00 1.88 0.00 0.859 0.00 0.394 
0.45 1.57 0.20 1.74 —0.05 1.96 —0.05 1.13 —0.10 0.642 
0.425 1.77 0.175 1.82 —0.10 1.99 —0.10 1.60 —0.20 1.05 
0.40 1.89 0.15 1.87 —0.30 1.47 
0.375 1.95 0.125 1.91 —0.40 1.82 
0.35 1.99 0.10 1.94 —0.50 1,99 
0.075 1.96 
0.05 1.98 
0.025 1.99 
0.00 2.00 


electrical conductivity. Band VI also is the 
widest band of those investigated. In general 
features the energy bands for face-centered iron 
resemble those of body-centered iron* more than 
any other metal that has been investigated. In 
general features the total n(e€) curve resembles 
that for copper," but since the energy units are 
different for the two cases, closer comparison is 
difficult. Band VI in copper is very much like 
Band VI for face-centered iron; that is the only 
single energy band shown for copper." 

Since there are eight valence electrons per 
atom in face-centered iron, reference to Table 
IIIb shows that the highest occupied level at 
absolute zero in face-centered iron occurs at 
e=0.117. The corresponding value of m(e) is 11.4. 
Thus the density of states at the edge of the 
Fermi distribution is 11.4 states per atom per 
Rydberg unit of energy. 


APPLICATIONS 
Energy Relations 


The average Fermi energy has been computed 
as a function of the total number of valence 


. electrons per atom. The average Fermi energy 


@ is defined by the equation 
‘n 
f n(e(e—ea)de, (1) 


where ¢€o refers to the lowest level occupied at 
absolute zero and €» to the highest. Equation (1) 


uJ. C. Slater, Phys. Rev. 49, 537 (1936). 


TABLE IIIb. Summation of N(e) values for all bands. 


€ € 
0.675 0.001 0.30 5.14 
0.65 0.009 0.25 5.90 
0.625 0.029 0.20 6.78 
0.60 0.057 0.15 7.46 
0.575 0.084 0.10 8.10 
0.55 0.126 0.05 8 59 
0.525 0.206 0.00 9.13 
0.50 0.399 —0.10 10.2 
0.475 0.850 —0.20 11.0 
0.45 1.57 —0.30 11.5 
0.425 2.43 — 0.40 11.8 
0.40 3.21 —0.50 12.0 
0.35 4.35 


can be integrated by parts, giving 


1 
f N(e)de. 
N (€m) £0 
This avoids the inaccuracies that arose when 
N(e) was differentiated numerically to give n(e). 
The results of these computations for the face- 
centered lattice are shown in Fig. 5; the average 
Fermi energy is plotted against the total number 
of valence electrons per atom. In the preceding 
paper, Manning has discussed the average 
Fermi energy for the two kinds of iron and the 
possible extrapolations to predict the structure 
of the other neighboring elements in the periodic 

table. 

For the face-centered lattice the lowest band 
begins at the energy for which s’=0 and the 
other five bands begin at the energy for which 
d'=0. The behavior of the energies for which 
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NUMBER OF STATES PER UNIT ENERGY 


“70 


ENERGY IN RYDBERG UNITS 


Fic. 4. Number of states per unit energy range vs. energy. The upper curve is 
the sum of the six lower curves. 


s'=0 and for which d’=0 as a function of 
interatomic distance is shown in Fig. 6. It should 
be noted that the abscissa is r, the actual 
interatomic distance, and not p, the log, of the 
interatomic distance. 


Magnetism 


It is interesting to see whether these calcula- 
tions can offer at least a qualitative explanation 
of why body-centered iron is ferromagnetic and 
face-centered iron is not. According to the 


discussion of ferromagnetism in the paper on ~ 


body-centered iron,’ when a substance changes 
from a non-magnetic state to a magnetic state 
there is an increase in the Fermi energy, and an 
increase in the binding energy due to the increase 
in the number of exchange integrals between 
electrons of parallel spin. A ferromagnetic sub- 
stance is one in which the change in Fermi energy 


is more than compensated for by the increase in 
the number of exchange integrals. Slater" showed 
that this change in exchange energy is propor- 
tional to the square of the number of uncompen- 
sated spins. The change in Fermi energy is also 
proportional to the square of the number of 
uncompensated spins if m(e) can be regarded as 
constant, the expression for the change in Fermi 
energy being 


Aer = p?/2n(em), (2) 


where yu is the number of uncompensated spins 
in Bohr magnetons, and n(e,,) is the number of 
energy levels of both spins per Rydberg unit of 
energy at the highest occupied level for the 
non-magnetic state, and Aer is the change in 
Fermi energy upon reversing 4/2 spins. The 
smaller Aer is, the more likely it is to be compen- 
sated for by the increase in exchange integrals; 
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Fic. 5. Average Fermi energy vs. total number of valence 
electrons per atom. 


hence, the larger m(e,) is, the more likely the 
metal is to be ferromagnetic. For body-centered 
iron n(e,) is 17 and for face-centered iron (én) 
is 11.4. Hence the body-centered metal would be 
more likely to be ferromagnetic. 

If the number of states per unit energy range 
is not constant, the change in Fermi energy is no 
longer given by Eq. (2), but must be found by 
using Eq. (1) for each direction of spin. In 
computing this change Bands I and II, being 
closed bands, are neglected, and Band VI, being 
concentrated relatively far from the nucleus, is 
also neglected. In the hope of applying these 
calculations to the neighboring ferromagnetic 
metals, cobalt and nickel, the change in Fermi 
energy as u increases was computed for the total 
number of valence electrons equal to 8.0, 8.5, 
9.0, 9.5, and 10.0, respectively. Equation (1) 
was used for each direction of spin, and it was 
applied only to Bands III, IV, and V. In general, 
the Fermi energy increased less rapidly than did 
uw? but there were some fluctuations which 
indicated the computational error was rather 
large. 

The saturation magnetic moment of a ferro- 
magnetic metal is determined by one of the 
following conditions :* 


(1) The number of negative spins per atom 
cannot be negative. 

(2) The number of positive spins per atom 
cannot be greater than three. 

(3) If the exchange contribution to the energy 
increases less rapidly with yu? than does the 
increase in the Fermi energy, the difference 
between the two may vanish at some value 
of u less than that determined by either of 
the other two conditions. 


The first condition never seems to occur in any 
atom or alloy involving iron group elements, 
Since, for face-centered iron, the Fermi energy 
did not increase as rapidly as y? for the total 
number of valence electrons equal to 8.0, 8.5, 
9.0, 9.5, and 10.0, it would indicate that condition 
(3) did not determine the maximum value of y, 
Then the number of positive spins per atom 
must be equal to three; i.e., all the levels with 
plus spins are filled. Table IV shows a comparison 
of the value of the saturation magnetic moment 
calculated on this assumption and the experi- 
mentally-observed value” for Ni and Co. 

The calculated magnetic moment for Ni is 
arrived at in the following way: There are 10 
outside electrons, of which 4 are in Bands I and 
II; from the results for face-centered iron there 
are 0.64 electron in Band VI. Thus there are 
5.36 electrons in the bands that contribute to 
magnetism. If there are 3 electrons with plus 
spin there are 2.36 electrons with minus spin 
and the magnetic moment is 0.64 spin per atom, 
This agrees fairly well with the observed value 


r IN BOHR UNITS 


Fic. 6. The energy for which d’=0 and the energy for 
which s’=0 as a function of interatomic distance. 


E. C. Stoner, Magnetism and Matter (Methuen and 
Company, Ltd., London, 1934), p. 366. 
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of the saturation magnetic moment, but the 
value calculated in the same way for cobalt is 
somewhat different from the observed value. 
This might be due to the fact that the results 
for face-centered iron cannot be applied to 
hexagonal cobalt; it might be due to errors in 
these calculations; or possibly, the levels with 
plus spin are not filled. If one uses the observed 
saturation magnetic moment and assumes the 
levels of plus spin are all filled, then one finds 
that the corresponding number of electrons in 
Band VI has to be 0.71. This value is much 
larger than the value indicated by the calcula- 
tions for face-centered iron. Manning*® found 
this same difficulty with cobalt-iron alloys. He 
suggested the possibility that the electrons in 
Band VI, by coupling with electrons in other 
bands, contribute to the magnetic moment. 


Electronic Specific Heat 


The electronic contribution to the specific hea 
may be computed from the formula™ . 


(Cv)e= (5.670/27.08) X10-‘n(em)RT, (3) 


where (C,), is the electronic contribution to the 
specific heat per gram atom, R is the universal 
gas constant, JT is the absolute temperature, 
and n(n) is the number of states of both spins 
per Rydberg unit of energy at the edge of the 
Fermi distribution. 

It is only at very high and very low tempera- 
tures that the electronic contribution to the 
specific heat is appreciable. The contribution at 
very low temperatures was not computed because 
there are no experimental data available for 
comparison. For high temperatures it was de- 
cided to compare the theoretical contribution to 
the specific heat with that computed from 
experimental data at a particular temperature. 
The temperature of 960°C is well above the 
transition temperature from a- to y-iron and 
there are very good data for the specific heat of 
iron at this point. 

Specific heat measurements are made at con- 
stant pressure, but for comparison purposes one 
needs to have the specific heat at constant 
volume, C,. This is generally computed from 

%E. C. Stoner, Proc. Roy. Soc. 154, 656 (1936); also 


H. Bethe, Handbuch der Physik (J. Springer, Berlin, 1933), 
Vol. 24, Part 2, p. 430. 


the equation 
C,—C,=9a? VT/B, (4) 


where 3a is the thermal coefficient of volume 
expansion per degree centigrade, V is the 
specific volume, T is the absolute temperature, 
and @ is the compressibility. All the quantities 
except 8 are available at the required tempera- 
ture. 8 is ordinarily measured at room tempera- 
ture and extrapolated to higher temperatures, 


TABLE IV. Saturation magnetic moments 
for nickel and cobalt. 


Calculated magnetic Observed magnetic 


Element moment moment 
Co 1.39 1.71 
Ni 0.64 0.61 


but face-centered iron is unstable at room 
temperature. However, if we make use of 
Griineisen’s law, 


vy=3aV/C.B, (5) 


where y has an approximately constant value 
for all temperatures,“ then, writing Eq. (4) in 
terms of 7, we get 


C,/C.=1+3yaT. (6) 


The value of y used'® is an extrapolated value, 
since there is no experimental one available for 
face-centered iron. The value of y does not 
vary much from metal to metal; for a-iron it is 
1.68; for copper, 1.96; for nickel, 2.2; and for 
cobalt, 2.1. Since copper and nickel both have 
face-centered structures, the value for face- 
centered iron should be about 2.0. 
Substituting the following values"* of C,, a, v, 
and T in Eq. (6): 
(Cy) = 7.94 cal./g atom °C, 
(a) 960% = 25 X10-* per °C, 
7 =2.0, 
T =1233°A, 
“See F. Seitz, Modern Theory of Solids (McGraw-Hill 
Book Company, Inc., New York, 1940), p. 138. 
wa of y are from J. C. Slater, Phys. Rev. 57, 744 
16 Value of (Cp)sso’c is an experimentally-determined 


value, reported by E. Lapp, Ann. de physique 6, 826-855 
(1936); Value of (a)sso°c m from Metals Handbook (Amer- 


ican Society for}Metals, Cleveland. 1939), p. 434. 
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one gets!” 
C,=6.70 cal./g atom °C. 


If one assumes that the specific heat due to 
the lattice vibrations at this temperature is the 
Dulong and Petit value of 3R, then the quantity 


C,—3R=0.73 cal./g atom °C 


represents the electronic specific heat at this 
temperature. The theoretical value of the 
electronic specific heat computed by Eq. (3) is 
0.59 cal./g atom °C, and the agreement between 
these two values is very good. 

It is to be expected that the theoretical value 
of the electronic specific heat would be smaller 
than the other value. If the lattice vibrations 
were not entirely harmonic the value of (C, —3R) 
would contain contributions due to the an- 
harmonicity as well as the contribution of the 

17 By a method involving an extrapolation of 8, Mme. 


Lapp has calculated a value of C, to be 7.66 cal./g atom °C; 
see reference 16. 


electronic specific heat. Furthermore, the method 
of calculating the electronic energy bands re- 
ported here did not consider such corrections as 
exchange and correlation, which would tend to 
make the total spread in energy smaller and 
hence n(e) larger. This would make the electronic 
contribution greater than that found here. 

Qualitatively one can see that the electronic 
specific heat of body-centered iron would be 
greater than that of face-centered iron because 
n(ém) is greater for the body-centered iron, 
Manning® has already pointed out that it is 
because of this high electronic specific heat for 
the body-centered phase that this phase again 
becomes stable as the temperature increases. 

In conclusion it should be pointed out that 
the same approximations have been made in the 
calculations for face-centered iron as were made 
for body-centered iron. A detailed discussion of 
these approximations and their accuracy has 
been given by Manning® in Appendix III of the 
paper on body-centered iron. 
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LATITUDE EFFECT CONSIDERATIONS 


S a preliminary, we notice that if B=v/c, 

the minimum value for vertical entry 
through the earth’s magnetic field H of a 
particle of mass m and charge ne is given by 


(m/n)[Bm/(1—Bm®)* ]=ef(H), (1) 


where f(H) involves only the field and terrestrial 
dimensions, and ¢ is the electronic charge. 

We note, therefore, that if subscripts (1) and 
(2) refer to magnetic latitudes ¢; and ¢2 


Bar 


We note further that if the primary particle 
splits into mesotrons which are born at rest in 


the frame of reference of the primary particles, 
then the value of 8 for the primary particle is 
equal to that for the corresponding mesotrons, 
and (2) applies for the minimum velocities of the 
mesotrons at the two latitudes. 

Now the writer has shown! that if F(E)dE 
represents the energy distribution of the verti- 
cally directed mesotrons at the point of produc- 
tion, then the number of mesotrons at a distance 
x below this point, and resulting from the group 
representing the integrated value of F(£)dE, is 
Nz where 


N.= {(1+A/(y+1)]e* 


(3) 
1W. F. G.!Swann,' Phys. Rev. 60, 470 (1941). 


4 
q 
| 
{ 


NATURE OF THE PRIMARY PARTICLES 


with 
and 


A =apo/Amyc?, 
y 


Mo, To Po, & Em being, respectively, the rest mass 
of the mesotron, its mean life at rest, the ratio of 
the atmospheric density at x=0 to that at sea 
level, the ionization energy loss per unit of mass 
(per unit area) traversed, the minimum energy 
for entry, and where ) is the constant in the ex- 
pression p=po exp [Ax] for the relative atmos- 
pheric density at the depth x. 

Now if Ris is the ratio of the values of N, at 
the latitudes g: and ¢2, we see from (3) that 


Riz at sea level depends only upon (1—B>1), 
(1 wie and the altitude of production. We can 


express (1—Bns) in terms (1—Bn1) by (2) and 
then, Ri2 depends only upon 8, and the altitude 
of production. 

If now we repeat the process for the latitude 


ratio Riz at another altitude, we shall find 


another set of values of (1—Bar) with corre- 
sponding altitudes of production necessary to 
fit the facts. There will be only one value of 


(1—Bm1) which fits the values of Riz and Riz 
for the same altitude of production. On sub- 
mitting the matter to calculation, using the sea- 
level data? for latitudes and 0°, and A. H. 
Compton's data* for an altitude of 13,000 feet, 


we find (1—fn:)=0.028 for the latitude 40°. 
Using this in Eq. (1), and realizing that f(H) is 
known for the latitude in question, we find that 
m/n equals the proton mass. Thus the conditions 
are satisfied by a primary particle of proton mass 
and single electronic charge splitting into 10 
mesotrons,‘ since the mesotron mass is approxi- 
mately one tenth that of the proton. Moreover 
there is no other value of ” which would lead to 
a particle having the mass of any known particle. 


. Johnson, Rev. Mod. Phys. 10, 194 


.¥ A. H. Compton and R. N. Turner, Phys. Rev. 
(1937); T. 
*A. H. Compton, Rev. Sci. Inst. 7, 71 (1936); R. T. 


Young and J. C. Street, Phys. Rev. 51, 386 (1937). 

‘In the writer’s previous paper, Phys. Rev. 60, 470 
(1941), a value 5 instead of 10 was reached. This resulted 
from considerations founded upon an erroneous value 
taken for the vertical magnetic cut-off energy. 
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CONSEQUENCES OF INTENSITY-ZENITH 
ANGLE CONSIDERATIONS 


The writer has called attention® to the fact 
that the flatness of the intensity-zenith angle 
curves at high altitudes points to a decay of low 
energy mesotrons, with the creation of electrons 
whose directions of motion show wide angles 
with the directions of the primaries. 

In an unpublished memorandum, the writer 
has shown that if electrons are born from 
mesotrons in such manner that they travel out 
on the average equally in all directions in the 
frame of reference of the mesotrons, then viewed 
from stationary axes they perpetuate the direc- 
tion of the mesotrons to such an extent that the 
ratio of the electron intensity in a direction 
perpendicular to a mesotron stream to that in 
the direction parallel to it is, in the case of 
mesotrons of constant velocity 8c, given by 1— £6”. 
If the ratio has to be as great as 50 percent at 
high altitudes, as experimental data demand, 
we must have 1—f? not much less than 0.5. 

Now the value of 1—§? which evolves from 
the foregoing latitude considerations, and appro- 
priate to the mesotrons from protons as pri- 
maries is 0.028. Such a value for the mesotrons 
would result in inappreciable horizontal electron 
intensity. If, however, we consider singly charged 
helium atoms as primaries, the value of 1—? 
for the minimum energy of entry at latitude 40° 
is 0.45. Such a value of 1— can provide, 
therefore, for an appreciable horizontal intensity 
in the resulting electrons. The inference is, 
therefore, that there may be two types of 
primaries—protons which are responsible for the 
generation of the mesotrons observable at sea 
level and at semi-high altitudes, and in addition, 
heavier particles, possibly singly ionized helium 
atoms which, through their offspring electrons, 
result in the special features pertinent to the 
broad intensity-zenith angle curves for high 
altitudes. 

I wish to thank Mr. Paul Weisz for his 
discriminating computation of the data in this 


paper. 


5W. F. G. Swann, Rev. Mod. Phys. 11, 242 (1939); 
Phys. Rev. 59, 770 (1941). 
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The Probable Accuracy of the General 
Physical Constants 


FRANK BENFORD 


Research Laboratory, General Electric Company, 
Schenectady, New York 


December 7, 1942 


HE general constants used in chemistry and physics 
are a closely knit group. Physical equations, derived 
from theory, connect most of these constants in numerous 
ways and hence there are different ways of mathematically 
deriving and experimentally determining some of the more 
fundamental constants. These various values have been 
closely studied for accuracy and the published lists of the 
weighted averages of these constants usually have attached 
to each constant a value for the assumed probable error. 
To illustrate, the 1929 value assigned to c, the velocity of 
light, was 
c= (2.99796 +0.00004) 10" cm sec.~!. 


From these limits it is to be assumed that a determina- 
tion of ¢ made with equal equipment and competence 
would have a 50-50 chance of giving a value between 
2.99792 and 2.99800 for the significant figures. The latest 
published value! however, is (2.99776+0.00004)10" and 
here the departure from the previous value is 0.00020 or 5 
times the original probable error. Consulting a table of 
probability functions one finds that the odds against this 
deviation, considered as a random event, is 1341 to 1. But 
this new value of c is based on too much data to be con- 
sidered as a single event, and certainly the forethought and 
care that go into a determination of the speed entitle it 
to a better adjective than random; therefore, it might be 
expected to be in better agreement. It thus becomes ap- 
parent that the assigned probable limits do not apply to 
the absolute accuracy of c, but only to its accuracy as 
judged by internal evidence of the experiments. If there 
are systematic errors of conception or execution the whole 
group will be shifted away from the true value without 
increasing the apparent probable error. The given probable 
error then is no longer a real measure of accuracy, but is 
reduced to being an indicator of the consistency of the 
data that went into the computation of the assigned value. 
The distinction between consistency and accuracy is im- 
portant, and it should not be overlooked. 
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The writer has been working in radiation and using the 
radiation constants ¢; and ¢z and also Wien’s displacement 
constant A and the Stefan-Boltzmann constant g¢. It is 
these particular constants that have suffered the greatest 
proportional change in Dr. Birge’s latest summary, as is 
illustrated in Table I. Here only significant figures have 
been retained, the various appended values of 10* and the 
zeros preceding the probable errors being dropped for 
better readability. The selection of constants was intended 
to cover those of personal interest and a few of the more 
important fundamental ones. 

The main contributing factor to the changes between the 
1929 and 1941 list is the new value assigned to the elec- 
tronic charge. In 1929 it was 4.7700+50 and in 1941 it 
was 4.8025+10, a difference of 325 as compared to the 
probable error of 50. The change is 6.5 times as great as 
the probable error, and the chances against such a change 
are 100,000 to 1, on the basis of the 1929 figures. It is here 
again evident that the assigned limits refer to the con- 
sistency of the data from which 4.7700+50 was derived, 
and the present value would, in 1929, have seemed im- 
possible from the internal evidence. 

One possible criterion of the accuracy of the constants 
would be to note in successive lists the number of the new 
values that fall within the preceding limits. If half are 
within and half without then the limits are becoming true 


Taste I. Probable values of the general physical constants (from Birge). 


Change 
1929 1941 1929 P.E. 1to 
Velocity of light ¢ 2.99796 +4 299776 34 50 13-16 
Gravitational constant 6.664 2 6.670 5 30 20-10 
Volume of ideal gas Vo 22.4141 8 22.4146 6 0.6 <1 
International ohm P 1.00051 2 1.00048 2 15 20-10 
Ice point To 273.18 3 273.16 1 07 <1 
Joule eq. (mech.) Jis 4.1852 6 4.1855 4 05 <1 
Joule eq. (elec.) Jis’ 4.1835 7 4.1847 3 #17) 
Faraday constant 
(chem) F 96.494 5 96.501 10 20-10 
Faraday constant 
F 96.489 7 96.514 10 3.6 65-10 
Electronic charge 4.7700 50 4.8025 10 865 10-16 
tro- 
e/m 1.7610 10 5 18 30-10 
(deflection) 3034 60 5.2736 15 50 13-10 
nck constant 5470 80 6.6242 24 
Ratio es to em units bb 99790 10 2.99776 4 14 2.0-10 
Rydberg constant Ry 109677.758 50 109677.581 7 3.5 54-10 
Avogadro's number No 6.0644 60 6.0228 7.0 40-16 
Boltzmann constant K 1.37089 140 1.38047 26 
Stefan- Boltzmann 5.7139 60 5.6728 37) 6.9 
Wien’s displacement A 2.8836 7 


Secondradiationconst. 2 1.43170 


indicators of accuracy, and the truth is being approached 
by steps similar to the sum of the series $+}+ +>" 
which after an infinite number of terms (lists) will approach 
1, the symbol of perfection. 

The debt of chemists and physicists to Dr. Birge is 
great, and I hope that he will, in future lists, increase this 
indebtedness by pointing out the distinction between 
consistency and accuracy. 


1 Raymond T. Birge, Rev. Mod. Phys. 13, 233 (1941). 
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LETTERS TO 


Comments on “‘The Probable Accuracy of the 
General Physical Constants” 


RAYMOND T. BIRGE 
Department of Physics, U niversity of California, Berkeley, California 
February 17, 1943 


N the preceding paper by Dr. Frank Benford, attention 
is called to the important distinction between con- 
sistency and accuracy. The recent large and wholly un- 
anticipated changes in the measured values of many of 
the general physical constants have been noted and dis- 
cussed in numerous papers and it would be superfluous 
to make further comments, were it not for the fact that 
Dr. Benford’s paper concludes with the remark, “I hope 
he will, in future lists, increase this indebtedness by point- 
ing out the distinction between consistency and accuracy,” 
For many years now I have attempted to emphasize 
just this distinction. Unfortunately the true value of any 
physical magnitude can never be known. Hence the abso- 
lute accuracy of a measured result cannot be determined 
and we can only note the difference, if any, between the 
internal and external consistency of the data. In 1932 I 
initiated a discussion! of these two important concepts and 
noted certain relations and facts that had hitherto been 
generally overlooked. 

Previous to the time that I began publishing critical 
values of the general physical constants, it was rather 
exceptional to attach a probable error, or other measure of 
reliability, to suggested “‘most probable” values. It seemed 
to me, however, desirable to attempt such an estimate, 
even in cases where the estimate was admittedly almost a 
pure guess. On the other hand, as I have noted,? F. W. 
Clarke, whose calculated atomic weight values were con- 
sidered authoritative, prior to the formation of national 
and international committees on the subject, regularly 
published probable errors. Unfortunately he ignored all 
discrepancies between internal and external consistency 
and as a result adopted a system of weighting that seems 
quite unjustified and misleading. 

The table published by Dr. Benford shows that the 
average change of value, between 1929 and 1941, is some 
five times the 1929 published probable error. But in 
Clarke’s work on atomic weights, I found that the probable 
error based on external consistency averaged fully ten 
times that based on internal consistency and in at least 
one case was one hundred times as large. In other words, 
atomic weight determinations prior to 1920 were afflicted 
with unknown and unsuspected systematic errors that 
averaged fully ten times the assumed experimental errors. 

It is unfortunate that even the most carefully conducted 
precision work on physica! constants seems, with rare 
exceptions, to contain errors considerably larger than any 
suspected by the experimenters themselves. No one has 
better reason to realize and to regret this fact than the 
present writer. I cannot, however, refrain from pointing 
out that I have always tried, in my estimates of probable 
error, to make allowance for systematic errors, and have 
also more than once noted and made correction for errors 
overlooked by the original authors. As a result of the 
diversity of methods now used in determining most of the 
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general constants, and of the greater attention being paid 
to the question of systematic errors, I feel that the values 
and probable errors of the 1941 list can be accepted with 
considerably greater confidence than in the case of any 
earlier list. One may hope that never again will it be 
necessary to base an adopted value of an important con- 
stant on the results of a single method, carried out chiefly 
by a single experimenter. 


'R. T. Birge, Phys. Rev. 40, 207 (1932). See pp. 213-223. Further 
pertinent material on the general subject matter of this note appears 
: pp. 122-123 of the paper by W. E. Deming and R. T. Birge, Rev. 


od. Phys. 6, 119 (1934), and in the Introduction (pp. 1-8) of my 1929 
paper on the general constants, Rev. Mod. Phys. 1, 1 (1929). My de- 

led discussion of the 1941 list of constants appears in Reports on 
iy) in Physics 8, 90-134 (1941). 

2 Page 221 of first reference of footnote 1. 


Cosmic-Ray Stars at 10,000 Feet 


Wayne E. Hazen 
University of California, Berkeley, California 
February 23, 1943 


PRELIMINARY survey of the stars observed in a 

cylindrical cloud chamber (30 X 30 cm) at Tioga Pass, 
Yosemite National Park, has been completed. The chamber 
contained eight 0.7-cm lead plates separated by 2.5-cm air 
spaces. Out of 9100 expansions, none of which was counter 
controlled, approximately 8500 useful photographs were 
obtained. 

A total of 58 stars were observed, two of which originated 
in the gas and the others in the lead plates. The number 
of penetrating particles that occurred in the same volume 
and about the same time was 19,000. At 14,000 ft. with 
similar apparatus Powell' observed 156 stars and 9840 
penetrating particles. A direct comparison indicates a de- 


Fic. 1. Phot phs of stars that occurred slightly before the cham- 
ber expansion. The star in the upper photograph contains both long 
and short range particles while the one in the lower photograph con- 
tains only short range particles. 
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160) 
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RANGE IN LEAD cu 


Fic. 2. Range Gpciation of the particles that occurred in 58 star 
showers at 10,000 ft. 


crease in the relative abundance of stars by a factor of five 
in going from 14,000 feet down to 10,000 feet, whereas the 
soft component decreases by only a factor of two. It is en- 
tirely possible, however, that the conditions of observation 
were sufficiently different in the two cases so that the factor 
for decrease in relative abundance of stars might also be 
nearer two. 

Thirty-one of the stars contained only particles with 
range less than 0.7 ¢m of lead and consisted of 3.3 particles 
on the average. Only one or two of these stars appeared to 
be initiated by an ionizing ray and none of them was asso- 
ciated with cascade showers. The remaining twenty-seven 
stars contained particles with range greater than 0.7 cm of 
lead in addition to shorter range particles and consisted of 
5.1 particles on the average. Eight of this latter group of 
stars apparently were initiated by ionizing rays but again 
none of them was associated with cascade processes. Photo- 
graphs of typical stars are reproduced in Fig. 1. 


The only previously reported data on the range distribu- 
tion of star particles have been obtained from tracks in 
photographic emulsions.* In Fig. 2 the range distribution of 
particles from the present observations is compared with 
the results of Widhalm’s observations with photographic 
emulsions. Most, if not all, of the long range particles 
observed with the cloud-chamber technique could not be 
detected in photographic emulsions for two reasons: First, 
since the ionization of the long-range particles was not much 
greater than the minimum, they probably could not be ob- 
served in the emulsion, and secondly, most of them would 
pass out of the emulsion before coming to rest. The last 


' four blocks in the diagram of Fig. 2 should be increased by 


factors of about 1.3, 2, 3, and 4, respectively, in order to 
take into account the decreasing solid angle subtended by 
the final plate as we proceed away from the star's origin, 
When this correction is included we find that nearly half 
of all the star particles have ranges exceeding 0.7 cm of lead. 

The number of star particles that might be interpreted 
as electrons is extremely small, certainly less than ten for 
all particles and less than five for particles of range greater 
than 0.7 cm of lead, while the number of particles in the 
latter group that are most likely to be electrons is even 
smaller. The infrequency of occurrence of particles in stars 
that are most likely to have been electrons is in essential 
agreement with the few previous observations of star 
showers employing a technique capable of ascertaining the 
nature of the particles. The present evidence indicates that 
most of the stars are probably initiated by neutrons, the 
others by penetrating rays (mesotrons or protons). 

The above work was made possible through the help- 
fulness of F. A. Kittredge, Superintendent of Yosemite 


National Park. 


1W. M. Powell, Phys. Rev. 61, 670 (1942). 
2A. Widhalm, Zeits. f. Physik 115, 481 (1940). 
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Proceedings of the American Physical Society 


“ANNUAL MEETING OF 1942,”” NEw YorK, JANUARY 22-23, 1943 


HE “Annual Meeting of 1942,’’ postponed 

from December, 1942 because of the con- 
gestion of trains expected to prevail during the 
holiday season, was held at Columbia University 
in New York City on January 22-23, 1943. 
The American Association of Physics Teachers 
(AAPT) and the American Society for X-Ray 
and Electron Diffraction (ASXRED), having 
postponed their meetings as we did ours, met 
on the same days and shared in some of our 
meetings, extending to us a reciprocal hospitality. 
Notwithstanding the sudden and belated change 
to a time normally regarded as much less favor- 
able, the meeting was surprisingly and delight- 
fully successful, ranking in size of attendance 
with the best winter meetings of the past. 
Nearly three hundred of our members registered, 
and another hundred who were members of the 
AAPT or the ASXRED or simply guests. The 
joint-dinner of the three societies was attended 
by nearly two hundred. Excellent arrangements 
were made by H. W. Webb as Chairman of the 
Local Committee. 

The special features of the meetings, attended 
by all three societies, comprised the Retiring 
Presidential Address of President P. W. Bridg- 
man of the Physical Society; the presentation 
of the Oersted Medal of the AAPT to G. W. 
Stewart, the response of Mr. Stewart, and the 
Richtmyer Memorial Lecture delivered by 
G. F. Hull, Sr.; and an address by P. Debye, 
arranged by the ASXRED, on the topic “Recent 
Developments in X-Ray and Electron Diffrac- 
tion.” Forty-three ten-minute papers were 
contributed, of which the abstracts are appended 
(five of them coming by way of the ASXRED); 
those numbered 5, 18, 39, 40, 41, 42, and 43 
were read by title. 

At the Business Meeting of the Society, short 
reports were made by or in the name of the 
Secretary (who gave the present membership of 
the Society as 4112, comprising 917 Fellows and 
4 Honorary Members) ; the Treasurer (whose full 
report will soon be sent to the membership); 
and the Acting Editor (who said that it is 


planned to continue issuing one number per 
month of the Physical Review for the present). 
The Secretary also announced that the member- 
ship list to be issued in 1943 will be confined to 
the names of members elected since the last 
previous list was issued. The Tellers of the recent 
election (R. T. Cox and F. C. Nix) reported 
through Mr. Cox the election of the following 
officers: President, A. W. HuLL; Vice President, 
A. J. Dempster; Secretary, K. K. Darrow; 
Treasurer, G. B. PEGRAM; Members of the Coun- 
cil: R. B. Linpsay, ALPHEUs W. Smita; Members 
of the Board of Editors: M. Putuips, B. Rosst, 
W. H. ZACHARIASEN. These entered upon their 
terms at the end of the meeting. 

The Council met on Thursday evening, Jan- 
uary 21, and again on Saturday evening. Two 
Members were advanced to Fellowship; two 
candidates not previously Members were elected 
to Fellowship; 60 candidates were elected to 
Membership. The names of all these are ap- 
pended. The Council appointed G. B. Pegram 
to succeed himself as Member of the Governing 
Board of the American Institute of Physics. 

A Joint Resolution originating from the 
American Institute of Physics was submitted to 
the Business Meeting of the Society by the 
Council thereof, and to the corresponding meet- 
ing of the AAPT by the authorities thereof, 
and passed by both Societies. The text of this 
resolution is appended. It has since been trans- 
mitted to the officials named in its text and to 
many others. 

The Society has lost through death the follow- 
ing Fellows and Members: H. G. Beutler of the 
University of Chicago, noted in spectroscopy; 
C. A. Hall of the Electrical Storage Battery 
Company in Philadelphia; James E. Ives formerly 
of the U. S. Public Health Service; G. B. Kare- 
litz of Columbia University. 

Members advanced to Fellowship: Arthur Adel, 
Frank Morgan. 

Other candidates elected to Fellowship: Pierre 
Auger, S. Chandrasekhar. 

Candidates elected to Membership: James Henry 
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Awbery, Arnold M. Bass, J. G. Bennett, Sey- 
mour Benzer, Daniel Bershader, William H. 
Billhartz, John Markus Blatt, Lyle B. Borst, 
Robert Stephen Burpo, Jr., Owen Chamberlain, 
Rev. Ernest Clarke, Alden P. Cleaves, Howard 
S. Coleman, Norman Davids, Pauline V. David- 
son, James J. Donoghue, Bernard Epstein, 
Herman Theodore Epstein, Vincent Fitzgerald, 
Elizabeth F. Focht, Edward M. Fryer, Donald 
N. Gideon, Gerhart Groetzinger, Harold Palmer 
Hanson, Rudolph G. E. Hutter, J. Wesley Kemp, 
Capt. Paul J. Kopp, Louis W. Labaw, Robert W. 
Leonard, John Holton McGinn, Arthur E. 
Middleton, Elliott W. Montroll, Thomas M. 
Moore, Herbert C. Muether, Raymond L. 


*Murray, Edward J: Nagy, J. Burton Nichols, 


Newton Hays Odell, William Elwood Ogle, H. 
Lowell Olsen, Francis Perrin, Fritz Pordes, H. 
Eugene Powell, Paul P. Reichertz, Dale W. 
Rinehart, Salomon Rosenblum, H. Gunther 
Rudenberg, Noel W. Scott, Gerald Lionel 
Simard, William Warner Sleator, Jr., John B. 
Smyth, A. O. Stanesby, Richard G. Stephenson, 
Rodolphus A. Swan, Jr., Vernon Thornton, 
Charles B. Vance, George Warfield, Richard 
Nicholas Work, Kenneth P. Yates, Ignace 
Zlotowski. 

Joint resolution of the American Physical Society 
and the American Association of Physics Teachers: 


WaerEas the Army and Navy College Training Programs 
as so far announced provide for the training of technical 


personnel for use almost solely, if not entirely, within the 
Armed Forces, and 


WuEREAs other war agencies, whose effective functioning 
is essential to the support of the Armed Forces in expediting 
the prosecution of the war, also need technically trained 
personnel in large numbers, and 


WuerEas these other war agencies, such as war industry, 
war research, and war training staffs are now confronted 
with the closing off of the usual avenues of supply of 
technically trained personnel from the same groups now 
being drawn upon through induction and enlistment. 


THEREFORE BE IT RESOLVED that the American Physical 
Society and the American Association of Physics Teachers 
in joint session in New York City on January 23, 1943, 
urge that effective steps be taken at the earliest possible 
date to provide an over-all War Training Program to 
meet all of the needs of the Army, of the Navy, of War 
industry, and of war research in scientific and technological 
fields. 


That in such a program adequate provision be made to 
insure such a flow, of personnel trained in the field of 
physics into war research, into war industry, and into 
training staffs, as will be necessary to prevent serious 
impairment of the critical war services to be rendered 
by these agencies. That copies of these resolutions be sent 
to the Chairman of the War Manpower Commission, to 
the Secretary of War, to the Secretary of the Navy, to the 
Chairman of the War Production Board, to the Director 
of the Office of Scientific Research and Development and 
the Chairman of the National Defense Research Committee. 


Kart K. Darrow, Secretary, 
American Physical Society 
Columbia University, New York, New York 


ABSTRACTS OF CONTRIBUTED PAPERS 


1. Decomposition of Matter Through the Magnet (Mag- 
netolysis). FELIX EHRENHAFT.—I have stated already that 
microscopic bodies, especially of ferromagnetic elements, 
do move in homogeneous magnetic fields, reversing their 
direction with the reversal of the field. In the spirit of M. 
Faraday and J. C. Maxwell one must therefore conclude 
that these test bodies carry an excess of north- or south- 
magnetic charge and consequently have to be called mag- 
netic ions. Therefore we have to establish the equation; 
Div B=0. A consequence of these facts is the phenomenon 
of magnetophoresis already published.* Lately I have dis- 
covered that the magnet causes chemical composition and 
decomposition of matter in several cases. The quantity of 
gas, originated in certain chemical processes occuring be- 
tween the magnetrodes and the magnetolyte (e.g., acidu- 
lated water) is considerably increased by the application 
of the magnetic field. The gas created is a mixture of 
hydrogen and oxygen. Herewith is proved that the poles 


of a magnet break up water through magnetolysis, the 
same as electric poles break up water by electrolysis. Now 
the way is open to measure the intensity of the magnetic 
current. 

It should be mentioned that the same electric potential 
of the magnetrodes is secured, be it by using one piece of 
metal for the magnet, or two pieces which are metallically 
connected. 

* Science 94, 232 (1941): 96, 228 (1942); Phys. Rev. 61, 733 (1942). 


2. The Effect of Oblique Incidence on the Conditions 
for Single Scattering of Electrons by Thin Foils. GeraLp 
GOERTZEL AND R. T. Cox, New York University.—In 
setting up and applying criteria for the single scattering of 
electrons by thin foils, it has been usual to assume that the 
principal deviations from single scattering are caused by 
the combination of small deflections with one large one. 
It is here shown that with electrons obliquely incident on 
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the foil much more serious deviations may be caused by 
the combination of two deflections of the same order of 
magnitude. It is concluded that the difference observed by 
Chase and Cox! between the scattering at 90° to the two 
sides of a foil on which electrons are incident at 45° may 
be ascribed to this cause. The bearing of this factor on the 
negative result of Dymond’s* experiment on the polariza- 
tion of electrons is also considered. 


Phys. Rev, 58, 243 (1940). 
Proc’ Roy. Soc. 687 (1934). 


3. The Small Angle Scattering of Electrons by Alu- 
minum. BowEeN C. DEES AND BERNARD HAMERMESH, New 
York University.—The angular distribution of scattering of 
46 kv electrons scattered by an aluminum foil of 1.3-10-* 
cm thickness has been studied in the angular range between 
0° and 10°, The mean projected angle of scattering was 
determined experimentally for the cases of normal and 45° 
incidence, and the resulting experimental values were com- 
pared with the values predicted by a theoretical treatment 
given by Goudsmit and Saunderson.' Agreement with 
theory was found in each case. Attempted comparison with 
the theoretical treatment given by Williams? led to results 
which showed that his treatment was inapplicable to the 
circumstances of this experiment, and forced the conclusion 
that the scattering observed in this experiment was multi- 
ple for a very small angular range at most. Finally, in- 
vestigation showed that the scattering found did not ap- 
proach single scattering even at the largest angles of ob- 
servation, and led to the conclusion that the observed 
scattering was mainly plural. 

1 Goudsmit and Saunderson, Phys. Rev. 57, 24 (1940); 58, 36 Use. 


asus Williams, Proc. Roy. Nay 160A. 531 (1939); Phys. Rev. 58, 29. 


4. Motion of the Copper Arc in Transverse Magnetic 
Field. Caartes G. SmitH,* Raytheon Mfg. Co.—An arc 
between equal parallel rings of copper is subjected to a 
transverse radial magnetic field. With lower ring as cathode 
and south pole at center, the arc viewed from above rotates 
counter-clockwise at atmospheric pressure, but below a 
critical pressure it rotates clockwise. The pressure for 
reversal varies approximately inversely as the distance be- 
tween rings. Critical air pressures ranging from 36 to 3 cm 
were observed. The peculiar clockwise rotation results from 
the thermomotive force or Righi-Leduc factor operative in 
the solid cathode, according to one emission theory' of the 
mobile arc cathode. Mass movement of ionized gas caused 
by electromagnetic force tends to drive the arc around 
counter-clockwise. This tendency becomes less at reduced 
pressure since viscosity factors increase due to spreading 
of the arc column. Below the critical pressure, phenomena 
in the solid cathode predominate. The observations ap- 
parently extend the emission theory to include the copper 
arc over a wide range of pressure. 


ormed at Harvard University. 
hys. Rev. 62, 48 (1942). 


5. The Adsorption of Thorium on Tantalum. CHARLES 
J. GALLAGHER, University of Notre Dame:-—lIt is known 
that thorium adsorbed on the surface of tungsten or mo- 
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lybdenum results in a decrease in the work function and 
also a decrease in the value of A, the emission constant.! 
The behavior of thorium on tantalum was investigated and 
similar effects were noted. Thorium was evaporated from 
thoriated tungsten onto a well-outgassed tantalum ribbon 
and the thermionic emission at 1390K studied as a function 
of the time of evaporation. The work function of the 
tantalum changed from 4.07 volts to 2.52 volts, while A 
decreased by a factor of thirty. The ribbon was then held 
at 1550K and the emission from the two sides studied 
separately. The emission from the activated side decreased, 
while that of the clean side increased until both reached 
the same final level after about twenty hours. This indicates 
that the adsorbed thorium had migrated around the edges 
of the ribbon and distributed itself uniformly over the 
emitting surface. 


1A. L. Reimann, zpeuteels Emission (John Wiley and Sons, New 
York, 1934) pp. 99, 


6. The Action of Oxygen, Nitrogen and Hydrogen in 
Lowering the Photoelectric Work Function of Zirconium, 
Titanium and Thorium. H. C. RENTSCHLER AND D. E. 
HENRY, Westinghouse Lamp Division.—Pure wire surfaces 
of zirconium, titanium and thorium were made by thermal 
decomposition of the iodides. Photoelectric surfaces of 
the pure metals were prepared by the method previously 
described' and the work functions of these metals were 
determined from the photoelectric threshold wave-lengths. 
Zirconium and titanium when heated in pure oxygen, 
nitrogen or hydrogen dissolve considerable amounts of 
these gases. The photoelectric work functions of surfaces 
similarly prepared from the metals with the dissolved 
gases were found to be appreciably lower than for the corre- 
sponding pure metals and the photoelectric response from 
such surfaces is as constant as that from the pure metal. 
Thorium reacts with these same gases to a lesser extent 
and the shift in the photoelectric response is similar to that 
for zirconium and titanium. 


1 Rev. Sci. Inst. 3, 794 (1932). 


7. The Photoconductance of Evaporated Bismuth Films. 
A. H. WEBER AND L. M. FRIEDERICH, S. J., St. Louis Uni- 
versity.—Observations were made of the dark conductance 
and photoconductance* of films deposited on Pyrex in high 
vacuum, continuously and discontinuously, both at room 
temperature and at liquid-air temperature. A photocon- 
ductance effect, observed in all the films deposited, dis- 
appeared in every instance at a critical thickness, 11 to 292 
atom layers depending on the conditions of deposition. 
Films deposited continuously, both at room temperature 
and at liquid-air temperature, attained a greater respective 
thickness (2 to 3 times) before the photoconductance effect 
disappeared than did those deposited discontinuously. 
Films deposited at room temperature lost the photocon- 
ductance effect at a much greater film thickness (8 to 15 
times) than did those deposited at liquid-air temperature. 
A film no longer exhibiting photoconductance at liquid-air 
temperature, when allowed to warm to room temperature, 
again showed a photoconductance effect. These observa- 
tions lead to the conclusion that bismuth films exhibiting 
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the photoconductance effects have a patch structure and 
that the increased electrical conductance when illuminated 
is not true photoconductance, but photoelectric emission 
from patch to patch. 

* Defined as the electrical conductance of the films when illuminated. 


8. The Compensation of Random Fluctuations in a D.C. 
Bridge Amplifier. J.C. M. BRENTANO AND E. R. SCHLEIGER, 
Northwestern University—In measuring small potential 
differences or currents it is sometimes necessary to use d.c. 
amplification. Instances are the measurement of weak 
ionization currents or photoelectric phenomena involving 
delay actions where a.c. amplification cannot readily be 
adopted. A compensated symmetrical bridge using two 
thermionic tubes has the advantage over single tube 
circuits that it is easy to control; with it all disturbances 
and fluctuations can be eliminated very effectively, except 
fluctuations due to random changes of filament emission 
which actually set the limit to the ultimate amplification 
attainable. In a bridge* comprising two vacuum tubes, 
essentially of tetrode type, the intermediate grids acting 
as control grids, compensation for such fluctuations can 
be obtained by a system of cross connections between the 
inner grids and plates by which an emission fluctuation of 
one filament produces similar changes in the symmetrical 
halves of the bridge. The unsteadiness level can thus be 
reduced by a factor of thirty. 

* Brit. Pat. 526869. 


9. Graphical Treatment of Non-Linear Circuits. P. I. 
Wo Lp, Union College.—With the increase in the use of 
circuits containing non-linear elements, it becomes im- 
portant to make calculations as to the behavior of any 
such circuit. Most non-linear elements have a volt-ampere 
characteristic which can be given by an infinite series only. 
The analytical treatment, therefore, becomes excessively 
burdensome unless all but a very limited number of terms 
of the series are neglected. This involves strong approxi- 
mations. Graphical methods of treatment will be presented 
which avoid approximations and permit the ready plotting 
of the current characteristics of any branch in simple 
networks. 


10. The Ionization Cross Section of the Silver L111 State. 
J. J. G. McCue,* Cornell University.—The cross section of 
silver atoms for excitation to the Ly state by cathode 
rays has been measured in arbitrary units as a function of 
cathode-ray energy, by using the thin target technique 
originated by Webster. The cathode-ray energies covered 
the range 1.5=>U=9, where U=(cathode-ray energy) /(ex- 
citation energy). The cross section was determined from 
the La doublet intensity, measured as a function of tube 
voltage and corrected for the effects of diffusion, rediffu- 
sion, retardation and fluorescence. The values of cross 
section computed by Burhop! by using the Born approxi- 
mation and hydrogenic wave functions, vary in approxi- 
mate agreement with the experimental results; the agree- 
ment however is not within the estimated experimental 
uncertainty of 3 percent. Comparison of the present results 
with those of Webster and his colleagues for the K state 


establishes that for both the K and the Lyz; states of 
silver, the ionization cross section is a maximum when JJ 
is about 3.5, and that for U>3.5 the decrease in ionization 
cross section with increasing cathode-ray energy is more 
rapid for the L;; state than for the K state. 


* Introduced by L. G. Parratt 
Burhop, Proc. Phil. 36, 43 (1940). 


11. Fine Structure in the X-Ray Absorption Edges of 
Potassium Chloride. J. W. TriscHKA, Cornell University — 
Variations in the absorption coefficient in the regions of 
the K edges of the potassium and chlorine in potassium 
chloride have been measured with a two-crystal vacuum 
x-ray spectrometer. The region studied extended from ten 
volts less than the energy of the first absorption line to 
forty volts greater. Absorbers were prepared by evapora- 
tion of the potassium chloride onto thin kodapak films, 
The lower limit of detectable variations in the absorption 
coefficients of the potassium and chlorine was three percent, 
Results were not corrected for the resolving power of the 
spectrometer. Using Sb La; and Pd La; as reference lines, 
the wave-lengths of the peaks of the first absorption lines 
were measured to be 3425.5 Siegbahn X.U. and 43788 
Siegbahn X.U. for potassium and chlorine, respectively, 
The separation of the first two maxima in the chlorine 
edge was measured to be 3.9 volts as compared with the 
theoretical 4.2 volts calculated by Landshof' on the 
assumption that the maxima correspond to excitation 
states in the solid. 

' R. Landshof, Phys. Rev. 55, 631 (1939). 


12. An Error Arising from the Use of Small Carbon- 
Lined Ionization Chambers in the Measurement of Low 
Voltage X-Rays in Roentgens. H. QuimBy anp 
F. Focut, Memorial Hospital—With 44-ky 
x-rays, small ionization chambers of lucite lined with 
Aquadag give lower values for ionization per cc than open 
air chambers. In a parallel plate chamber of the lucite- 
Aquadag type, with electrodes perpendicular to the x-ray 
beam, ionization per cc increases steadily with electrode 
separation from 0.2 mm to 2.0 mm; after that, the incre- 
ment of ionization per unit increase in volume is constant. 
The effect is attributable to the difference in the photo- 
electrons emitted by the carbon-coated electrodes, and by 
the air between the electrodes. When one electrode is 
made of aluminum the effect is reversed, the ionization per 
cc decreasing from 0.2 to 2.0 mm spacing. Theoretically, 
the ratio between the readings in carbon and aluminum 
chambers of this type, with very small spacing, can be 
calculated on the basis of Kramer’s formula, the accepted 
formulae for .r and .a, and published absorption data. 
The ratio thus obtained is in satisfactory agreement with 
the experimental results. It is concluded that, for the 
measurement of x-rays in the voltage range in which photo- 
electric absorption is important, small carbon-lined ioniza- 
tion chambers do not satisfy the requirements imposed by 
the definition of the roentgen. 


13. Electromagnetic Method for the Determination of 
Velocity Distribution in Fluid Flow. ALEXANDER KOLIN, 
Columbia University—A method has been developed which 
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utilizes the principle of electromagnetic induction for 
determining the velocity distribution in a liquid. An 
homogeneous magnetic field is maintained at right angles 
to the direction of fluid flow. A potential gradient is 
established in the liquid by induction. Its magnitude at 
any given point is: dE/dz= 4H Xv. The local instantaneous 
velocity component in a direction perpendicular to the 
magnetic field and to the Z axis is determined by measur- 
ing dE at a given point. This is accomplished by introduc- 
ing a probe consisting of two fine wires, insulated except 
at the tip (dz, the separation of the wire tips, is in the 
order of a few thousandths of an inch). The induced e.m.f. 
picked up by the probe is amplified and photographically 
recorded. The probe can be placed at any point where the 
instantaneous local velocity is to be measured. Since the 
method has no lag, it furnishes the velocity distribution in 
time as well as in space. To avoid difficulties due to polariza- 
tion, an alternating magnetic field is utilized. This work 
was sponsored by the Department of Civil Engineering, 
and the author wishes to thank Professor B. A. Bakhmeteff 
for guidance and suggestions, and Professor J. K. Finch for 


his interest and generous support. 


14. Molecular Absorption of Vibration. C. W. CHAMBER- 
Lain, Michigan State College—The solution of vibration 
problems by means of springs is complicated by limitations 
introduced by resonance, load carrying capacity, stability, 
and dampening factor. Springs constructed of organic 
materials like rubber in shear, sponge rubber, cork, and 
highly compressed felt obey, to a greater or less degree, 
the same law as springs made of metal. A gas-filled cylinder 
with movable piston obeys the gas law and possesses high 
load carrying capacity. Its dampening factor, due to the 
fact that cooling of the gas on expansion is less than the 
heating on compression, is small. It has been found feasible 
to construct a series of gas-filled cylinders so short that 
the low heat conductivity of the confined gas introduces no 
limitation. Work done on the confined gas by the moving 
piston is directly transmitted to the cylinder base and 
rapidly disposed of by conduction. In operation although 
the initial temperature of the gas is high the temperature 
of the assembly of cylinders is comparatively low. Highly 
efficient absorption cushions composed of 100 cylinders in 
series and only } inch thick have been successfully in- 
stalled under machines varying in weight from a few 
pounds to 500 tons. 


15. On the Excitation Functions of Proton-Induced Re- 
actions in Lithium. E. D. Courant, University of Rochester. 
—Evidence for the hypothesis that the ground state of 
sLi’? is odd has been obtained by comparing the experi- 
mental excitation functions of the reactions 3Li’ (p, y)«<Be* ! 
and 3Li’(p, 2a)? to the computed penetration probabilities 
for protons of angular momentum 0 and 1. If sLi’ is odd, 
the y-reaction is produced by protons of even angular 
momentum and the a-reaction by protons of odd angular 
momentum, and vice versa. When the effect of the Cou- 
lomb potential barrier is eliminated by dividing the ob- 
served excitation functions by the penetration factor, the 
result corresponds to theoretical expectations only when 
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the a-excitation function is associated with the pene- 
trability for ]=1, and the y function with /=0, and not 
the other way around, leading to the conclusion that the 
ground state of ;Li’ is odd. The penetrability has been 
taken to be proportional to 1/G*p in the notation of Yost, 
Wheeler, and Breit.* 


1 Gentner, Zeits. f. Physik 107, 354 (1937); Curran and Strothers, 
Proc. Roy. Soc. A172, 72 (1939). 

? Haworth and King, Phys. Rev. 54, 38 (1938); Rumbaugh, Roberts, 
and Hafstad, Phys. Rev. 54, 657 (1938). 

* Yost, Wheeler, and Breit, Terr. Mag. Atmos, Elec. 40, 443 (1935). 


16. The Isotopes of Cerium and Rhodium. ARNoLpD A. 
CoHEN,* University of Minnesota.—An analysis of the 
isotopes of cerium and rhodium has been made with a 
mass spectrometer in which ion currents are detected by 
an electron multiplier of the type described by Allen." 
Dempster’s* discovery of Ce and Ce™* was confirmed, 
and the abundance of these relative to Ce™® measured as 
1/500 and 1/400, respectively. No isotopes other than the 
four known ones were found, but the following upper limits 
of abundance relative to Ce™® were set: Ce™, Ce', 
1/20,000; Ce'#7, 1/10,000; Ce, 1/1000; Ce, 1/1500; 
Cel, Cel, Cel, Cel, 1/25,000. Rhodium was studied in 
an attempt to verify the existence of Rh™, reported by 
Sampson and Bleakney*® as being present to 1/1300 of 
Rh'®, Although 1/30,000 of Rh** could have been de- 
tected, no Rh™ peak was found. Other upper limits 
relative to Rh’ were determined as follows: Rh’, 
1/60,000; Rh’, 1/6000; Rh, 1/5000. Columbium was 
also investigated, but a zirconium impurity prevented a 
search for Cb™. An upper limit for Cb® was set at 1/4000 of 

* Present address: RCA Manufacturing Co., Harrison, New Jersey. 

'J. S. Allen, Phys. Rev. 55, 966 (1939). 


2A. J. Dempster, Phys. Rev. 49, 947 1936). 
3M. B. Sam W. Bleakney, Phys. Rev. 50, 732 (1936). 


17. Disintegration Scheme of Co®. L. G. ELLiott AND 
M. Deutscu, Massachusetts Institute of Technology.— 
Radioactive cobalt was separated from a nickel target 
which had been exposed to fast neutrons from the M.I.T. 
cyclotron. Examination of this source in the gamma-ray 
spectrometer! shows that it emits a gamma-ray of 0.810 
+0.015 Mev. No other gamma-rays were observed. The 
beta-ray spectrometer shows that it emits positrons with a 
maximum energy of 0.470+0.015 Mev. A counter sensi- 
tive to iron x-rays shows that a number of the disintegra- 
tions are by K capture. The number of observed coin- 
cidences between beta- and gamma-rays and a knowledge 
of the efficiency curve for the gamma-ray counter used 
show that each positron is on the average accompanied by 
one 0.810 Mev gamma-ray. Observation of coincidences 
between the x-rays and gamma-rays show that at least 80 
percent of the disintegrations by K capture are also to the 
excited state 0.810 Mev above the ground state in the 
Fe®* nucleus. 


1M, Deutsch and A. Roberts, Phys. Rev. 60, 362 (1941). 


18. The Mg II 24481 Line in the Spectrum of Algol. F. 
T. Rocers, Jr., The University of Houston ——A rather 
simple equation has been found for the profile of the \4481 
line in the uneclipsed spectrum of Algol. It is a representa- 
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tion of data taken by Struve and Elvey!' from seven spec- 
trograms (10A per mm dispersion at 44500), and was de- 
veloped for a further analysis* for rotation by the method 
of Carroll.* If \ represents wave-length in cm, A» being the 
proper wave-length of the line, then the measured absorp- 
tion in percent is ; 
i(x) = 25.8e71-2127, (1) 
where 
x= |A—Ao| X 108. (2) 


This representation differs from the actual data by an 
average deviation of 0.3 unit of absorption, which is only 
about one-third the A. D.’s of previous representations. 
10. Struve and C. T. Elvey, M. N. R. A. S. 91, 663 (1931). 
2 J. A. Carroll and Miss Ld. Ingram, M.N. R. A. S., 93, 508 (1933); 


also reference 1. 
3 J. A. Carroll, M. N. R. A. S., 93, 478 (1933). 


19. Thermal Diffusion Experiments with Ammonia and 
Argon. W. W. Watson AND D. WoeERNLEY, Yale Uni- 
versity—Calculations from viscosity data indicated that 
even at 450°K the thermal diffusion constant a of NH; 
might be negative. In view of the nature of the intermo- 
lecular forces for NHs, its critical temperature, and the 
known trend in Rr values with varying temperature 
(Rr=ratio of a to the value for hard spheres), it was 
thought that with. decreasing temperature Rr for NH; 
might change from + to — at some available T. With a 
“two-bulb” experiment, and NH; of about 15 percent N™ 
content to increase the accuracy of the mass spectrometer 
analyses, which were made after converting to Nz, we have 
detected this effect. The shift from.+ to — values for a 
occurs, however, at about room temperature. At mean 
temperatures of 377°K, 285°K, and 247°K, a@ is +0.011, 
—0.004, and —0.009, respectively, indicating a rapid de- 
crease in the repulsive force index and marked ‘‘softening”’ 
of the NH; molecule with decreasing 7. With tank argon 
at 100 cm*/hr, two 3-meter hot-wire columns in series have 
been set into operation, the procedure being based upon 
the results of R. Simon. With the wires at 1100°K the 
separation factor is observed to increase at the average rate 
of 2.56 per day. After four weeks no decrease in this steady 
rate of rise has been noted. 


20. Determination of Chlorine Vibrations in Dichloro- 
benzenes from Intensities and Polarizations of their Raman 
Lines. GERTRUD NORDHEIM AND H. Sponer, Duke Uni- 
versity.—To facilitate the assignment of observed Raman 
lines of the dichlorobenzenes to motions in which mainly 
chlorine atoms are partaking, the intensities and polariza- 
tion properties of these vibrations are studied theoretically. 
It is assumed that the change of polarizability in these 
vibrations is connected primarily with displacements of the 
chlorine atoms and that contributions from the individual 
atoms are additive. The Raman tensors belonging to the 
chlorine vibrations are then obtained by either adding or 
subtracting the single C—Cl contributions with reference 
to a chosen system of axes rigidly connected with the mole- 
cule. The following results were obtained: (1) Valence vibra- 
tions. The intensity ratio between the symmetrical and 
unsymmetrical vibrations is 8 : 1 (meta and ortho, the un- 
symmetrical vibration is forbidden in the para-spectrum). 


Lines representing symmetrical vibrations are highly polar. 
ized, lines representing unsymmetrical vibrations have the 
depolarization factor 6/7. (2) Deformation vibrations. Our 
simplified theory gives complete depolarization for all] 
allowed deformation vibrations, i.e., for symmetrical and 
unsymmetrical ones and also in and out of the plane. Lines 
from planar symmetrical vibrations are three times stronger 
than from unsymmetrical ones. Similar relations hold for 
nonplanar vibrations. These results are discussed in con. 
nection with existing experimental data. 


21. Transition Probabilities of the Asymmetric Rotor, 
Paut C. Cross, Brown University anD G. W. Kine, 
Arthur D. Little, Inc-—The successive convergents of the 
continued fraction used in calculating the energy levels of 
the asymmetric-rotor are closely related to the elements of 
the transformation matrix expressing the asymmetric-rotor 
wave-functions as linear combinations of symmetric-rotor 
functions.! This transformation was evaluated and applied 
to the direction-cosine matrices for the symmetric rotor 
in order to determine the Einstein transition probabilities 
for rotors having various degrees of asymmetry. The values 
determined enable interpolation to be carried out with 
sufficient accuracy for any molecule. The change in in- 
tensity of a given transition with increasing asymmetry is 
most marked for levels of high or low r-values, i.e., where 
the difference in the two limiting K values is highest. 
Strictly speaking, the intensities given apply only to the 
pure rotational spectrum of a rigid molecule. In electronic 
and vibrational transitions the moments of inertia are 
markedly changed, but the tables can be used to give a 
satisfactory picture of the rotational structure to be 
expected in any given case. 


ess W. King, R. M. Hainer and P. C. Cross, J. Chem. Phys. 11, 27 


22. Raman Spectra of Symmetrical Top Acetylenes. 
Forrest F. CLEVELAND AND M. J. Murray, Illinois Insti- 
tute of Technology—Raman data have been obtained for 
(A) 1-iodo-1-propyne (CH;C =ClI), (B) 1-bromo-1-propyne 
(CH;C =CBr), (C) 1-chloro-1-propyne (CH;C =CCl), and 
(D) 3,3-dimethyl-1-butyne (C(CH;);C =CH) in the liquid 
state. Relative intensities and depolarization factors were 
obtained by use of a Gaertner microdensitometer. C2H;Br 
was present as an impurity in the sample used in obtaining 
the spectrum of (C) and depolarization factors were not 
obtained for this compound. Calculation shows that the 
CX;C=CY molecule should have five A; frequencies and 
five E frequencies, all of which are active in the Raman 
effect. The A; frequencies should be polarized, the E fre- 
quencies depolarized. The observations confirm these ex- 
pectations, although there are some cases of doubling of 
lines due to resonance interaction. Since there is only a 
small increase in frequency for the gaseous state for such 
molecules, it is possible to make approximate thermo- 
dynamical calculations for molecules (A) and (B). In the 
triple bond region a line of medium intensity and smaller 
displacement than the strong triple bond line was observed 
for compounds (A), (B), and (C). This may be due to 
resonance interaction with the difference tone of the 2900 
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and 700 cm~! fundamentals. For compound (D), two lines, 
rather than the usual one for monosubstituted acetylenes. 
were observed at 2106 and 2137 cm™. 


23. A Dynamical Approach to the Condensation Prob- 
lem. M. DRESDEN AND WILFRED Kaptan. University of 
Michigan.—A number of attempts have been made to 
explain the condensation of a gas on the basis of sta- 
tistical mechanics, starting from the classical or quantum- 
mechanical partition function. From this one proved essen- 
tially the existence of a singularity in the ~, V curve for 
sufficiently low temperatures. However this proof does not 
elucidate the physical mechanism of condensation. It might 
be possible to understand the condensation phenomena by 
linking it with the stability conditions for dynamical sys- 
tems. Consider the gas as a classical N-body problem, then 
by generalizing the concept of dynamic stability we can 
show that the stability of the system is essentially deter- 
mined by the local structure of the energy manifold in the 
3N dimensional phase space. So in order to determine where 
transitions can occur we must find out where the topology 
of the energy manifold changes, as we vary the total energy. 
For simple cases it has been shown that such changes 
occur;! dynamically, they correspond to the transition 
from bounded to unbounded orbits. The extension to the 
general case and the possibility of computing transition 
temperatures will be discussed. 


1 Wilfred Kaplan, Am. Math. Monthly, No. 5 (May, 1942). 


24. The “Génies’’ of Poincaré, and Carnot’s Principle. 
PrerreE Demers, Montreal.—Poincaré has proposed that 
perhaps a being that he calls a ‘“‘génie,’’ would see the in- 
dividual molecules in a fluid mass; he could manipulate a 
trap to allow the faster molecules to accumulate in one 
compartment, where temperature would rise, against Car- 
not’s principle. Maxwell has proposed the same, with 
“demons.”” This scheme was supposed to involve no 
energy expense. The theoretical possibility of this experi- 
ment is examined with the following conclusions: (1) Ma- 
nipulating the trap involves a minimum energy ~kT/2, 
from kinetic considerations. (2) To locate the molecules we 
can detect photons, diffused or emitted in a thin layer of 
molecules. A microscopist, or an automatic device, could 
manipulate the trap. It is seen that anyway some energy 
must be spent from the outside. (3) We can enclose the 
molecules to avoid loss of radiation, keeping them in 
thermal and radiative equilibrium; then we have a black- 
body, in which details and individual molecules are indis- 
tinguishable; the ‘‘génie” would be unable to locate the 
origin of any one quantum perceived. Not seeing the mole- 
cules, he would be unable to sort them. 


25. On Strong-Coupling and Weak-Coupling Theories 
of the Meson Field. W. Pau, The Institute for Advanced 


Study —With the assumption that the meson has spin 0 
and is described by a pseudo-scalar field, the model of an 
extended heavy particle (nucleon) leads to the strong 
coupling condition f*>>a*, where f is the coupling constant 
with the dimension of a length, and a is the size of the 
source. In connection with the prediction of excited states 
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(isobars) of the nucleon, this theory leads to serious diffi- 
culties, since Dancoff and Serber' showed that the theory 
is in contradiction to the well-known saturation properties 
of the nuclear force, and Pauli and Dancoff? showed that 
it gives for the magnetic moments of the proton and the 
neutron values equal in magnitude but opposite in sign. 
On the other hand purely classical theories of the meson 
field which are based on the picture of a point source and 
which fulfil the condition of relativistic invariance were 
developed by Bhabha? and Fierz.* The singularities of the 
field due to the point character of the source are avoided 
by a suitable subtraction formalism. The author has shown 
that the classical point-source theories can be translated 
to quantum theory by the use of the so-called \ limiting 
process of Wentzel and Dirac. The application of this 
method to the theory of meson field produces automatically 
the value 0 for the spin inertia, K, and gives convergent 
results for all quantum theoretical problems where the nu- 
cleons can be considered to be at rest. Such a weak coupling 
theory with point sources avoids all difficulties of the strong 
coupling theory with an extended source. For the scattering 
cross section of mesons by nucleons (due to the non- 
electromagnetic interaction), the results of this theory 
agree completely with that given by Heitler and Peng,* 
which is in fairly good agreement with experiment. No 
stable isobars of the nucleon exist in such a theory as long 
as (fx)*X1, a condition which is certainly fulfilled. In the 
problem of the interaction between two (or several) nu- 
cleons at rest, the perturbation theory holds for distances, 
r, of the nucleons for which r*>>f? which can be derived by 
the method given by Stiickelberg and Patry.* 

1S. M. Dancoff and R. Serber, Phys. Rev. 61, 394A (1942). 

?W. Pauli and S. M. Dancoff, Phys. Rev. 62, 85 (1942). 

37H. J. Bhabha, Proc. Roy. Soc. A178, 314 (1941). 

4M. Fierz, Helv. Phys. 257 

5 W. Heitler and H. W. P. Proc. Camb. Phil. Soc. 38, 296 (1942). 

Helv. Phys. Acta 13, 167 


*E.C G. Stiickelber = . F. C. Patry, 
(1940). E. C. G. Stiic jelv. Phys. Acta 13, 347 (1940). 


26. On the Theory of a Mixed Pseudo-Scalar and a 
Vector Meson Field. S. KusAKA AND W. Pau, The Insti- 
titute for Advanced Study.—The authors have investigated 
in more detail Schwinger’s' modification of Mgller-Rosen- 
feld’s theory of nuclear forces which is based on a mixture 
of a pseudo-scalar and a vector meson field. With the as- 
sumption of an extended source of the nucleon, the strong 
coupling condition is fulfilled, and if the coupling constants 
of the two fields are taken to be equal, the excitation energy 
of the isobars turns out to be } of its value in the pseudo- 
scalar theory. The interaction between two nucleons is 
described with sufficient accuracy by a Hamiltonian which 
consists of three parts: (1) the isobar excitation energy of 
the nucleon; (2) the kinetic energy of the nucleon; and (3) 
the potential energy of the interaction between the nu- 
cleons. In the most interesting case of the symmetrical 
theory, the last expression is the same as the corresponding 
expression of the weak coupling theory with the spin and 
the isotopic spin matrices replaced by the more general 
ex* vectors of Pauli and Dancoff. A classification of the 
energy levels of the two-nucleon system has been made, and 
a preliminary estimate indicates that a suitable choice of 
the coupling constant and the masses of the mesons can 
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be made in such a way that the binding energy and the 
quadrupole moment of the deuteron agree with the observed 
values, and such that the size of the source is small com- 
pared to the Compton wave-lengths of both mesons and 
that the excitation energy of the isobar is large compared 
to the binding energy of the deuteron. However, this theory 
gives for the magnetic moment of the deuteron a value only 
a few percent of the observed value, and the difficulties for 
higher nuclei described in the preceding paper are also 
present here. Hence we propose to consider this theory in 
the weak coupling scheme with a point source of the kind 
described in the preceding paper. More accurate calcula- 
tions for the wave functions of the ground state of the 
deuteron are in progress. 


1 J. Schwinger, Phys. Rev. 61, 387A (1942). 


27. On the Meson Field Theory of the Magnetic Mo- 
ment of Proton and Neutron. J. M. Jaucn, Princeton Uni- 
versity.—It is shown that the application of the \-limiting 
process together with Dirac’s introduction’ of a second 
type of field oscillators with negative energy leads in the 
case of the pseudo-scalar theory to a convergent result for 
the magnetic moment of the heavy particles in the approxi- 
mation of the order of the square of the coupling constant. 


1P, A. M. Dirac, Proc. Roy. Soc. A180, 1 (1942). 


28. Ionization Yield of Radiation Energy. U. FANo, 
Carnegie Institution of Washington.—According to experi- 
ments with ionizing radiations acting on various sub- 
stances, the average energy absorbed per ionization pro- 
duced is approximately 30-35 ev and bears no apparent 
relation to the ionization potential of the particular sub- 
stance. Bagge’s theoretical calculation of this energy proved 
correct for atomic hydrogen, but for other substances in- 
correctly indicated a positive correlation with the ioniza- 
tion potential.'! This discrepancy is accounted for qualita- 
tively by giving proper consideration to the fact that most 
energy absorptions are due to electric dipole oscillations. 
Available theoretical information shows that the oscillator 
strength of ionization in certain characteristic atoms is 
directly and critically related to their ionization potential: 
the greater the energy actually required to ionize individual 
atoms, the smaller is the average energy ‘‘wasted” in pro- 
ducing excitations. A crude theory complementary to 
Bagge’s has been applied to H and He. It indicates that 
the factor considered here is essentially sufficient to balance 
the great difference between the ionization potentials of 
these two atoms. This accounts for one particularly puzzling 
feature of the known experimental data, but does not 
positively establish a yield of one ionization per every 30-35 
ev as a theoretical result valid for all substances. 


1 E. Bagge, Ann. der Physik 30, 72 (1937). 


29. A High Pressure Cloud Chamber of New Design. 
T. H. Jounson, R. P. SHutt, AND S. DE BENEDETTI, Bartol 
Research Foundation.—A cloud chamber is described which 
operates at a pressure of 200 atmospheres without the 
difficulty of having to subject its transparent walls to high 
pressure differentials. The chamber consists of a glass ring 


of }-inch wall thickness, 12 inches in diameter, and 3 inches 
deep covered on the front with a }-inch thick glass plate of 
the same diameter and on the back with a synthetic rubber 
diaphragm. The chamber is mounted in the center of a 
large steel tank. The tank is filled with transparent white 
oil which subjects the outside of the cloud chamber to a 
hydrostatic pressure equal to the internal gas pressure. The 
chamber is compressed by pumping additional oil into the 
tank, and it is expanded by releasing a valve connecting 
the tank with an evacuated vessel of adjusted volume, 
Ample space is provided inside the tank for magnetic field 
coils and the heavy walls of the tank will serve as the mag. 
netic conductors. The apparatus has already been operated 
at a pressure of 110 atmospheres giving excellent cosmic- 
ray tracks. At this pressure its efficiency for observation of 
rare events such as mesotron decay is about 1000 times 
that of a chamber of the same size operating at atmospheric 


pressure. 


30. Some Results Obtained with a New High Pressure 
Cloud Chamber. S. DE BENeEpetT!I, T. H. JOHNSON, AND 
R. P. Sautt, Bartol Research Foundation.—460 pictures 
have been taken in a newly designed high pressure cloud 
chamber, most of them in an atmosphere of n-propyl 
alcohol, water vapors, and argon at a pressure of 110 atmos. 
pheres. The volume expansion ratio decreases almost 
linearly with pressure, and at 110 atmospheres it becomes 
as low as 1.05. Owing to small diffusion at such pressures, 
the sensitive time of the chamber is increased by almost 10 
times that at normal pressure, making counter control 
unnecessary. At 110 atmospheres each random picture 
contains the equivalent of 3 full length tracks, correspond- 
ing to more than 10‘ cm of air at normal pressure. A total 
length of 510* cm of air equivalent has been observed, 
which equals the track length obtained in 500,000 random 
pictures, or 150,000 pictures controlled by counters, in a 
cloud chamber of the same dimensions but at normal pres- 
sure. Among the results obtained are several cases of 
heavily ionized tracks attributable to the proton compo 
nent, one case of a mesotron disintegration, and several 
pictures of “‘knock on” electrons of energies above 20 Mey 
produced in the gas. Typical pictures will be shown. 


31. Experiments Concerning the Origin of the Second 
Hump in the Rossi Curve. W. F. G. Swann, Bartol Research 
Foundation of the Franklin Institute —Former experiments 
at the Bartol Foundation have confirmed the existence of 
the second hump. These experiments have now been varied 
in a manner calculated to give further information on the 
significance of the hump. The minimum between the two 
humps occurs at a thickness of about 10 cm of lead. After 
carrying observations out to this thickness, the thickness 
was further increased by the addition of tin. This caused an 
immediate rise in the curve of an amount comparable with 
that for the first hump in tin itself. On the other hand, 
when a similar experiment was performed starting with a 
thickness of tin corresponding to the minimum and it 
creasing the thickness by the addition of lead, a diminution 
in the intensity of shower production with increase of 
thickness was observed. These experiments were suggested 
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by an attempt to account for the second hump on the 
hypothesis of the creation of electrons by the primary ra- 
diation after a certain degree of penetration, and their 
significance is discussed in relation to this view. 


32. Production of Mesotrons by Ionizing and Non- 
Ionizing Radiation. MARCEL SCHEIN AND JULIUS TABIN, 
University of Chicago.—A cosmic-ray balloon experiment 
was carried out in the stratosphere using Geiger-Mueller 
counters in 6 different 4- and 5-fold coincidence arrange- 
ments. The equipment reached an altitude corresponding 
to 2.5 cm of Hg pressure. The following phenomena were 
investigated: (a) single mesotrons produced in 2 cm of Pb 
by non-ionizing rays, (b) multiple production of mesotrons 
by ionizing and non-ionizing rays, (c) cascade showers 
below 2 cm of Pb produced by electrons and photons, (d) 
high energy electron showers. Each individual event was 
registered separately on the same film and it was possible 
to obtain all the correlations between the different events. 
The results obtained can be summarized as follows. At a 
pressure of 3.cm Hg (first radiation unit from the top of 
the atmosphere) 47 percent of the penetrating particles is 
accompanied by multiples (mesotron showers) produced by 
ionizing rays in 7 cm of lead. The producing radiation 
probably consists of primary protons. At altitudes higher 
than 14 km the number of mesotrons produced by non- 
ionizing rays in a lead thickness of 2 cm is as high as 28 
percent of the intensity of the penetrating component. The 
mesotrons produced in this process are of relatively low 
energies (1 to 2X 10® ev) and most of them (87 percent) 
occur in singles. The other 13 percent probably consists of 
mesotron pairs or mesotron-proton pairs. The non-ionizing 
rays producing these mesotrons appear to be low energy 
photons. Practically none of these processes was accom- 
panied by cascade showers registered below 2 cm of Pb. 


33. X-Ray Diffraction Investigation of Sodium Stearate 
from Room Temperature to the Melting Point. ALEXANDER 
DE BRETTEVILLE, JR. AND JAMES W. McBain, Stanford 
University—X-ray photographs were obtained of the 
y-form! of sodium stearate by means of a heating camera* 
starting at room temperature and continuing through the 
melting point at 288°C to 301°C. The photographs show a 
gradual coming together of the two most intense side 
spacings to form a single broad line, somewhere between a 
temperature greater than 145°C and less than 170°C. Since 
a phase change’ exists at 167°C it is assumed that this is the 
temperature at which the radical change in structure takes 
place. Orientation of the long spacings is exhibited at 301°C. 
Complete interpretation of the photographs is not possible 
because of a lack of crystal structure data. However, as- 
sumptions can be made about the structure from other 
sources that enable one to show that sodium stearate 
exhibits phenomena analogous to phenomena exhibited by 
the paraffins.‘ 
10h, de Bretteville, Jr. and J. W. McBain, Science, 96, 470 (Nov. 20, 

* A. de Bretteville, Jr., Rev. Sci. Inst., 13, 481 (1942). 


*M. J. Vold, J. Am. Chem. Soc. 63, 163 (1941). 
‘A. Miiller, Proc. Roy. Soc. 138, 514-530 (1932). 


34. The Shape of the Debye-Scherrer Lines. C. C. 
Murpock, Cornell University.—The equation for the shape 
of the Debye-Scherrer line has been derived for crystals of 
uniform size for the following cases: the (hk1) line for 
tetrahedral and cubic crystals, the (hk k 0) and (h & k) lines 
for octahedral crystals and the (100), (110), and (111) lines 
for rhombic-dodekahedral crystals. When the crystal has a 
center of symmetry, the lines associated with Bragg planes 
which are approximately parallel to crystal faces have 
secondary maxima. Because of these the estimation of line 
breadth from experimental data in terms of “integral 
breadth” is less reliable than in terms of “half intensity 
breadth,” By. Let K =(Byv! cos @)/d in which » is crystal 
volume. For the dodekahedron, (100) and (110), K =0.87; 
(111), K=0.91. For the octahedron, K ranges from 0.82, 
(100), to 0.94, (111). For the cube, from 0.83, (110), to 0.91, 
(311), (411), and (511). For the tetrahedron, from 0.73, 
(110), to 1.03, (100). The equations give the same values 
for “integral breadth” as those computed by Stokes and 
Wilson.' For measurement of particle size, a line should be 
used for which X is insensitive to crystal shape. The results 
so far obtained suggest the (210) reflection, K =0.89+0.02, 
for this purpose. 

1 Proc. Camb. Phil. Soc, 38, 313 (1942). 


35. An X-Ray Diffraction Study of the Age-Hardening 
of a Cu-Al Alloy. CHartes H. TowNEs AND WHEELER P. 
Davey, The Pennsylvania State College—The age-harden- 
ing of Cu—Al alloys has been studied by many investiga- 
tors using density, elastic limit, elongation, hardness, 
electrical conductivity, and lattice parameters. Results are 
summarized elsewhere.! The present work (using balanced 
filters and Geiger-Mueller counters) shows the time rate 
of change of the diffracting powers of the atomic planes of 
the Cu—Al solid solution and of its decomposition products 
(6’ and @). Pre-strained 25S alloy was heated for 24 hours 
at 540°C (to dissolve all the Cu in the Al) and then water- 
quenched at 20°C. The Cu was then allowed to migrate at 
309°C. It was found that: (1) In the Al “solid solution” 
phase the diffracting power of the (111) planes decreased to 
a constant value in 120 hours; the (100), (110), and (311) 
planes showed a decreased diffracting power followed by a 
considerable increase. (2) The (112) (2) plane of 6’ de- 
creased to a constant diffracting power in 120 hours. (3) 
The (001) (3), (101) (2), (112), (110) (2), and (211) planes 
of @ increased in diffracting power to a constant value in 
approximately 120 hours. Definite interpretations of the 
above are offered. 


“Age Hardening of Metals,” Am. Soc. for Metals (1940). 


36. X-Ray Measurement of Short Range Order and 
Long Range Order in Cu;Au. Z. W. Wiccuinsky, Massa- 
chusetts Institute of Technology—From a powder pattern, 
the value of long range order, S, can be obtained from 
intensity measurements of the superstructure lines. When 
S is zero it is possible to obtain the value of short range 
order by applying Fourier integral analysis to intensity 
measurements of the diffuse background. From diffraction 
theory the background intensity for the alloy A;B is 
proportional to (fa—fs)? (3— 2; sin kri/kri) where fa 
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and fg are the atomic scattering factors, k= 4m sin the 
r;’s are the interatomic distances, c; is the number of atoms 
at distance 7; from a given atom, and a; is the generalized 
form of the short range order parameter as defined by 
Bethe. The dominant parameter, o:, was found to be 0.19 
for sample quenched from 750°C; 0.55 for 540°C; and 0.67 
for 400°C. The intensity of a superstructure line is pro- 
portional to S*(fa4—fz)*. Values of S obtained are 0.842 for 
sample held at 250°C; 0.814 for 370°C and 0.804 for 380°C. 


37. A Simplified Electron Microscope. C. H. BACHMAN 
AND SIMON Ramo, General Electric Company.—An electron 
microscope has been constructed to test a number of new 
attacks. The aim has been to obtain a resolving power in 
the range of the greatest present usefulness of electron 
microscopy (say ten times better than light microscopy), 
operating at a voltage range sufficient to obtain this resolu- 
tion for the usual thin specimens, in as simple a manner as 
possible. The instrument is a combined electron and light 
microscope and utilizes unipotential electrostatic electron 
lenses. It has the following desirable characteristics: (1) 
The electrons are accelerated and focused by a single 
crudely regulated voltage. (2) The unipotential lenses can 
readily be used in a multiple-stage magnifying system of 
compact construction so that the chamber to be evacuated 
becomes very small. (3) A short electron path results from 
the use of a multiple lens electron-optical system in combi- 
nation with a stage of light-optical magnification. This 
small cathode to fluorescent screen distance and the me- 
dium resolving power aim make it possible to align the 
entire electron-optical system permanently. (4) The cham- 
ber can be re-evacuated to operating conditions within a 
few minutes after insertion of a new specimen. (5) The over- 
all weight and dimensions, if proper advantage is taken of 
the above conditions, are such that it is feasible to design 
a portable instrument. 


38. Compact Electron Microscope with Magnetic Lenses. 
V. K. ZworyKIn AND J. HILLIER, RCA Manufacturing 
Company.—To meet the demand of a small, simply oper- 
ated and inexpensive electron microscope, a desk-type 
instrument has been developed, with a total length from 
electron source to observation screen or photographic plate 
of only 16 inches. The objective and projector lenses are of 
a new design, and the magnetomotive force is provided 
either by a coil or by permanent magnets. Focusing is ac- 
complished by varying the stabilized operating voltage, 
normally maintained in the neighborhood of 30 kilovolts. 
The illumination is controlled by a bias voltage applied to 
a grid cylinder surrounding the cathode. The total volume 
of the instrument being only about 1 liter, airlocks are 
superfluous; the pumping-out time required for complete 
evacuation is only about 2 minutes. The microscope is 
mounted so that the observer views the final image directly 
on the rear of a transparent fluorescent screen. 


39. Distance to Which the Attractive Field of a Solid 
Extends into a Liquid. Wi_L1aM D. HARKINS AND GEORGE 
Jura, University of Chicago and the Universal Oil Products 
Company, Chicago.—Various estimates of the distance to 
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which the attractive field of a crystal extends into an 
adjacent liquid vary from the thickness of one to several 
thousand monolayers. A method has been developed for the 
determination of the variation of the energy with the dis. 
tance. A definite amount of the vapor of the liquid js 
adsorbed on the surface of the solid, which must give a 
zero contact angle with the liquid. The powder is then 
immersed in the liquid in a sensitive calorimeter. The de. 
termination of area is described in the following paper, 
From the integral curve the energy of desorption per 
molecular layer, is calculated. With titanium dioxide 
(anatase) and water the excess of the energy of desorption 
over that for vaporization (10,480 cal. mole) for the suc. 
cessive molecular layers is; 5360, 480, 110, 60, and 20 cal, 
mole~!. Quartz would give slightly higher values. Although 
the energies of attachment of the higher layers are very 
small, films of slightly over nine layers of water or of 
nitrogen have been obtained before saturation is attained, 
The accuracy of the values can be improved by more ex. 
tensive experimental work. 


40. A Method for Determining the Area of a Crystalline 
Powder without the Assumption of a Molecular Area, 
GEORGE JURA AND WILLIAM D. HarKINs, University of 
Chicago and Universal Oil Products Company.—Emmett, 
Brunauer, and Teller have developed a gas adsorption 
method which gives the volume of gas adsorbed to fill 
the first monolayer on the surface of a solid. If the molecu- 
lar area were known for each case the area could be calcu- 
lated. In the application of this method they used porous 
solids of extremely large area, while we employed non- 
porous crystalline powders. We have developed a new 
method, independent of any assumption of a molecular 
area, based on a proof that a solid surface which gives a 
zero contact angle with a liquid is covered by a duplex 
film of that liquid when in equilibrium with saturated 
vapor. A duplex film is one just thick enough to exhibit 
the same surface energy as the liquid itself. At 250°C 
fine crystals of titanium dioxide, in equilibrium with 
saturated water vapor, and therefore covered by a duplex 
water film, were immersed in liquid water. Since the surface 
energy was that of water (118.5 erg cm), the total energy 
of immersion, divided by this value, gave the area (140 
meter g~'), while nitrogen adsorption gave exactly the same 
area when Emmett’s molecular area of 16.2A? was assumed. 


41. Raman Spectra of Methylheptenes. Forrest F. 
CLEVELAND, Illinois Institute of Technology.*—Continuing 
work on olefinic hydrocarbons containing eight carbon 
atoms, the Raman frequencies, relative intensities and 
depolarization factors of 2-methyl-1-heptene and 6-methyl- 
1-heptene have been obtained. The relative intensities and 
depolarization factors were obtained by use of a Gaertner 
microdensitometer. The strong lines characteristic of the 
RHC=CHsz group, previously observed at 1293(94)0.88, 
1415(42)0.77, 1641(134)0.15 and 3078(23)0.93 for 1-octene, 
appear at 1295(212)0.45, 1417(112)0.56, 1644(342)0.19 and 
3079(70)D for 6-methyl-1-heptene. For 2-methyl-1-heptene, 
in which the group is R(CH:)C=CHz, the corresponding 
lines were observed at 1303(52)0.64, 1414(108)0.81, 
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1653(203)0.20 and 3077(47)D. It may be noted that the 
1653 frequency is higher for the R(CH:)C=CH: group 
than the corresponding value of 1642 in the RHC=CH; 
group. Furthermore, the frequency which was constant 
at 3000 cm~ for the four octenes previously studied and 
for the 6-methyl-1-heptene drops to 2978 cm™ for the 
2-methyl-1-heptene. The line 1334(48)0.55 in the spectrum 
of 6-methyl-1-heptene was not observed for the 2-methyl-1- 
heptene or for the four octenes previously studied while the 
line observed at 1383(69)0.64 for the latter compound was 
observed in the previous investigation only for cis+trans 
2-octene and trans-3-octene. 


* To be called for after Paper No. 22 if time permits. 


42. A Method for the Accurate Determination of the 
Adsorption of a Vapor on the Surface of a Solid. WiLL1am 
D. HARKINS AND GEORGE JuRA, University of Chicago.*— 
While it is easy to obtain a moderately accurate isotherm 
for the adsorption of a vapor on a porous solid by the 
ordinary volumetric method, or by a McBain-Bakr quartz 
spiral, these methods are not accurate if non-porous sheets 
of powders are under investigation. Our method is ex- 
tremely accurate, not affected by time, and as many 
temperatures as are desired may be used. (1) As large a 
sample of sheets or powder as possible is put in a glass 
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bulb in which it can be outgassed. During the adsorption 
process this bulb is kept in a thermostat. (2) The samples of 
liquid are placed in sealed glass bulbs filled with proper 
amounts of the liquid. The empty and full bulbs are weighed 
on an analytical or a microbalance, and are held in a 
horizontal tube which is attached by a vertical tube to the 
adsorption bulb, and also to a small sensitive manometer. 
The horizontal tube is supplied with a mechanism to break 
any single bulb as desired. (3) In the application of this 
method it is often possible to adsorb on the solid 99 percent 
of the mass of liquid in the bulbs. The manometer is used 
to obtain the vapor pressure in order that the mass of 
vapor unadsorbed may be determined. 


* To be called for after Paper No. 40 if time permits. 


43. The Lift Function for an Oscillating Airfoil.* W. C. 
RANDELS AND S. SiLverR, University of Oklahoma.—A 
method used by Biot! is applied to an airfoil undergoing 
torsional oscillations. The usual approximations of the 
linearized theory for thin airfoils have been made. The 
boundary conditions are shown to be satisfied by a function 
consisting of the sum of three functions corresponding to 
elementary singularities. The total lift and moment are 
given. 


* To be read by title. 
1M. Biot, J. Aer. Sci. 9, 185 (1942). 
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